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A USEFUL DEVELOPMENT IN THE CHEMISTRY OF 
FLAVONOIDS AS NATURAL PRODUCTS 


by T. R. SEsnapnr, F.R.S., E.N.I., President, National Institute of 
Sciences of India 


(Delivered on 2 January 1969) 


IMPORTANCE Or Ginkgo biloba (MAIDEN-HAIR TREE) 


Ginkgo biloba (maiden-hair tree) is a tree of great botanical and chemical 
interest (Fig. 1). It is the only surviving species among Ginkgoales which 
were widespread in the very distant past between 60 and 250 million years 
ago. It has apparently remained unchanged over geological ages and it 
has therefore been called a ‘living fossil’. It has long been cultivated in 
temple gardens of Japan where it is considered sacred. Not only is the tree 
beautiful but the green fruits are eaten after frying. The flesh of the ripe 
fruits is found to contain toxic phenols. The tree has been introduced into 
European countries and the U.S.A. 

The extract of the leaves of this tree (Fig. 2) has been reported to have 
medicinal properties and it is now being marketed in West Germany in forms 
suitable for oral administration and injection. It is considered to possess 
blood flow increasing and vascular widening properties and is recommended for 
use in hormonal, nervous, vascular, blood flow and nutritional disorders. 

The special chemical components of these leaves have been studied 
intensively during the past ten years and their structures have been established. 
They are known as ‘biflavonoids’ in whose molecules two Ci, units are linked 
together. In view of the important medicinal properties attributed to the 
leaves we have recently examined Indian plant sources for compounds of 
this type and have obtained some useful results. 

Before proceeding further, it may be useful to mention earlier discoveries 
on the medicinal uses of flavonoids. These polyphenols were originally 
valued as mordant dyestuff and this use has become obsolete after the intro- 
duction of synthetic dyestuffs. But their importance has continued for other 
reasons. One of these is their vitamin P property. "This was discovered by 
Professor Szent Georgyi more than 30 years ago. At the present time Rutin, 
which is obtained on a large scale from buck-wheat and eucalyptus leaves, is 
used in medicine. Similarly the bioflavonoids, hesperidin and naringin ob- 
tained as by-products of the citrus industry are employed for the maintenance 
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of capillar y integrity, control of spontaneous abortion and to aid 1n su gery 
and W ound healing. 


Finkgetin and related compounds 

A. leaves of G. biloba were firs& examined chemically by Furukawa a 
1932. He proposed the structure of 5, 8-dihydroxy-4’-methoxy A 
methoxy primetin) for the chief component which he labelled as l P : es 
B’. This was later on named ‘ ginkgetin’ by Baker who also showed t 2 i 
was different from synthetic—4’-methoxy primetin. ‘The study Was progr rl 
ing for over 25 years without conclusive results. The main difficulties were 


Cd z 


Fig. 1. Ginkgo biloba. Fia. 2. Ginkgo biloba—Leaves, 


as follows: (i) the leaves yielded a mixture of three closely related compounds 
and their separation was originally difficult; (ii) much of the earlier study 
was carried out with mixtures; (iii) there was also difficulty of correctly deter- 
mining the molecular size of these compounds. Eventually it was concluded 
that they were biflavonoids containing two flavonyl units and their consti- 
tutions were arrived at by the very extensive work of Japanese and Br 
workers. 

The crude material contained at least three related n 
which two major ones were ginkgetin and isoginkgetin. 


itish 


atural products of 
It was later shown 
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that these compounds could be comparatively easily isolated by taking ad- 
vantage of Nakazawa's observation that ginkgetin formed a rather insoluble 
erystalline potassium salt. Treatment of the mixture with aqueous potas- 
sium carbonate gave the sparingly soluble potassium salt of ginkgetin and 
acidification of the mother liquor gave isoginkgetin. Later the surprising 
observation was made that the mixture could be easily separated by fractional 
crystallization from acetone. Since by the use of the Rast method for the 
determination of molecular weight satisfactory results were not obtained, 
the modified Menzies- Wright ebullioscopic method was used and the molecular 
formula of ginkgetin was found to be CggH2201 and it contained two methoxyl 
and four hydroxyl groups. 

The UV spectra of ginkgetin tetraacetate and tetramethyl ether indi- 
cated that they are derivatives of apigenin. Moreover, the intensities of the 
absorption bands were approximately double those of the monomeric apigenin 
derivatives. This demonstrated the presence of two isolated flavonoid units 
per molecule. Further studies of spectra particularly IR and colour reactions 
gave more information about the location of groups. But definite informa- 
tion was obtained by chemical degradations. The most important of these 
are given in Chart 1. 


H20 
A o OCH, + TOR 
CINKGETIN 


COOH 
je KOH 9 


HCO. OH H¿cOC Y OCH 
3 
Ho) coc so qs 3 
HO o OQ OH 
OH | 


OH O 
HOOC « 2 OCH, 
ng NC H¿CO OH 
CINKCE TIN [0cH4) H202, H,co OH 
3% 3 ur COOH 


COOH H4CO 
HCO 


Cuart 1. Chemical degradations of Ginkgetin. 


They led to the conclusion that the flavonyl units had C-C linkage and 
were linked at 3’ and 8”-positions. The position of the methoxyls was 
not only indicated by the nature of the fission products but also confirmed by 
the study of UV spectra in the presence of bases. The constitution proposed 
by Baker and co-workers in the year 1959 was 4', 7-O-dimethyl-3', 8"-bi- 
apigeninyl (Chart 1a). 

Nakazawa effected a synthesis of ginkgetin tetramethyl ether in 1959 by a 
mixed Ullmann reaction between 3'-iodo-5, 7, 4'-trimethoxyflavone and 
8-iodo-5, 7, 4'-trimethoxyflavone using copper in dimethyl formamide medium. 
One of the products was found to be identical with the tetramethyl ether of 


3 
Ai 
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natural ginkgetin n akaz y i i in 1 lows : 
atural gink ti In: 1963, N kazawa S nthesized ginkgetin itself as fol 
i n. Y c j I i à f 
He used for the n of the free hydroxyl groups in the is Ms 
b 1 and benzo 1 units and th se were conveni n l R ed later o 10) 
enzy. y e V ent y remov 1 t n $ 


yield ginkgetin (Chart 2). 


HO 0 


CHART la. Ginkgetin. 


As already mentioned two more related compounds (partial methyl 
ethers) have been isolated from G. biloba leaves. They are isoginkgetin 


+ ULLMANN y 


[Naxazawa - 1959) 


OCH,” 
ey Qr 
i 
0 H3CO 
OR | 
+ Il. ULLMANN | 


2.R—H 


[ NAKAZAWA — 1963] 


OR 2. R=COC6Hg . 


CuarT 2, Synthesis of Ginkgetin and its methyl ether. 


and bilobetin. The former is an isomeric dimethyl ether and the 
latter a related monomethyl ether. Soon many other methyl ethers of this 
group were isolated. The parent compound of all these ethers is called 
“amentoflavone”. Originally it was considered to be rare; now several sources 
have been found. The trimethyl ether “sciadopitysin” is of importance; it 
was isolated by Kariyone and Kawano in 1956 and their extensive work on it 


was complementary to the work of Baker et al. (1959) on ginkgetin (Chart 3). 
IB 
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AMENTOFLAVONE 
BitogeriN 


PovocarPusrLavone A 


SoTETSUFLAVONE 


Ginkcetin 


Tsocinkee TIN 


Povoc ARPUSFLAVONE B 


SciaooPirYsiN 
Kavaf Lavone 


ÁMENTOFLAVONE 
TETRAMETHYL ETHER 


CHART 3. Amentoflavone series. 
Hinokiflavone 


A different type of biflavonoid was isolated by Kariyone and Sawada in 
1958 from Chamaecyparis obtusa. It was named ‘hinokiflavone’ based on the 


Japanese name of the plant “hinoki” (Chart 4). The constitution of this 


Sa 


OH O 


CHART 4. Hinokiflavone. 


O 


compound was arrived at by Kawano in 1960. It was also a dimer of api- 
genin, the molecular formula being C3,H,50;,; it had no methoxyl groups. 
It formed a pentamethyl ether and pentaacetate. The UV pattern was 
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i i ls but the UV of the complete 

i f the earlier known biflavony. ; pee 
E er E with the UV of ginkgetin E ET E 
B T cyls were inferred to be 

'eacti IR data two hydroxyls : 

E i Nonus degradations of hinokiflavone and its pentamethyl ether 

ositions. g pe 
s diphenyl ether constitution was arrived at (Chart 5). 
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OCH3 
Carr 6. Synthesis of 2, 4, 6-trimethoxy diphenyl ether. 


Among the products of degradation compound (III) having a diphenyl 


ether structure was important. It could be converted by methylation followed 
by oxidation and decarboxylati 


This could be synthesized by the Ullm 
bromobenzene and potassium phenoxide. 
was proved (Chart 6). 


į 
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Since compound (III) gave positive reaction with Gibbs’ reagent there 
was indication that the diphenyl ether link was ortho to the hydroxyl and the 
alternative possibility was eliminated. We attempted the synthesis of hinoki- 
flavone pentamethyl ether through the Ullmann condensation of 8-iodoapigenin 
trimethyl ether and potassium salt of 4'-hydroxy-5, 7-dimethoxyflavone. 
There was partial demethylation in this reaction and remethylation was done 
to obtain the pentamethyl ether. The resulting product was identical with 
hinokiflavone pentamethyl ether. However, Nakazawa has recently claimed 
that he has synthesized both 4’, 8" and 4’, 6"-isomers by unequivocal methods 
indieated in the formulae given below and only the 4', 6"-compound agrees 
with the natural sample. This is explicable if in the above Ullmann synthesis 
the partial demethylation that takes place leads to isomeric change of the 
biflavonoid system. This is supported by his observation that synthetic 
t, 8"-biflavonyl compound gives natural hinokiflavone when boiled with 
hydroiodie acid. On this basis, the structure of hinokiflavone and the partial 
methyl ethers require revision (Chart 7). 


I 
CH3O o ES CH30 

o n 

QC rer ES 


CH30 0 CH0 O 


Utumann |a Mernvrariou 
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HinoxirLAvoNe PENTAMETHYL ETHER 


Da ENS Rs 
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[xz NO5 3  ArrER CONDESATION —N0;— —NH?2 —H] 


CHART 7. Synthesis of Hinokiflavone methyl ether. 


Hinokiflavone is now found to occur fairly widely in nature; several of its 
mono and dimethyl ethers are also found (Chart 8). 


Cupressuflavone 

While looking for Indian sources of biflavonoids, we isolated from the 
leaves of Cupressus torulosa and Cupressus sempervirens (Figs. 3, 4) a new 
biflavone which belonged to a C-C type and called it cupressuflavone. In 1964, 
it was assigned the constitution of 8, S'-biapigeninyl. Its UV spectrum 
was similar to that of apigenin (peaks at 226, 274 and 330 mu). The intensity 
of the UV maxima suggested that cupressuflavone was a biflavone. IR 
spectrum indicated the presence of chelated hydroxyl (s). Its molecular 
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formula was found to be CgHis0j0. It did not contain any methoxyl or 
O-methyl groups. It gave a hexamethyl ether and hexaacetate. On oxi- 
dation cupressuflavone hexamethyl ether 6 anisic acid as the only 


Crvetomerin A 


LsocrYPTOMERIN 
NeocrYPTOMERIN 


C HAMAECYPARIN 


Caveromerin B 


CHarT 8. Hinokiflavone series. 


Cupressus torulosa, Fro. 4 


—  . 


ANNIVERSARY ADDRESS 605 


isolable entity. The NMR spectrum of the complete methyl ether of cupressu- 
flavone was symmetrical. From the coupling pattern and the chemical 
shifts of the aromatic protons, it was inferred that oxygenation was at 5, 7, 
and 4'-positions. On these considerations the following constitution was 
proposed (Chart 9). 


Cuart 9. Cupressuflavone. 


The decision between these two structures was sought to be made by the 
synthesis starting from 8-iodoapigenin trimethyl ether and subjecting this 
to Ullmann condensation. The synthetic hexamethyl ether was identical 
with the natural sample and from it cupressuflavone itself and its hexaacetate 
were obtained. But the possibility of rearrangement during these reactions 
has not been fully eliminated (Chart 10). 


I 
Gino) O v IN OCH; 


Untmann (Cu) 
ewe © ETHYLATION 
Ac¿0 a 
15 / CgH 
Prmoine s/ Coto 
CupressuFLAvONE 


HEXAACETATE 


CHarT 10. Synthesis of Cupressuflavone. 


Recently cupressuflavone tetramethyl ether (7, 7", 4', 4''-) and dimethyl 
ether (4’, 4'"-) were isolated by Rahman et al. (1968) from the leaves of 
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ii (Fi unds are im- 
Araucaria cunnighamii and A. cookii (Fig. 5). These compo 8 


portant because they exhibit optical activity arising from restricted rotation. 


Fig. 5. Araucaria columnaris (Syn. A. cooki). 


Biogenesis 


The biogenesis of the biflavonoids seems to involve oxidative coupling of 
polyphenols. It has been well known for a very long time that phenols are 
readily oxidized by many different reagents. The products are complicated 
mixtures of dimeric, polymeric and quinonoid types of compounds. This 
method is useful for the synthesis of simple and complex compounds when 
coupled with efficient methods for the separation of the complex mixtures, 


ANNIVERSARY ADDRESS 607 


Phenol coupling has been suggested as an important step in the biogenesis of 
different classes of natural products and has stimulated intensive research. 
In general, the dimeric and polymeric oxidation products can be formulated 
by abstracting one hydrogen atom from the starting phenol and coupling 
together these intermediates by C-C and C-O bonds at active nuclear positions 
o- and p- to the hydroxyl. We owe mainly to Barton (1956) the development 
of these ideas leading to the synthesis of usnie acid and many others by oxi- 
dative phenol coupling. 

The above considerations make it probable that the biflavonyl structure 
could arise by oxidative coupling of the monoflavonoid precursor apigenin. 
The oxygenation pattern of the precursor restricts the interflavonyl linkage to 
6, 8 and 3'-nuclear positions (Chart 11). Jn the formation of ethers, each of 
these positions may be linked with any of the three oxygens of the hydroxyl 
groups. However, among the large number of possibilities only a few are 
so far known and this may be due to some factor of preferential activation. A 
later stage of methylation can yield the known methyl ethers. Support for 
this scheme is provided by the fact that in a number of plants apigenin is 
also present along with the biflavonyl pigments. Some experiments in our 
laboratory using one electron transfer reagents confirm the scheme. 


pe vr OOR 


| J 


[9] HO O HO 0 


(0) 
HO , O HO, O 


HO 0 


CHART 11. Free radical intermediates of oxidative coupling. 


Study of Casuarinaceae 


The study of the chemical components of the gymnosperms, particularly 
the conifers, will be useful for their possible application as drugs. There is 
also interest in regard to the taxonomy of the plants since they contain bifla- 
vonoids as special components. 1 may mention the study of Casuarina as a 
typical example. This genus has about 40 species chiefly found in Australia 
and in Indo-Malayan region, Of these, Casuarina equisetifolia is commonly 
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cultivated in India (Fig. 6). There is some difference of opinion regarding 
the taxonomic position of the genus. Sawada examined the leaves of Casuarina 
stricta and found that hinokiflavone is present in them. Earlier investigations 
had shown that the leaf-wax of this plant gives on hydrolysis juniperic acid 
(16-hydroxy palmitic acid) another typical gymmosperm constituent. The 


Fic. 6. Casuarina equisetifolia. 


occurrence of these two gymnosperm compounds in the leaves of C. stricta 
nt su to the view that Casuarinaceae is close to the gymnosperms. 


atesmi 
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revealed the presence of cupressuflavone. We also examined Casuarina 
suberosa and it contained hinokiflavone. These results may help to settle the 
taxonomic position of the genus (Chart 12). 


Cupressaceae 


C. FUNEBRIS Co [Us Ap H-70CH,. 
C.sEMPERVIRENS ES ¿AS 
C.TORULOSA C, À, H. 


Biota ORIENTALIS 


CASUARINACEAE 


(C. CUNNINGHAMIANA 


C-superosa 


5H 
C—CurressurLAvONE H-70CH3—1socayeromerin 
A— AMENTOFLAVONE H=—HinokifLavone 


CmarT 12. Chemical components of some species of Cupressaceae and Casuarinaceae. 


Study of mass-spectral fragmentation 


During our study of the biflavonyl pigments of Cupressaceae we have 
found the presence of partial methyl ethers of the three known types of bi- 
flavones and also new types of biflavones. They are present in very minor 
amounts and chemical degradation cannot be carried out for arriving at 
their structures. Mass spectra can be utilized in these cases and we have 
therefore made a critical study of the three types of biflavones known. Since 
the hydroxy compounds are very high melting the spectra of their complete 
methyl ethers have been recorded. From such a study, we can distinguish 
diphenyl and diphenyl ether type of biflavones and the position of the me- 
thoxyls can also be fixed. 


Other biflavonyl compounds 


Other modes of biflavonyl linking has been shown to be possible. Re- 
cently Jackson et al. (1967) isolated three compounds of a new dimer that 
could be called flavanonyl flavanone (biflavanone) from the heartwood of 
Garcinia buchanii. This involves the linking of two flavanone units at 3 and 
8 positions. Their constitutions were arrived at by UV, IR, NMR and mass 
spectra. In these compounds both the units consist of flavanone units. 
Morelloflavone, a flavonyl-flavanone, isolated by Venkataraman ef al. (1967) 
from the pericarp of the seeds of Garcinia morella contains a linking between 
a flavanone and flavone. Xanthorrhone and 14-hydroxy xanthorrhone which 


A m 
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from Australian xanthorrhea contain a 


ay 


were isolated by Pelter ef al. (1967) 
flavanone and flavan units (Chart 13). 


DXANTHORRHONE , R=H 


14 — HyoroxyxantHorrHone , R=0H G8—1, ROH; R¿=H 


GB—IA , R¡=R¿=H 


GB-2, R¡ =R2=0H 


MORELLOFLAVONE 
Cuart 13. Recent biflavonyl compounds. 


More widespread and important is another group of this type which has 
also been studied with some clarity during the past ten years. They are the 
proanthocyanidins (Chart 14). Under the name leucoanthocyanidins, they 


TRI-PROFISET'NIDIN 


"TRi- PROCYANIDIN 


CuarT 14. Proanthocyanidins, 


had earlier attracted considerable attention as 
terials and as astringent principles in foods and 
of readily forming polymeric resins and 


components of tanning ma- 
drugs. They are also capable 
à number of them have been isolated 


ANNIVERSARY ADDRESS 611 


in recent years. They are made up of either two leucoanthoanthocyanidin 
units or of leucoanthocyanidin and catechin. But the linking is easily 
broken by mineral acids, yielding eventually the corresponding anthocyanidin 
and catechin. Some typical examples are given in Chart 14. Frequently 
we obtain trimers and higher polymers. 
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NOTE ON RESPONSES IN A PIEZOELECTRIC CRYSTAL 
WITH DIVIDED ELECTRODES 


by P. Roy, Department of Mathematics, Victoria College, Cooch Behar 
(Communicated by B. Sen, F.N.I.) 
(Received. 17 April 1968 ; after revision 8 July 1968) 


The present note is an attempt to obtain the electrical responses in a piezo- 
electric crystal plate with divided electrodes with given electrical input 
employing the method of Laplace transform. 


1. INTRODUCTION 

The piezoelectric problems, particularly those concerning the electro- 
mechanical responses in a piezoelectric transducer, have long been treated from 
the standpoint of equivalent circuit theory, vide Mason (1948, 1950). It is 
only in recent years that such investigations are carried out using the princi- 
ples of mechanics of continuous media and electromagnetism and the pioneer- 
ing efforts in this regard have been made by Redwood (1961a, b). It is perti- 
nent in this connection to refer here to similar attempts made by Sinha (1963, 
1965, 1967a, b), Giri (1966, 1967) and Das (1967, 1969) on the determination 
of responses in piezoelectric crystals. The present note is an analogous attempt 
on a problem as yet unsolved and it seeks to work out the responses in a piezo- 
electric crystal with divided electrodes, some time-dependent inputs being pre- 
scribed. The analysis presented here proceeds on the lines of the mechanism 
of a erystal provided with electrodes as considered by Van der Veen (1956). 
It has been found that the Laplace transform serves as an effective tool for 
the solution of the problem. 

2. PROBLEM 

We take a bar whose centre is chosen as the origin of coordinates (Fig. 1). 
The three axes X, Y, Z are in the direction of thickness (= 2a), the length 
(= 2b) and the width (= 2c). A first set of electrodes x — +4, m <y<0 
and —c < z < c is connected to the terminals J, and J i while a second set of 


electrodes with coordinates x = + a 0 y <b and —c <z < c is connected 
to a second set of terminals Jo and J A 
A longitudinal disturbance is s 


et in piezoelectrica A EE 
dependent voltage VA piezoelectrically by applying a time 


t) across the terminals where 
Vi(t) =0, t<0 


t 
= in 0<t<t (1) 
=Vo, L2 tg 
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and the current induced /, is given by 
I, = H(t) 
A(t) being Heaviside’s steep function defined as 
H(t) =1, t>0) 
=0, t=0) Ga) 


Our object here is to find out the electrical responses represented by the voltage 
Vs (say) and the current J, (say) across the second set of terminals arising out 
of the mechanical coupling between the two halves of the bar. 


Fic. 1. Piezoelectric crystal with divided electrodes. 


3. FUNDAMENTAL Equations AND BOUNDARY CONDITIONS 
The fundamental equations of the problem are evidently those which 
represent the interacting fields. Thus the equations are the constitutive 
relations of the piezoelectric material, equation of motion and the boundary 
conditions. The piezoelectric equations, as in Mason (1950), are given by 


S=sT+dE us io ae .. (2a) 
D=dT+eE an Be as .. (2b) 


where S = 0n/dy, y being the displacement in the direction of length; T, normal 
stress in the length direction; Æ, electric field strength in the thickness direction; 
D, electric displacement in the thickness direction; s, electric compliance; 
d, piezoelectric strain constant; e, permittivity. 
The equation of motion is 
or 02» $ 
By) = Pas EN us ae em 
p being the density of the bar. 
Since the electrodes are equipotential surfaces 


dE 
— we 
ao for y20 be EN ot seo (4) 
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From the continuity of y and 7’ in the plane y = 0, we get 


(1 = (e | 
(7); = (D)o j 


ony=0 (5) 


where the suffixes indicate the respective electrodes. A third set of boundary 
conditions is obtained if, following Van der Veen (1956), we take into account 
the influence of various causes of damping of a vibrating bar by assuming a 
fictitious acoustical medium which exerts on the boundary surfaces (y = + b) 


a stress 


0 è 
(L) x0 = Foa) ds aU DE ane (6) 


+b 
o being a positive constant. 
Again eliminating T from eqns. (24) and (2b), integration over the first set 


- of electrodes J,, J, leads to the total charge Q on the electrodes J,, J, given 


by 
d c 0 d? c 0 
e-i[ a | says (E) | ef E dy 
E -b -c -b 


= {y= -0— =- v} YO 6D ane oo o o (7) 
where $ = 2cd/s. 
The second integral is ¿CpV, where Cp is given by 


Differentiating eqn. (7) with respect to time we get the current through the 
terminals J,, J 4 where 


d d 
I, = $g -o-ta (EnV) oc T so (8) 
The second half of the crystal gives in the same way, 
alld dE vs 
: Is =g g Ub eo) 3; (è Coo) E "> (9) 
Equations (4) to (8) constitute the boundary conditions of the problem. 


4. SOLUTION OF THE PROBLEM 


To solve the problem we take Laplace transform of eqns. (3)-(9) where 


all the quantities with bars indicate the transf i 
ele ranstorm of the corresponding quan- 


de frena [Re (p) > 0] 
0 
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Therefore, all equations mentioned above yield a second set of equations as 
follows: 


oT 2s 
mcum X 6 T T A .. (10) 
Gu] — fo y20  .. "s è (11) 
oy 2 : 50 
mi = (F) fory=0 .. o5 og 2. (12) 
(7); =(D) for y=0 .. a T .. (12a) 
(Duo = Fap() eo (13) 
or 
9n e d = ; 
(55) eme zT. (13a) 
07 $ d= 
(52) tev» = 3, Va (135) 
I, = $p(5-o—5-0)H-32C5V1 be DD - (04) 
and 
I» = $p(55—54o)H-3pC5Vs e ao - (15) 
Solving eqn. (10) we get 
y -8 
5-—4e"-FBe " fory<0 
(16) 
by 2 
=Cev+De " for y>0! 
where p = 1/2, 
The boundary conditions (12) to (13b) yield respectively 
AFB=C+D REL M uu (m) 
B dv > p 
A— —C+D Soma) .. .. se (170) 
AT VA 
Ae A on a ae (17c) 
pb p 7 
Tr ar ÍA 
Ce K — De Kı— 2ap =0 .. n .. (17d) 


where 


` 


Kı =1—avs, Ke =1+avs 
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From eqns. (17a)-(17d) using eqn. (14) A, B, C and D can be obtained as 


follows: 


N 


#7 ph hy pho 
o= ixfar (1—Kae P) (Kee? + Kae *)- (ie "—Kze")} 
p 


SUAM Se n = ee 
ee Y, (1=K0 >) le 0) (o +es—2))] . (8) 
D y 


dv = 
A 
_ bb pb bb 2 -2 
D= malar] (Ke "- et )- (Kye? —1) (iot +0 Jl 
do = HA ph bb a 2 | dv 
D HA ?-pe" —2— | Koer —1 6 v—e f o 
where 
E pb 
L-1 =2—K;e ¥—K,e" 
=P) pb QU _ pb 
Mx =e "i (1e | + Kae (1=e ai Tc (1:9) 
= 2 pb 
4g = ke 1) ae) 
and 
I Y 
M —-——105—. 
pp “74 
and from (1) and (1a) 
= V Z 
ar e o (200) 
and 
se 20b 
=> . (205) 


knowing the values of O and D, Vs and I, and hence V, and J, can be com- 
puted. Here we only compute Vs. 
From eqn. (17d) 


509, po =D 
LES [cerne "x; 


Substituting the relevant values of C and D we get V.(p). This inversion is 
very complicated. "Therefore, proceeding as in Redwood (1961a) we suppose 
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t small, i.e. p large, such that higher powers of 1/p can be neglected and we get 
on simplification 


y, 2 Ll Ko =| 

2 dw $|220—K$,) Kı 
2ap | Ks T E 1) E 
dv y p — K, (2ap K m 


e elo. as 3] 
7 A (eese 


vd$ép|2(2—Ks.) Ko 
e Pace (1) ele 
dv plo 2 2—Ko Kəf 2ap 
Ks Kil >| 
tax (0 
— Pu i Koa S ol 
do óp|9(2—Ks;) Ko 
2a Vo al) Ke Se ee 
^ dv ty |[2al VK IK 0 Ke p2 


p 
Inverting this we have, vide Churchill (1958), 


2a Kı Kı 
E TE r 
Vat) = dod E 50 
2a Vo 2 2 ) Ks _& uf 
-5 ^ E (En K. K tot 
( Ko Ks P 
+01 Kop 21 A 


Z3 sl Ky El () 
— dv) |9(9— Ke) 


A(t 
2a Vo dv (e+ ae Co F 
-pe Po] ge ial ISI ot Ve oM 
Thus we find that the electrical voltage response has different characteristics 
for different ranges of time, namely part of it is step in character throughout 
while another part of it is linear in 0 < t < to and this latter part continues 


to be constant beyond time to. 
This expression V(t) enables us to calculate Za(t) from eqn. (15). 
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SOME SPECIAL CASES OF SPHERICALLY SYMMETRIC 
GRAVITATIONAL COLLAPSE WITH REFERENCE TO 
NEUTRINO EMISSION 


by K. P. Sinan and M. C. Srivastava, Department of Mathematics, 
Banaras Hindu University, Varanasi 


(Communicated by B. Dayal, F.N.I.) 
(Received 17 April 1968; after revision 1 July 1968) 


Misner's field equations (1965) for spherically symmetric gravitational collapse 
with escaping neutrinos have been considered in some particular cases. It is 
shown that spherically symmetric gravitational collapse with neutrino 
emission is not possible if the congruence associated with the matter field 
satisfies any two of the following properties: (i) it is geodetic, (ii) it is 
expansion-free and (iii) itisshear-free. Also the field does not permit neutrino 
emission if the null congruence associated with neutrino emission is expansion- 
free. A case has been worked out with geodetic matter congruence where the 
possibility of neutrino emission is demonstrated. Similarly it is shown that 
neutrino emission is possible in the case of shear-free matter congruence. 


1. INTRODUCTION 


It has been suggested by Colgate and White (Misner 1965) and also by 
Chiu (1964) that in case of gravitational collapse of supernovae and quasi- 
stellar radio sources the emission of neutrinos is possible. Misner (1965) for- 
mulated the field equations for a spherically symmetric case of gravitational 
collapse with escaping neutrinos. It has been assumed that the neutrinos 
after emission travel in the radial direction. In the absence of neutrino emis- 
sion the field reduces to that obtained by Misner and Sharp (Misner 1965). 
The internal field of collapsing matter is continued with the field of pure radia- 
tion such as the one obtained by Vaidya (1953). In this paper we discuss 
some of the special features of the field equations of spherical gravitational 
collapse with reference to the possibility of neutrino emission. It is shown 
that spherically symmetric gravitational collapse with neutrino emission is not 
possible if the congruence associated with the matter field satisfies any two of 
the following properties: (1) it is geodetic, (ii) it is expansion-free and (iii) it is 
shear-free. Also the field does not permit neutrino emission if the null con- 
gruence associated with neutrino is expansion-free. A case has been worked 
out with geodetic matter congruence where the possibility of neutrino emission 
is demonstrated. Similarly it is shown that neutrino emission is possible in 
the case of shear-free matter congruence. 
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2. "THE FIELD EQUATIONS, KINEMATICAL QUANTITIES AND 
OPTICAL PARAMETERS 


We give below the system of field equations and define the various sym- 
bols involved. The meaning of the symbols is explained subsequently. 
The field equations are as follows (Misner 1965): 


A = (nu), 4 —0 uet e. an at ao — (XB) 
B= —no-- ue Ts. ,=0 N, 2 jo .. (2,2) 
Di = R,—485;R-4-82(T,--NP) 20 .. fe .. (2.3) 
where 
E =uwua+l=0 m de 25 59 .. (2.4) 
Ti = (e--p)utup--97p .. E Sa o so (AG) 
N; = gk@ky 00 og 0-0 om o dre (2.6) 
and 
F = kika = 0. x a X tu me (27) 
Next we consider the following kinematical quantities associated with «^: 
VS Ua +. n eS t ae F y T (2.8) 
wt = Anrep, d 00 Oi dio v5 0m .. (2.9) 
— Way = F(a; b— Uv; a) HA (UUp—ÚDU A) .. T o .. (2.10) 
Gap = F(a; Up; a) +A (aus gus) — H (Gap+ tatty) .. (2.11) 
where 
tha = Ua; pu? 09 9.6 2 "S as 66 +. (2:12) 
and 74'c? is the pseudo tensor density with $1284 — sedet ’ 
=g. 
Also the ‘optical parameters? are 
v= Mia. (2.13) 
DA ooa (24) 
Z= Ha; lei... T S Si: 2 .. (2.15) 
Q = Zka; plime .. m v ats M. .. (2.16) 
where 
Ja = fe” kb = 
(o c IE (P= oa a .. (2:17) 
t= 
Po (2.18) 
ala = 
PO san e er .. (2.19) 
E: L0 on $3 do T 2: B ». (2.20) 
af, = 
; a=0 m o0 se 55 <. (2.21) 
ap, = 
em oo (222) 


R 
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loma — 0 eo o£ GHI MS 
lama = 0 X E: e E re a .. (2.24) 
lla—l =0 WOT, —— 3 
k*mg—1 = 0 T ox AC 3: T. T .. (2.26) 
and 
Hay = Jar—(laly+lola)—(kamotkyma) =0. .. .. (2.27) 


Equation (2.1) is the ‘equation of continuity’ for matter, eqn. (2.2) 
correlates the cooling rate of unit volume of matter with the rate of decrease of 
internal energy due to neutrino emission and eqn. (2.3) represents the 
space-time structure in terms of total stress-energy. Here n is the baryon 
number density and c the cooling rate of unit amount of matter. The assump- 
tion of perfect fluid distribution form of the stress-energy tensor Tay for 
matter is inherent in (2.5) while (2.6) gives the ‘geometrical optics’ form of 
the stress-energy tensor Nay for neutrinos. Here e, p and q are respectively 
the matter density, pressure and neutrino energy density. 

The kinematical quantities (Witten 1962) for an observer following one of 
the curves of the time-like congruence and using Fermi-propagated axes are 
described as follows: vr, way, cay are respectively the velocity of expansion, 
rotation and shear of the neighbouring cloud of particles. The vector w is 
the angular velocity in the infinitesimal rest-space of the observer. The 
“optical parameters’ (Witten 1962, p. 58) vrn w, 2 and Q are called the ex- 
pansion, twist, shear and rotation of the null congruence with respect to the 
Observer. A semi-colon (;) preceding a suffix indicates covariant derivative. 
Here we have 37 equations: 


1(A)+1(B)+10(D$) +1(E)-+ 1(F)-+4(G2) +9fegns. (2.18)-(2.26))+10(Ha») 
in 35 unknowns: 


1(n) +4(u2) 4- 1(c) J-3(e; p, q) + 16(10, 12, ka, ma) + 10(gay). 


3. THE CASE or SPHERICAL SYMMETRY 
We consider the spherically symmetrio metric 
ds? = —e dt? + ef dr2+-y2(d02-+sin? Odee) — .. — .. (31) 


where «, B, y are functions of r and t. The coordinates r, 0, ¢ and t will be 
denoted by 2}, 22, x3 and x4 respectively. 


u=0 fora=1,2,3; ut =e, T. «+ (3.2) 


Next we choose the orthonormal tetrad given by 


(6) = diag (V/gu, Vg, /g8, y —git) 
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and pseudo-orthonormal tetrad (1%, la, k^, d is defined as 
e= (CS 1» w=, z 2:9) 
(3.3) 
po ger m= pe 
where a suffix followed by a vertical stroke indicates the particular vector. 


Equation (2.27) can equivalently be written as 
Hay = Jab— E riii Ay adso =0. 
For the metric (3.1) the non-vanishing components of (12, ja, ka, ma) are 


vi 


I; = lg]i sin 8 = Ty = —Ig/i sin 0 = 


| ebl2 2 ex re) 
bi =m =—=, k Ms = — —= 
1 1 Y? 4 a4 NE: 
and due to (3.2) and (3.4) eqns. (2.1) and (2.2) reduce to 
e-Al2 
n = P(r) yi (3.5) 
= e-% Ma -« $4 
c = e-*(e--p) up s dt T .. (3.6) 
where F(r) is an arbitrary function and n, , = cem 
The non-vanishing components of Day in (2.3) are 
rn AA y 
Dig = —8m ¿+2e7*-PÍ2 (#2 y ee > s) =0 (3.7) 
8 q “oa | 2Y44 274 75 hs Y 2y 1 
DiS (s) i era. y utet esty — 7 =0 
(3.8 
= Dis emp ect y aa +% Bs, Bas Bi 2 
s Y: Yı B De 
—e (- ELE "i S E o, o Bt -d-«) =0 .. (8.9 
2 
Dj = &r(«+-4) — e Ya ]- (- £u * Y1 ) ges 
: pty Pale y yty 73 =0 
UR .. (8.10) 
oy 
ME Or yam 2" ote. 
Equation (2.17) reduces to 
gn oy, = = 0, "E & s .. (3.11) 


21 
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Here we observe that for spherical symmetry «^, way and 2 vanish 
identically, whereas w and £2 vanish due to (3.11). Consequently, the con- 
gruences associated with u and k@ are necessarily normal. The non-vanishing 
kinematical quantities and parameters are given by 


: -al Ba, 24 j PED 

vp =e (E? ae Js x: 5% Sí .. (3.12) 
9eB MEA CAS 2y, 

oy = ye G92 = — y2 sin? 0 = 3 (8. y Qo .. (3.13) 

Uy =o, lg=Ug¿=U¿=0 .. m T me .. (3.14) 
1 LA Ya * 

ee lo HB tease), o» Ws a6 A 

VII Và ( jr z) (3.15) 


Thus we have seven equations, viz. (3.5) to (3.11) in eight unknowns, viz. 
n, C, €, P, q, % B, Y, along with vr, vrm 911, 1. In the next section we study 
the consequences when some of the quantities vr, vrr, 911 and ù, vanish. 


4. SPECIAL CASES 
The eqns. (3.7), (3.8), (3.9) imply 


pa ra Ta EE Yaa Ys Ys Ba Bag Bi B4 
a (s 70 Fs Pas pone 4 y ty tty oT Hs pty "ug 


1 
hen - Yang a Bg % ata T ER A=0 (4D 


Now we consider the following cases: 


Case (2) —v, = 011 = 0. 
Equations (3.12) and (3.13) give 


Ba 27. 
gp co 
2Y 

B oo 
which imply 

B4 —Y4-—0 
and hence from (3.7) we get 

q — 0. 


Since e-*Y, is the velocity of the fluid surface, Y, — 0 implies no collapse 
while q = 0 implies absence of neutrinos. Hence expansion-free and shear- 
free (or rigid) congruences associated with matter field imply no spherical 
gravitational collapse with escaping neutrinos. 


Case (13) —ú, = 011, = 0. 
Equations (3.14), (3.13), (3.11) and (3.7) imply 
=f, =Y4=g=0. 
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Hence geodetic shear-free congruences associated with matter field do not 
permit spherical gravitational collapse with escaping neutrinos. 

Combining cases (i) and (ii) with ù% = 0 we have: 

Spherical gravitational collapse with escaping neutrinos does not admit a 
space-time which is the direct product of a time-like line and a 3-space (Witten 
1962, p. 58). 

Case (44) —ú41 = v; =0. 
Equations (3.14), (3.12), (3.11) and (3.7) imply 
4 =f=Y=9=0. 
Hence geodetic expansion-free congruences associated with matter field imply 
no spherical gravitational collapse with escaping neutrinos. 
Case (10) —U;1 = Me 
Equation (3.15) gives 


O wo ss (&2) 
From (4.1), (4.2) and (3.9) we get 
1 
cof -u eere t -&u) E. ACE 6 1 (43) 
and 
ox te Y, 
eo Ha) +e-o(— mun S. .. (4,4) 
Equations (4.2), (3.7) and (3.10) give 
5 n Y 
8re =e «(ra "a teef- mn A .. (4.5) 


which along with (4.4) implies 
€ — —P. 
Since e and p both should be positive we have 
e= p = 0; 
then field equations (3.7) and (3.8) give 


-o A Y 
—9eg-« (rus — fl Ba = Yam) = eta Pas Pin) 3002 1 


which by virtue of (4.2) reduces to 
1 
7 = 0: 

This is not compatible with the metric. Hence geodetic expansion-free null 


congruences associated with neutrino emission im i itati 
e i; ply no spherical 8 
collapse with escaping neutrinos, 7 1 DAE 


OASES OF SPHERICALLY SYMMETRIC GRAVITATIONAL COLLAPSE 


Case (v)—34, = 0. 
Equation (3.14) gives 


and (3.11) gives 


Then (3.1) can be written as 
ds? = —f(t) di2+ g(r) dr2--Y?(r, t) (d02-- sin? 0 dd?) 
which can be transformed to 
ds? = —di?--dr?--Y2(d62--sin? 0 d$). 


For this metric we have 
wm=w=w=0, ut=1l 


lg 3 ls = Y A " =) = L 
Cno  isnó 9/5” By m 24-5 
F(r l 
n = ZQ > c= nF (Y14(Y1 +Y4) -Y(Y1147-Y144)) 
Y14 
2 = — 
ag y 
q Y y? y? 1 
AV 44 4 1 
~8x(+4) => typ piz 
Yas Yu 
—8 = ET 
TD 3 y 


vr = 2Y4/Y 

011 — —2oss[Y2 = —2¢33/Y2 sin? 0 = — 274/37 
_ ll ftir, % 

on = VAT 


and (4.1) reduces to As 
Y(Yii-F2Y14-4-Y44) +y Y) +1 = 0. 


One solution of (4.7) for which neutrino emission is possible is 
y =f(lr+mt) 
where l, m are constants and function f satisfies the equation 


(L4- m)2ff" + (m8—12)f'2--1 — 0 
where 


jy e et 
f =a u =lr4mt. 
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(4.6) 


. (4.6a) 


. (4.6b) 


. (4.6c) 


. (4.63) 


(4.7) 


(4.8) 


ll TUNE 
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For this value of y, q, p and e are given by 
27q = lmf" if 
Sap = (2—m?)f"[f 


Sre = — 2)? 7 —2lm E+ mt) Ets f 
If p = 0, which implies a distribution of discrete particles, we have from 
| (4.60) 
| Ya —Y11 = 9 


for which the most general solution is 


y =f(r+t)+9(r—2). 
Equation (4.7) reduces to | 
Af (f--g) 4f'g--1 — 0 UNS L5) | 

where 


Differentiating (4.9) with respect to (r—t), we get 


p 
AO 

or 

e = y = k (constant). 

yg 
Then 

ern) +41} efk(r—t)-+k3} 
=< +h, gh. 


Equation (4.9) becomes 
Lehr +N H1) -- A (y+ kaert — 0 | 


which implies exp (k(r--£)--À;) = constant. 

The above equation on further examination leads to the conclusion that 

discrete particle distribution with geodetic congruence associated with matter 

is not consistent with spherical gravitational collapse with escaping neutrinos. 
In particular, when g=0, e and q can be calculated directly from 

(4.62, d). In this case it is found that q turns out to be negative. 


Case (vi)—o,, = 0. 
Equation (3.13) gives 


which on integration gives 
ofl? = g(r)y(r, t). 
Now (3.1) becomes 


ds? = — e% di2-4 y2 2(92 dr2 + d02+sin2 6 dd?) 
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which can be written as 
ds? = — e« di24+-Y2 (dr?-]- d02-- sin? 0 d$?). s .. (4.10) 


For this metric we have 
wm=w=w=0, ut=e 


led = Ll =l = 2,3 =m, =M3 = 0 
ls - la Y 
ag eS i sin 0 imu cUm 
ex 
kg = —Ma = — — 
v3 
n= 2 


CxS 9, 9 = oes 2 

c= SaP [vo —14 tety Geared de *y14—26-9*yy. 
Dana 3Y11 dv; 9, 

PA c y tyk — 20 %Y y ,--2e YY 14] +2Y 114 


2yiY 
— PA aue dta 


meu m y yi y 
2 ? 9 Z 
—8e(+4) = cera tr 2 Exec L (% = 2% ge) a 
2 MIS Y Uo MIT Y Y" 
-ox| 2 Y 2y 1 CAM ÍA a 
7 8mp =e «( ALL Se Al ESAME 
E y Uy c y Pasce alu tue 
2 2 
va, ul 2 Y 1 
Sm[«+2) CIEN | Yn 5 
a EAU > PY Ya 
Vr = 3e-*y,[v 
VO PREDA 
vae y) 
and (3.11) and (4.1) become 
Yat ex = (0) ou Is cO .. (4.11) 
2 
a A a : 
€-*y14--2e Ud =0. bu .. (4.12) 


One solution of (4.11) is 
yz—Ff e«=Pf 
where 


SEX 
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Equation (4.12) reduces to 
go go Ta Ae 


m ON 
and 
1 
anq = mp 
ME Eu Up ge 
SIS ptem Fo Ta) tapre 
Sze 


1 1 9p" p 1 
= pt paa (- p = E 


which demonstrates the possibility of neutrino emission in this case. 
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COUPLED VIBRATIONS OF A SLENDER BEAM IN A 
CENTRIFUGAL FORCE FIELD 


by J. S. Tomar, Department of Mathematics, University of Roorkee, 
Roorkee 


(Communicated by Jai Krishna, F.N.L.) 
(Received 4 May 1968) 


The differential equations for the coupled bending and torsional vibrations of 
a slender beam in a centrifugal force field are obtained. A method based on 
the Rayleigh's quotient is used to obtain an upper bound for the fundamental 
frequency, and the accuracy of this upper bound is investigated by considering 
tho application of the method to uncoupled vibration. 


NOTATION 
x = radial co-ordinate 
Y = tangential co-ordinate 
z = axial co-ordinate 
v = deflection of elastic axis in Y direction 
e = axis of minor moment of inertia J min 
v' = deflection perpendicular to axis € 
l = length of the blade in radial direction 
m = mass per unit length of the beam 
E = modulus of elasticity 
I = moment of inertia of the cross-section about the axis perpendicular to 
the plane of bending 
G = shear modulus of rigidity 
J =a constant depending upon the cross-section of the beam such that GJ 
is the torsional rigidity 
3 = distance between the elastic axis and the centroidal axis 
0 = rotation of the beam cross-section about the elastic axis 
S = area of cross-section of the beam 
Ip = mass moment of inertia about the elastic axis per unit length 
t = time 
M = moment of centrifugal force at the point x 
£ = dimensionless variable £ = xjl 
w = frequency of vibration 
INTRODUCTION 


The analysis presented in this paper considers vibration of a slender beam 
that could represent a turbine blade of simple geometry. The beam is attached 
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to a disc of radius rg as indicated in Fig. 1 and the disc rotates with angular 
velocity 2. The shear centre of each cross-section does not coincide with the 
centre of gravity, consequently the torsional and bending oscillations are 
‘coupled’. 


Fic. 1. Cantilover beam attached to a circular disc. 


For the analysis it is assumed that S, J and J are constants. This, how- 
ever, is not an essential assumption. If S, J and J are functions of a, the 
amount of computing work is greater than for S, J and J independent of x. 


THE DIFFERENTIAL EQUATIONS 
The differential equation for the deflected form of the neutral axis of a 
bar according to the elementary theory of bending is 


Q4v 
EL" Wp ya a os (Qi) 


where w is the intensity of the distributed load. 
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If the load is distributed along the centroidal axis, the given load can be 
replaced by the same load distributed along the shear-centre axis, and a torque 
of intensity w . ò distributed along the same axis. 

Let the z-axis coincide with the shear-centre axis. Since the torsion is 
not uniform, the relation between the variable torque T' and the angle of twist 
0 is given by (Timoshenko 1955) 


90 030 


where GJ is the torsional rigidity for uniform torsion and C; is the warping 
rigidity. Differentiation of this equation with respect to x gives 


020 ~ 030 ; 
GI ga gri UO Sus Us go we (2) 
with the positive torque taken as shown in Fig. 2. 


Z 


Elastic Axis 


(b) 


le) 


Fra. 2. Deflection and rotation of a cantilever beam from and about elastic axis. 
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For a vibrating bar the intensity of the inertia force is 
: as 
-mp (v--8 . 0) 


and the intensity of the inertia moment about the v-axis is 
020 . 
=Ig E 
The following differential equations for the coupled bending and torsional 
vibrations are obtained by replacing the statical load in eqns. (1) and (2) by 
the inertia forces, 


0*v 02 ; 
A rem 5. zn a (8a, 
EI opi my (w+. 0) (3a) 

020 040 02 029 ; 
AN ONE m AE : == 3b 
da 01 P m.835 (v4-8 9) 3s (35) 


Referring to Fig. 1 the centrifugal force dc exerted by a mass element of 
the radial length dy at x = y is 


dc = mR? (1 +17) de 


If the deflection of the blade at time £ is v'(v, t), the tangential and radial 
components of dc are 


v' sin x 
decy = d 
Cy 7 Em 
and 
D: r- 2714 
v' sin x 
d= ao j- (03) 
fz a| e 


respectively, where X is the angle between Y and e axes. 
If it is assumed that sin X < 1, dez = de. The components cy and c; 
applied at x follow by integrating from y = x to y = I and are 


Cz = mQ?[ro(L—2) 4-3(02 —32)] * d .. (4) 
i 
cy = MR? sin xf v'(x, t) dn ae zs 25 (6) 
Referring to Fig. 3 the moment dM about the minor principal axis of 
inertia is 


dM = cxdv'(x, t) —cy sin x dz 5s zs .. (6) 
Eqns. (4)-(6) give 


02M Q2y' f 
Ou m? Mrt- A —(ro+2) oe sin? xv' i 
If x is a small angle and v’ = v this may be written as 


02M 32 
E ms? [rot—2) +4029 a — (rotz) e| of Sm) 
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The second derivative of the moment M with respect to x can be considered 
as a lateral load depending upon the centrifugal effect. 


M *dM Cc 


dv 
dx 


Fic. 3. Showing moment about the minor principal axis of inertia. 


The R.H.S. of eqn. (3a) represents an T loading, and eqn. (3a) 


must be modified by adding the lateral load, if the centrifugal force 


x dU , 
effect is to be considered and the governing equations then become 
Di 02M 
BI = -m zz > ta. 0+=> ETE 
020 010 0*0 Q 
a ici T Ee S 
Sl z =C1 JA im. Z(t. 9) +19 s (8) 


OM. . 
where E given by eqn. (7). 


DETERMINATION OF NATURAL FREQUENCIES 


The equations are now put in terms of a dimensionless variable = x/l 
and the substitutions, 


are used. 
Eqns. (8) become 


0?v 0?0 — 3 9 
ag tan —° > 
„20 g 28 0% (59, go 20 
ie p aa O +10) 73 =9 
The solutions of eqns. (9) are of the form: 
e] = Aree | A E 
6(£, t) = Belge 


where A and B are constants which are not independent; and f(£) and 4(&) 
are functions of £ only. 
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The functions f(£) and 4(£) satisfy all the boundary conditions of the 
beam wbich are as follows: 


«(0) =3 0 = 6(0) = 53 (0) =0 | 
a2y 83v 00 930 i e (LL) 


20 = aga (D = ag (0 = 3 (1) = o] 
For an approximate determination of the fundamental frequency F(€) is chosen 
as the shape function for the fundamental mode of uncoupled bending vibration 
and ¢(é) as the shape function for the fundamental mode of uncoupled tor- 
sional vibration of a uniform cantilever beam. These shape functions satisfy 
the boundary conditions (11) and are 
f(E) = cosh A£— cos A£—0-73410(sinh A£—sin A£) 
pnd B (2) 
plé) = sinz é 
where À = 1:87510. 
Substitution of eqns. (10) in eqns. (9) gives: 
ER cally Dg l JE 2 (70, 4f 
[eee —Q (2; 7 é 56 qp? mant dé —w?f | A—89?4 B = 0 


2 


(13) 
Sw2f, A+ eas — 0,544 regen B=0 

Equations (13) can be solved for w? but the result is a function of £, since f 
and ¢ are not the exact shape function. This difficulty can be overcome (Fung 
1955) by multiplying the first and the second of eqns. (13) by f and ¢ respec- 
tively, and integrating with respect to £ from 0 to 1. The method results in 
the familiar Rayleigh quotient (Collatz 1960) when applied to uncoupled prob- 
lems, and is an extension of Rayleigh’s method to the coupled problem. The 
following equations are obtained: 
(ay — ao +43 — 49?) A —ago?B = 0 
—(4W024 + (47, +ag—agw?)B = 0 j a) 


1 


ro, l 
sce tote iet oa 


where 
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1 
"=a | $.f.d£ 
0 


ER ,| d$ e dd 
m=) Gad a J. a es 
i 
_o,| 29$ 
cm 02] Tta 


1 
(lg = (425 $ p2dé 
0 


For a nontrivial solution A and B must not both vanish, consequently the 
determinant of the coefficients of eqns. (14) must be zero: 


(1— (5-3 — (gw? — (gw? 
— gw? (iz +dg—Agw? 0 
i.e. Pot — Qo?-- R =0 .. = o6 (15) 
where 
P = aqüg—agte, Q = (a1—42+43)a9+ (a7 +0g)a, 
and 


R = (01— 02-48) (07 +45) 
The solution of eqn. (15) is 


QV (Q3 — 4PR) 


2P 


m 


(16) 


The right-hand side of eqn. (16) is positive since it may be shown that 
Q?—4PR > 0. 

The smaller of the two values of œw? given by eqn. (16) is an upper bound 
for the frequency of the fundamental mode of vibration. The larger of the 
two o? values is an upper bound for the next higher mode of vibration. 

Also, for the larger value of w?, from eqns. (14), A and B have opposite 
signs and for the smaller value of œ? they are of the same sign. The cor- 
responding two configurations are shown in Fig. 4. 


NUMERICAL EXAMPLE 
A numerical example for the coupled torsional vibrations of a slender 
rotating beam is now presented. The frequencies are computed from eqn. 
(16), and the cross-section of the blade is taken as a semi-circle of radius do 
and thickness ¢,. The physical constants of the blade except the additional 
assumed value of 3 are taken from Grammel (Biezeno and Grammel 1954). 
= 4-41 in E = 29-20 x 108 Ib/in? 
S = 0-1488in? I — 0:001845 in* 
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dy = 0:23885 in 


tı = 0:19706 in 


m = 0:00011 Ib 


ro = 29:10 in 


G = 11:53 x 109 1b/in? 


2 
J =0-00193 in = rs 
ô= 0-15in 
—314sec^i 


9 


Ig = 0:14834 m lb[in? = m ie (204015) | 


T7 


3 8 
C, = 0:00014 x 109 Ib/in* = Zag 5t, (5 - B 


(a) 


(b) 


Fie. 4. Configuration corresponding to two values of w2. 


With these values, the various parameters used in the frequency equation are 


a8 follows: 
44 = 16:20711 x 109 
45 = 0:30724 x 109 
as = 1-44909 x 108 
a=1 


and from eqn. (16) 


as = dg = 0-10151 
à; = 12-75067 x 109 
Gg = 0:01021 x 106 
ay = 0:08542 


t; = 170830 x 107 soo”? 
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and 
w = 172511 x 108 sec"? 


which are the upper bounds of the frequencies corresponding to the first two 
modes of the torsional coupled vibrations of the rotating blade. 


Discussion AND CONCLUDING REMARKS 
For the uncoupled case (i.e. 3 = 0), eqns. (13) are two separate eigenvalue 
problems that satisfy the conditions of self-adjointness and full definiteness 
(Collatz 1960). The first of eqns. (13) when 8 = 0 is 
2 Uf 2 Das UE n 2) - To Jsl- 2f 
PASE (; pi E NS ce 
with boundary conditions f(0) = f'(0) f'(u)2f"0)29 
This can be written as 


Ulf] = oV [f] 

where U and V are differential operators and the boundary conditions as 
B[/]=0. 

A sequence of functions fi, f2 fa... ean be obtained from an arbitrary fo and 
application of the boundary value problems 

Ol fel = VUr-11 

BL fe] = 0. 
If the ratio fx_a/fx tends to a constant as k increases, the sequence converges, 
and the function fy converges to the form of the first eigen-function f. The 
quotient U[fx]/VL fe] should provide an approximation to w,, Where ws is 
first eigenvalue, but as it is still a function of € the numerator and denomi- 
nator are each multiplied by fr and integrated with respect to $ from 0 to 1 to 
give Rayleigh’s quotient which is an upper bound for «,, that is 


1 1 
A Uscuddé 
_ 0 — 0 
a 1 
| nvisad | AVLE 
0 0 


inem tls 25 asc 


A 


The form of the integrals appearing in the above expressions for 4| may be 
denoted by 


1 E 
basa = f Se V feild 
box = | AVLE (Em 1,200 


asl F 
oki = NISUS 
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and are called Schwarz constants (Collatz 1960) and 4 is the quotient of two 
successive Schwarz constants 


Dorx-1 
A = Hop = — 
Por n 
The quotients por) = Vor_o/Dor_1, per = Dor-1/Dox, pony = Dor/Dor yr ... ave 


known as Schwarz's quotients (Collatz 1960) and the even numbered quo- 
tients are identical with Rayleigh's quotients. 

Bounds for the first eigenvalue o. can be obtained from the quotients 
Hx, k+, and la where lo is a lower bound for the second eigenvalue and such 
that la > rg. The bounds are given by (Collatz 1960) 


Py — Ek 2 9 « 

Brg =I Su < Mey (Eo 1,2, 3...) 
m | 
Prt 


For numerical computation of the quotients fı is taken as the shape function 
for free vibration of a uniform cantilever, that is 
fı = cosh 1:8751€—cos 1-8751£—0-7341 (sinh 1-8751£—sin 1:8751€) 
J and is a non-zero function and satisfies all boundary conditions and possesses 
continuous derivatives. 
Also, a function f; such that 


Ui fi] = VLfo] 
or 
Uf aonn A a) def, i; a 
poo ei Seca) Sh (22) a} — 4, 


can be readily obtained. 
In this special eigenvalue problem with V[f] the condition 


[i foU) Ifa] ae = o 


is satisfied and according to Collatz (1960) fọ need not satisfy any boundary 
conditions. 


Using the same numerical values for the physical constants of the beam 
already considered but with 8 — 0 Schwarz constants and. quotients are 


1 

bg = | fo d£ = 3:02437 x 104 fy = Eo = 1-74325 x 107 
0 1 
; b 

b= | fofi dé = 1-73490 x 107 po = 7" = 1-73490 x 107 
; 


1 
by = i fidé =1 
0 
and an upper bound is 


Hy 2 Fz 2 w? ie, 1774325 x 107 > 1-73490x 107 zw? 
3B 
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For calculating lower bound from the expression 
[yee 
~  Is/Ho—1 
second eigenvalue and even a rough value for la is justified since changes in la 
have little effect on the lower bound calculated from the above expression 
when l is appreciably greater than py. 

Further, it is quite reasonable to assume that for flexural vibrations, the 
second eigenvalue for free vibration of the uniform beam is a good enough 
approximation for the lower bound of the second eigenvalue of the rotating 
cantilever beam. Hence, the calculated value for la is 6:36515 x 108 sec" ?. 

Thus, the frequency for fundamental mode of uncoupled vibration is 


bounded as 


< o, it is not essential to have a close lower bound /5 for the 


1-734067 x 107 sec”? < w? < 173490 x 107 sec? 
The above result shows that the method provided very narrow limits for funda- 
mental frequency of uncoupled vibrations in the presence of a centrifugal 
force field. 

The fundamental frequency for coupled vibrations could not be bounded 
at present, but the above computation indicates that Rayleigh's quotient 
obtained from the shape function f; is a very close upper bound for the first 
eigenvalue w for the uncoupled problem. Consequently, it seems reasonable 
to believe that the method results in a close upper bound for the fundamental 
frequency for the coupled problem, especially if è is small. 
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PETROGRAPHY OF THE TAL SERIES, SOUTH OF THE VILLAGE 
KOTIDHAMAN, SIRMUR DISTRICT, HIMACHAL PRADESH, 
WITH A NOTE ON ITS PROBABLE ENVIRONMENTAL 
CONDITION OF DEPOSITION 


by S. KansiLaL, Department of Geology, Banaras Hindu University, 
Varanasi 5 


(Communicated by Dr. Rajnath, F.N.I.) 
(Received 22 February 1968; revised 9 September 1968) 


Petrography of the Tal Series south of the village Kotidhaman (77° 33': y 
30? 38:5") has been studied here in detail. On the basis of this study they A 
have been subdivided into three units as against two by Audon (1934). The 
study also revealed that the rocks originally described as “graywacke” and 


‘quartzite’ are actually siltstone and a variety of sandstone with a very high ii 
detrital quartz percentage. A probable environmental condition of deposi- | 
tion of these rocks has been suggested on the basis of characteristic lithologic j 
association and the various sedimentary structures present in the area. 


INTRODUCTION 


The area under study forms the southernmost part of the “Western Basin’ | 
of the “Krol Belt’ described by Auden (1934) in his classic memoir, “Geology of j 
the Krol Belt’. There is a beautiful development of the Tal Series in this area. 
In the present paper a petrographic study of the Tal Series, as occurring to the 
south of the village Kotidhaman, Sirmur District, Himachal Pradesh, and lying 
between longitudes 77° 32' to 77° 35’ and latitudes 30° 36’ to 30° 39’, has been 
made. Besides the Tal Series some formations of the Krol Series also occur in 
the area, but they have not been included in the discussion here. 


Ra mM o 


GEOLOGY OF THE AREA 


The topography of the area is mountainous and highly rugged with very 
prominent relief, characterized by escarpments, high hills dissected by deep 
valleys and spurs (Fig. 1). Except a few, most of the streams are dry during 
the summer season. The height of the hills gradually rises toward north. 
The Tal Series, which is a very thick pile of partly argillaceous and partly 
arenaceous sediments, occurs throughout the area except in the south-w 

. portion (Fig. 2) where the Krol Series is found outeropping. The rocks 
comprising the Tal Series are folded into a syncline. The beds are found 
dipping at 25°-35° due south-west, in the northern side of the area, while on 


the western to southern sides the dip direction changes from N.N.E. to almost | 


estern 


e how” RE ee n 
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north and the amount of dip varies from 40? to 50%. The Krol Series, too, 
correspondingly changes its dip direction from N.E. to north. 


FIELD RELATIONSHIPS 


The Tal Series of the area consists of a great variety of argillaceous and 
arenaceous rocks. Often they create illusion with regard to their position 
owing to repetition of sediments of similar nature. All the rock types pass on 
or grade into one another and thereby create much difficulty in demarcating 
sharply the contact planes of the different lithic units. The author has, 
however, included each transition bed into the younger formation, because the 
transition beds indicate a beginning of a change in lithology from the older 
sediments. 


Fic. 1. Photo illustrating the intricate system of bifurcating valleys—a characteristic 
feature of the hills where Upper Tals are exposed. Photo taken from the hill near Kuwal 
village and looking due south. 


SUCCESSION 


Auden (1934) has proposed a twofold division of the Tal Series of this 
locality. The Lower Tal formation of Auden, however, can be further sub- 
divided lithologically in view of the presence of different and distinct rock 
types. The base of the Lower Tal consists of soft black shales with chert 
but the upper portion is made up of hard and tough siltstones. Although 
there is no sharp line of junction between the two, and the shales gradually 
change to siltstones, there seems no justification in combining them 
together, particularly when they strikingly differ in their physical 
properties and composition, which are more elaborately discussed ahead. 
Therefore, the present author has further subdivided the “Lower Tal’ of Auden 
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into Lower Tal (in restricted sense) and Middle Tal, leaving Auden’s “Upper 
Tal’ as such. 

The whole Tal Series rests somewhat unconformably over the Krol C 
unit of the Krol Series, as evidenced by the difference of dip direction and 
the presence of a thin conglomerate bed between the two series. This 
unconformity, however, is not persistent throughout and is found only in the 
western side of the present area. The Krol C unit in the area consists of 
sandstones at the bottom and bluish limestones at the top. 

After a careful study the following stratigraphic sequence has been 
established : 


Calcareous sandstone. 
Protoquartzite and orthoquartzites 
with intercalations of red and green 
Upper Tal shales. 

Arkose and thin shale beds. 

Interbedded siltstone and siliceous 
bands of varying thickness. 

Thick series of dark coloured siltstone 
passing upwards to light coloured 
siltstone, with ripple marks. 

Banded, black, carbonaceous siltstone 
passing upwards to dark coloured 
siltstone. 


| Massive black, carbonaceous silty 


Tal Series (Jurassic) 
Middle Tal 


shales gradually merging upwards into 
banded, black, carbonaceous siltstones. 

Soft, black shales interbedded with 
very thin and irregular chert bed. 


Lower Tal 


—— —_ siicht unconton(y—<————— 
Krol Series 
(Permo-Carboniferous or Triassic ?) 


STRUCTURES 


The structural features of the area can be divided into two broad heads, 
namely ‘tectonic structures’ and ‘sedimentary structures’. Among the former 
are included the fault, folds and joints while in the latter mention may be 
made of bedding, ripple marks, mud cracks and vermicular tubes which may 
be worm borings. 

The Krol C unit is faulted against the other older formations, i.e. Infra 
Krol, Krol A and Krol B. The Tals, however, have not yielded to stress due 
to folding. Very small faults, however, which are of the nature of minor 
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slips along the bedding as well as joint planes could be seen at places. As 
mentioned earlier, the Tals are folded into syncline. Joints are highly deve- 
loped in the Tals, particularly in the Middle Tal formation (Fig. 3). Very 
good dip and strike joints are developed due to which, at times, big rectangular 
blocks are found. Often very good slickenside markings can be seen due to 
slipping along the joint planes in the Upper Tal formation. 


SL CKOTIDHAMAN .' 


[Ju TAL 
EXZIMTAL [TAL SERIES 
EA. TAL 


AN SLIGHT UNCONFORMITY 
EEkroL c 
E=IkroL B¡KROL SERIES 


ESSkroL a 


Fie. 2. Geological map of the area south of the village Kotidhaman, Sirmur District, 
Himachal Pradesh (a part of survey Topo sheet No. 53 F/10). 


As regards the sedimentary structures, nature of the bedding is the most 
striking feature. The most common type is the massive bedding ranging from 
0-5 to 1:5 m in thickness, displayed by the middle and upper ldivisions of 
Tal Series. Cross-stratification or cross-bedding is frequently seen in the 
Upper Tal formation (Fig. 4). These cross-beddings consist of “sets” of beds 
that are similarly shaped and approximately parallel. The sets are of limited 
and generally moderate thickness (0-1 to 0-5 m) and are separated from other 
sets by stratigraphie planes which are, however, comparatively much thinner 
than the cross-bedded units. The cross-beds are of tabular type (Potter and 
Pettijohn 1963) with the individual foresets having dips varying from 20? to 
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25? with respect to the true bedding plane. Their dip direction confirms with 
that of the bedding plane. The foresets have planer basal to slightly tangential 
contact with the underlying stratigraphic plane and are truncated upwards 
against the overlying sets. This feature clearly indicates that the Upper Tal 
sediments are in a normal sequence of deposition and not inverted. Among 
some other less frequent sedimentary structures are the ripple marks in the 
uppermost beds of the Middle Tal and also in Upper Tal formations. These 
are asymmetrical current ripples with broadly rounded crests. Their certain 


Fic. 3. Photo illustrating jointing in Fic. 4. Beautiful large scale cross-bedded 
Middle Tal rocks. Upper Tal sandstones. 


Fie, 5. Photo illustrating mud cracks in Fra. 6. Polished rock surface showing 
red. sandy shale beds of Upper Tal formation. vermicular tube, nodules and turned up edges 
of siliccous bands in the red, sandy shalo beds. 


properties were noticed in the case of Middle T'als only ; for example, they have a 
fairly persistent strike in a direction ENE-WSW, their stoss side direction is 
NNW, having a very small amplitude as compared to wavelength (ripple 
index 17-3). Mud cracks are essentially found in the red, sandy shale beds 
of the Upper Tal formation which indicate that these sediments were exposed 
to sun (Fig. 5). Peculiar vermieular $ubes (which may be worm borings ?), 
whenever seen, are found loealized in these red, sandy shale beds (Fig. 6). 


z 


SI. No. 


10. 


eG 


13. 


14, 


16. 


Modal analysis of the Tal sediments 


TABLE I 


Cement Dod 
Locality and speci- 6 e Mica f fragment 
ano! Quartz % Felspar % % e. OEC. 
E o, chert % 
Near the hill 6903 ft. ; 
57/817 72-88 2-22 1-34 22-42 1-09 
Near the hill 6893 ft.; 
51[819 70:34 1:82 1:90 23-70 138 
South-east of 
Kotidhaman; 57/829 93-70 1:32 0-27 3:66 0-86 
Near the hill 6903 ft.; 
57/831 93:17 nil 0:51 5:23 0-85 
South of the hill 
6893 ft.; 57/841 89-57 4-38 0-25 9:83 2-85 
North-east of the hill 
6893 ft.; 57/843 86:61 1:83 0:68 6:84 3-34 
North of Sakhauli; 
57/845 91-16 2°23 0:52 3:66 176 
West of Pab; 
57/865 82:37 5-90 0-73 54l 4:57 
South-west of Kota; 
57/869 90-98 2-78 0-38 3-72 144 
S.S.W. of Kota; 
57/880 81:59 4:17 0-95 5:61 7:57 
W.S.W. of hill e 
6903 ft.; 57/821 75:66 18-17 0-26 2:52 3:13 
West of the hill 
6903 ft.; 57/822 77:33 16:49 0-28 2:83 
South of Rajaina; 
57/815 60:30 0-40 6-20 32:S0 
North-west of 
Rajaina; 57/859 , 48:90 nil 8-10 42-90 
South-east of Kuwali; 
57/876 50:70 0:20 4-70 43-60 


Heavy 
minerals 


% 


0-04 


0-75 


0-19 


0-24 


0-13 


0-26 


Dark green; medium 
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PETROGRAPHY 


For the sake of clarity the petrography of each formation is described 
below separately : 


Lower Tal 


The Lower Tal formation consists of chert and black shales. In the 
lower portion soft black shales are found interbedded with thin irregular 
chert bands, while in the upper portion massive, black, carbonaceous, silty 
shales are found. Both the shales and associated cherts are found to be 
phosphatic. The phosphatic shales are softer and soil the fingers and, there- 
fore, can be easily distinguished from the overlying barren, black shales in 
the field. The phosphatic cherts occur in two ways; firstly, as nodules which 
are rather elliptical in shape and range in size from 1 to 3 cm, and secondly, as 
friable granular variety. The chert beds are very fragmentary and the shales 
are so soft that neither of them could be sliced properly and examined under 
microscope. 

Overlying the cherts are barren, black, carbonaceous shales. Under 
microscope nothing more than silt sized quartz grains with some slender 
flakes of mica in a dark, almost opaque ground of carbonaceous material is 
visible. 


Middle Tal 


The black shales, which constitute the uppermost part of the Lower Tal 
formation, merge very gradually into the banded, black, carbonaceous silt- 
stone, which is a transition zone, and normally it is not easy to separate one 
from the other. However, considering the lithological characters on the whole, 
the Middle Tal sediments can be distinguished clearly from those of the Lower 
Tal. The Middle Tal formation corresponds to the upper part of the Lower 
Tal of Auden (1934), i.e. graywacke. In the field the mode of occurrence of 
these rocks can tempt anybody to name them as graywackes. ‘The rocks are 
hard and massive, virtually structureless, dark gray coloured and break with 
a sub-conchoidal fracture. All these megascopic characters can be seen in 
any typical graywacke but the mineralogical composition, as described below, 
does not fulfil the essential conditions to call it graywacke (Table I). The 
term graywacke was originally used by Auden with the following explanation: 
‘... No ferromagnesian minerals like augite or hornblende has been found, 
such as are probably required by the definition of these rocks, but the abun- 
dance of chlorite suggests the formation from minerals derived from basic 
rocks.’ But for such assumption evidences such as relies of augite or horn- 
blende in chlorite or pseudomorph of chlorite after the ferromagnesian minerals, 
etc., are necessary which was never seen in a single instance. Moreover, 
felspar is rarely seen while rock fragments are totally absent (Kanjilal 1966), 
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both of which are essential for a rock to be called graywacke (Pettijohn 1957). 
Due to this mineralogical paradox, size phenomena are used by the present 
author in giving a new name to these rocks. Since the maximum size of 
about 80 per cent of the detritus is never more than 0:06 mm, the term ‘silt- 
stone” seems to be the most suitable one. 

These siltstones are massive, well jointed, dark coloured, micaceous and 
fine grained. The lower part is black, carbonaceous siltstone marked by fine 
bandings of lighter shades of gray but of similar composition. They merge 
upwards into dark coloured siltstones, the upper part of which gradually 
becomes light coloured and more fine grained. There is a characteristic 
phenomenon invariably true for the entire Middle Tal formation of the area 
and that is a decrease in the amount of carbonaceous matter and size of detrital 
quartz grains, together with an increase in argillaceous matter, when followed 
upwards. The rocks also respond to dilute HCl, indicating the presence of 
carbonate. 

Under mieroscope can be seen angular to sub-angular quartz grains, small 
flakes of micas (generally ferrimuscovite) showing a parallel to sub-parallel 
orientation which may be due to slight reworking in the basin of deposition 
(Fig. 7), some chlorite which is undoubtedly detrital and not of secondary 


Fic. 7. Dark coloured siltstone, showing somo parallel arrangement of micas. 


à Crossed 
nicols; 38:88 x . 


origin as advocated by Auden, and a very high percentage of argillaceous 
matrix. The matrix is a homogeneous mixture of very finely divided pulp of 
chlorite together with some clay minerals, sericite, mica and extremely fine 
grained quartz. Carbonate dust is invariably present in the matrix. In the 
lower horizon the matrix is flooded with carbonaceous material. 
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Only a few grains of plagioclase felspar are seen scattered in the matrix. 
Small erystals of pyrite and rarely rounded tourmaline, zircon and iron oxide 
are seen. 


Upper Tal 


The Upper Tal formation consists of different types of sandstones. They 
show a wide range of colours, e.g. white, shades of green, yellow, brown and 
grey. Pinkish and reddish beds are at places found interbedded with them 
(Kanjilal 1966). They also have a wide range of size grades, ranging from 
silt size to pebble size. The upper portion of this division shows some alter- 
nate coarse- and fine-grained bands with intercalations of silty, red and green 
shales. 

The microscopic study reveals that the sediments of the Upper Tals 
should be broadly named as sandstone only. The author wonders if the 
term 'quartzite', used at several places by Auden for different beds in the 
Upper Tal formation, stands for a metamorphic rock. (Please see the sequence 
given by Auden in his memoir, ‘Geology of the Krol Belt”.) It is certainly 
not possible to imagine a metamorphic rock sandwiched between two sedi- 
mentary beds and showing a full conformity among them. The rounded 
grains, overgrowth on the detrital grains and the presence of cement verify 
that they are still a sedimentary rock and not metamorphosed to quartzite. 
On the basis of modal analysis, four different rock groups have been made, 
namely arkose, protoquartzite, orthoquartzite and calcareous sandstone. 
The terms ‘protoquartzite’ and ‘orthoquartzite’ may not be acceptable to all, 
but the author prefers them as indicating a particular type of sandstone. 
Krynine (1948) defined orthoquartzite as “pure sandstone' with 90 per cent or 
more detrital quartz and protoquartzites as those with 75-90 per cent of 
quartz and cement below 15 per cent. This mineralogical definition is follow- 
ed throughout in this paper. Besides those sandstones with more than 12 
per cent felspar are considered as arkose (Pettijohn 1954) and those with 
less than 75 per cent of quartz as only sandstone (e.g. calcareous sandstone). 

Roundness of all the above-mentioned four groups has been measured 
by Powers’ method (1953). Sphericity was measured by drawing outlines of 
grains as seen in thin sections with the help of a camera lucida. The areas of 
different grains were measured along with the largest diameters of the grains. 
The sphericity for each grain was calculated by the following formula (Krum- 
bein and Pettijohn 1938): 


$1 (sphericity) = x 


where 
dc — diameter of circle of same area as magnified grain projection 


= Vide 


vti pyra * E K 


» 


o wh A is the magnified area of the grain and De the diameter of the smallest 
circle circumscribing grain projection. - ; 

. Actual roundness and sphericity data for four representative thin sections 
representing all the types of Upper Tal rocks are given in Tables II and IIT. 


TABLE II 
Roundness data 


Specimen No. Numerical mean Roundness/size 


~ 


Sb Ney and name roundness Borm ofgrain relationships [ 
1.  57[819 0:32 Sub-angular Smaller grains aro moro 
> Cale. sandstone angular. 
2. 57/831 0:39 Sub-rounded Both smaller and larger 
Orthoquartzito grains are equally 
rounded. 
3. 57/843 0-41 Sub-rounded Roundness varies in the 
Protoquartzite grains of same sizo. 
E 4. 57/822 0:37 Sub-rounded Both smaller and larger 
Arkose grains are equally 
: rounded. 
—————————— —Re— 2 EA A 
TABLE III 
Sphericity data 
——————— 
96 of grains in specimen No. 
Sphericity 2 
grado 57/819 57/831 57/843 57/822 
A IA eS i 
1:0 -0:95 = 7 3 z 
0:95-0:90 4 10 5 2 
0:90-0:85 6 18 12 6 
0:85-0:80 30 25 18 32 
0-80-0:75 35 15 26 26 
0:75-0:70 20 12 18 14 
0-70-0-65 5 8 12 10 
0-65-0-60 = 5 6 7 
060-055 — — = » 3 
ie S eS 
Petrography of each rock group is dealt below one by one in detail: 
(a) Arkose—After a few tens of metres of interbedding, there is an arkose 
E" i g . B y 
bed of varying thickness.’ The rock is very coarse grained, hard and compact 


and pinkish in colour. The major constituents are quartz and pink felspar. 
- The grains seem to be sub-rounded with naked eye and range in size from 

"5 mm. _ Under microscope the rock consists of Sub-rounded to rounded | 
rS  felspar grains. The qu grains are mostly inclusion 


A 


| eg ia ax 
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free and have straight to slightly curved contacts. More than 90 per cent of 
the detrital quartz show undulatory extinction and rest of them are poly- 
erystalline (8 per cent) and non-undulatory type (2 per cent) (Conolly 1965). 
It is seen that only the highly rounded grains show a non-undulatory extinc- 
tion. On some detrital grains of quartz, silica overgrowths in the same optical 
orientation are seen. Two types of felspars are seen, microcline and un- 
twined felspar, whose species could not be determined due to its weathered 
state. The latter type is much less in amount than the former one. Often 
rhombie overgrowths are seen on the felspar grains (Fig. 8). Among other 


Fie. 8. Photomicrograph showing an authigenic felspar (F) grain at the centre. Two 
authigenic quartz (Q) grains and cement are also seen. Crossed nicols; 155:55 x. 


constituents are insignificant amounts of muscovite, granitie rock fragments 
and cement which is primarily a clay-chlorite mixture. 

(b) Protoquartzite and Orthoquartzite—The protoquartzites and ortho- 
quartzites are of different colour and grain size. Comparatively fine grained 
rocks are more micaceous. Thin shale partings are always present after 
about 30 to 40 m of these rocks. In the upper part these rocks contain red 
coloured sandy shale beds with frequent mud cracks and vermicular tubes. 
Under microscope this group of rocks consists of medium to coarse grained 
sediments (0:31 to 2:00 mm) with individual grains having a sub-angular to 
rounded form. Quartz is the most dominant mineral, most of which is clear 
and inclusion free. However, few grains show inclusions of zircon, rutile, 
muscovite, quartz and iron oxide which are either scattered at random or 
arranged in a linear fashion. The grains show different types of contact 
(from curved to straight) due to diagenetic effects. Silica overgrowth, which 
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is always in optical continuity with the nucleus, is very common (Fig. 9). | 
As in arkose, maximum of the detrital quartz grains (90-97 per cent) are j 
of undulatory type. Among rest of them, 0-4 per cent are non-undulatory 1 


i 

i 

1 

5 

i 

Y 

D: 

Hi 

y 

3 

Fic. 9. Orthoquartzite with alternate fine Fic. 10. A typical protoquartzite from the i 

and coarse bands, Also showing rounded area showing a poly-crystalline quartz (Pq) at i 

grains and secondary overgrowth on quartz. the centre, an undulatory quartz (Uq) and a i 
Some quartz grains show a linear inclusion of granitic rock-fragment (Gr) in the middle left 
iron oxide. Crossed nicols; 33:33 x. hand side of the photo. Crossed nicols; 

33.33 x. 
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Fic. 11. Secondary silica cemont (S) Iie. 12. Anotber protoquartzite showing 1 
(derived from primary clay-chlorite cement-C) clay-chlorite cement. Note how the cement 
forming a wedgo shaped overgrowth over has corroded the quartz grains. Ordinary j 
quartz grains. Crossed nicols; 33:33 x. light; 33-33 x. x | 
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and 0-10 per cent of polyerystalline types (Fig. 10). The second most pro- 
minent mineral is felspar which is chiefly microcline. Orthoclase is altogether 
absent and plagiclase members are rare. Some mierocline grains show in- 
clusions of accicular muscovite. Mica is present in very small amounts except 
in a few specimens. Muscovite is the chief member which is actually an 
alteration product of biotite by the removal of iron. 

The matrix and cementing material of these rocks are very interesting. 
They are always less than 15 per cent (Kanjilal 1966). Three types of cement 
have been noted: (i) primary clay-chlorite cement, (ii) secondary silica cement 
and (iii) secondary ferruginous cement. 

The first group consists of an extremely fine grained mixture of clay and 
chlorite which appears to be primary and deposited in the basin together 
with other detritus. This type of cement dominates over the other two types. 
It is, however, difficult to say that which one of the two (i.e. clay or chlorite) 
is greater in amount. The silica cement appears here to be secondary and 
derived from the basin of sedimentation itself by diagenetic changes. Two 
mechanisms may have operated in the derivation of this cement: 

(1) “Silica leached out from the clay minerals due to diagenetic changes in 
the basin of deposition.' (Towe 1962). This is evidenced by silica aggregates 
near or around the clay-chlorite cement (Fig. 11). Where leaching is maxi- 
mum, relies of clay-chlorite are found within the silica cement. 

(2) Silica removed in solution from the point of contacts between two or 
more than two quartz grains and from the detrital quartz grain itself attacked 
by the primary cement (Fig. 12). 

The ferruginous cement has developed only locally. The probable 
source may be the iron released from biotite during the process of alteration to 
ferrimuscovite. 

Among other constituents are the rock fragments af schists, shales, 
granite and chert. These occur as rounded to sub-rounded grains. Several 
types of heavy minerals including both opaque (like pyrite, iron oxides and 
ilmenite) and non-opaque (like tourmaline, zircon and rutile, some apatite, 
garnet and chlorite) varieties are seen. 

(c) Calcareous sandstone—This is the topmost member of the Upper Tal 
formation as occurring in the area. It is brownish white in colour, hard and 
medium grained. The weathered surface is deeply stained with limonitie 
matter. The term calcareous sandstone is given to ib because the quartz 
content of this rock type is less than 75 per cent and the cement is chiefly 
calcareous. Under microscope the texture is medium grained, with sub- 
angular to sub-rounded grains of quartz and felspar, cemented in a calcareous 
matrix. The quartz grains are mostly of 0-125 to 0-250 mm size. Besides 
the calcareous cement which is definitely primary, some limonitie matter also 
oceurs in the matrix. 
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ENVIRONMENTAL CONDITION Or DEPOSITION 


The Lower Tal formation consists of interbedded black shales and cherts 
which are about 200 ft. thick. This may have been deposited in an euxinic 
condition (Pettijohn 1957). The association of this black shale with car- 
bonate-sandstone (Krol C) rocks below and the immature rock, siltstone 
(Auden's greywacke) above it, records a transition between the open sea 
deposition and craton derived sands and the influx of sediments from a rising 
island arc. The initial rise of the future are results in barriers which lead to 
partial isolation of the locus of deposition and restricted circulation conditions 
conductive to stagnation of bottom waters and deposition of reduzates and 
related sediments (James 1954). 

When the basin of deposition was almost filled up, a delto-continental 
condition was probably reached and the mode of deposition of sediments was 
attended by a change in provenance, during which time the interbedded 
siltstone and siliceous bands may have been deposited. Interbedded red shale 
beds of the Upper Tal formation indicate a strongly reducing condition. 
Mud cracks and the vermicular tubes in them also point to a land or near- 
land condition under hot and arid climate. Sedimentary parameters of the 
Upper Tal formation indicate a beach environment. 
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In this note the motion of light rays in the stationary gravitational field of a 
rotating star is discussed. It has been shown that rotation drastically 
modifies the behaviour of the ray. However, for stars with small angular 
velocity the escape of ray is similar to a non-rotating case. 


1. INTRODUCTION 

Recently, interest in the problem of gravitational collapse has been 
revived by the discovery of quasi-stellar radio sources and the possibility of 
ascribing their enormous energy output to gravitational origin. Under 
continued gravitational collapse the energy of the star is given away in the 
form of electromagnetic radiation. It is interesting to study the effect of the 
gravitational field of the collapsing matter on the escaping radiation. The 
effect of the gravitational field of a stationary spherically symmetric body on 
escaping radiation has been extensively discussed by Hagihara (1931), Darwin 
(1959, 1961) and Synge (1966). In this note we wish to discuss the escape of 
photons from a spherically symmetric massive body which is rotating about 
an axis passing through the centre of the body. 

The rigorous solution of Hinstein’s field equations corresponding to the 
above problem is, unfortunately, not known and, therefore, we make use 
of an approximation method. We assume that the angular velocity of the 
rotating body is small and consider the effect of rotation as a small per- 
turbation away from the static Schwarzschild metric. A general perturbation 
away from the Schwarzschild metric will give rise to various modes of gra- 
vitational radiation (Brill and Hartle 1964). Here we confine ourselves to the 
case when resulting metric is stationary. 


2. THE METRIC AND THE PATH OF THE Rays 


Let the axis of rotation of the body be taken as z-axis of the coordinate 
system. Then we may take the line element of the external field of the 
rotating star as (Kerr 1963) 


A 2m|-' PS 2m 
ds? = | 1— aa dr24+12(102+sin? 8 d$2)— | 1— = di? 
Oma meu ate 2m|-' , 
neri sin? 0 d$ fas (1—22) den xd 5c s BY 
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Here m is the mass and a is a constant such that ma is the angular momentum 
of the body. The form of the line element is suggested from the linearized 
theory. It is seen that (1) differs from Schwarzschild metric only in the 
addition of the last term. (The metric given in eqn. (1) is, however, an 
approximation of Kerr's metric.) 

It is convenient to introduce the dimensionless coordinates 


t 
= — Tam 
P ES 2m 


so that eqn. (1) reads 
ds? = 4m?[(1—p71)71 dp?+p2(d02+sin2 0 d$2)— (1—p71) dz? 
+Kp7! sin? 0 dd(dr--(1—p71)-1dp)] .. t T a5 (2) 
where K = afm is a characteristic constant of the rotating body and for a 
given body of spherical shape depends on its angular velocity. The above 
equation shows that for all stars without rotation the behaviour of the ray of 
light may be discussed in a single argument because m appears as a common 
factor if the last term were absent. Thus we conclude that in case of rotating 
stars the path of light rays will be determined by their mass as well as angular 
velocity. 
In order to discuss the path of a given light ray we choose a coordinate 
system in which the ray is emitted initially in the plane @= 7/2. Now it has 
been shown in the Appendix that it will continue to move in the same plane 


in the approximation we are using here. Then we obtain the equations of 
motion of the null geodesic as 


d$ Ap"? dr 


du > du 


a (DY E on (dr\? 
(=p) (2) +p? (3) —(l—p7}) (5) 


_, d$ Ts d 
CKD Qr ya SO 
tKp du an 9-5) PED 


E) 
a 


where A and B are constants depending upon the initial conditions and u is 
an affine parameter along the ray. The constant B gives the inverse of the 
distance of the light path from the point of origin of the coordinate system 
(Bergmann 1960). In order to avoid Schwarzschild singularity we have to 
Yo 


discuss the above differential equations for Po <p<0 w BU 
2m ` 


here 1 «py = 


Suppose now that we follow a particular ray emanating from the surface of 
the star. "Then, from (3), the differential equation of the ray is obtained as 


d o = r2 
(Eus Jj = p ES pac n E 


(4) 
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Thus for a star of given gravitational intensity the path of the ray will be 
given by the above equation. 


3. [DISCUSSION AND CONCLUSION 


To discuss the escape of rays in general we draw curves with equation 


KDE 
pr plc cum 

p? p3 4pS 
with B and p as coordinates for different values of K. A point in this diagram 
represents a point on a ray. In Fig. 1 we have drawn curves for K = 0-5, 


1, 1:5, 2. For the sake of comparison we have also drawn the corresponding 


0 PR . (6 


1 2 3 4 


P 

Tic. 1. Escape of photons from a rotating star (curve with K = 0 corresponds to a non- 
rotating star and curves for K = 0:5, 1:0, 1:5 and 2:0 correspond to stars with increasing angular 
velocity). 


curve for a non-rotating star. Since the right-hand side of eqn. (4) cannot 
be negative, no representative point can lie below the curve for a given value of 
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K. Further since p lies in the range 1 <p < © no point can lie to the left | 
of the line p = 1 (shown shaded in Fig. 1). 

The interesting point is that for small values of K the curves have maxima, 
similar to the case of non-rotating stars, which shift towards left for increas- 
ing values of K. But for K > L5 (approx.) there are no maxima. This 
leads to the conclusion that as long as the angular velocity of rotation is 
small, the escape of the ray from a rotating star is similar to that of a stationary Y 
star. "Thus, three typical points, say A, B and C, may be seen in the diagram 
for a particular curve, say for K = 0:5, having a maxima. In each case, as 
we move along an outgoing ray, the representative point moves horizontally, 
starting to the right and continuing unless it runs into the curve. For those ] 
points, e.g. A, for which the horizontal line runs into the curve the ray has an | 
apse. But if the representative point lies above or to the right of the maxi- 
mum (e.g. B or ©), it may be moved horizontally to infinity without meeting 
the curve. It is clear from the diagram that the curvilinear triangle A, As 
A; is the domain of recapture. : 

In case of stars with large angular velocity we have a large value of K 
and corresponding curves have no maxima. It means that none of the rays 
has an apse and a ray once emitted from the surface of the star escapes to 
infinity. 

The above considerations apply to all rays, whether incoming or out- 
going. If we wish to discuss the behaviour of an incoming ray, the representa- 
tive point will move horizontally from right to left. All the previous argu- 
ments equally apply to this case as well. 

Thus, we see that the escape of a ray emitted from the surface of a star is 
sufficiently modified due to rotation. But there are certain limitations to 
the conclusions which we have drawn here. Firstly, we have used an approxi- 
mate form of Kerr's exact metric (1963) retaining the first power of the cons- 
tant a. In doing so we intended the metric to be valid to first order in the 
angular momentum parameter. But to this order we obtain the Schwarzschild 
metric plus a term in dødt. Although many authors have used this procedure 
(Brill and Cohen 1966, and Landau and Lifschitz 1962) but there is no 
reasonable justification for adding a term proportional to dødt. The con- 
clusions drawn here are, therefore, not rigorous.* Secondly, the behaviour of 

_ the rays in only equatorial plane is discussed.+ Thirdly, the constant K in 

egn. (2) will generally have lower values than the one we have assumed here. 

— Nevertheless, the conclusions drawn in this note are interesting. It is 
ected that more exact treatment of this problem will lead to similar 


roof: For the full discussion of Kerr metric, see Boyer, R. H., and 
hys. 8, 265 (1967). 
all the rays should be discussed. 
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APPENDIX 


The equation of motion of the ray is obtained from the geodesic equation 


dun u da* dak 
du? «B du du 


When the above equation is specialized for the metric (2) for a? = 0, 33 = $ 
and a4 = ( we obtain the respective equations for the plane 0 = 7/2 


d20 
X390 eo us VS NEED) 
EMI 


du? du du 2mp \du 
p p 


a (dry? 1 2 a zn 
Ee a (ali a E 
4mp® (5) ar 3p | isset) | 2mp* (E 


a Hoy? | dr dp a? dr d$ 


dp d$ _ 


—p. NY 429 
ALMA MM uc 


dup 
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dèr |l dr dp a dp d$ 


(1—7?) het p? du du 2mp? du du 


a (dry? o e hee] 
sa (E) t-e T) xs Ip} du 


a? E a 5 Ee e dp 


7 4m2p3 \du a du du 


4m?p? | p 


E Qu CR eae do cem (43) 
8m3p5 du du  8m3p^ du du 


Equation (41) shows that a ray initially moving in the plane 0 — z/2 will 
continue to move in the same plane. The other two equations cannot be 
integrated as such and we, therefore, make use of the approximation referred 
to earlier. We expand different terms in decreasing powers of p and retain 
terms of the order O(p73). Thus, we obtain 


d 2 dp dġ a (2) =o id 3 Eo 


du? "p du du ` 2mo? du 


and 
dr 1 dz dp a de dọ 
1=p7! == : = c 
Cay dui pi du du 2mp? du du 


Equations (44) and (45) are completely integrable but for the presence of 
a a dp dj 
terms —— >. 
2mp? 2mp? du du 
consider the effect of rotation as a small perturbation away from the stationary 
body. The corresponding equations for spherically symmetric stationary body 
are completely integrable and if we compute the order of magnitude of 
2 

me. (2) ein) cs ee Ce a 
2mp? \du 2mp? du du 
we are interested in. Therefore, we are justified in neglecting even these 
terms. Thus we obtain eqn. (3). 


(dp/du)? and — But as we have already remarked, we 


this is found to be much smaller quantity than 


ON IMPROVING AN APPROXIMATE FACTOR OF A 
POLYNOMIAL WITH REAL COEFFICIENTS BY 
AN ITERATIVE PROCESS 


by R. D. Brarcava and R. N. SoRAL,* Indian Institute of 
T'echnology, Bombay 76 


(Communicated by R. S. Varma, F.N.I.) 


(Received 1 June 1968) 


Bairstows’ (1914) method of improving an approximate quadratic factor of a 
polynomial of degree n > 2 with real coefficients is generalized to improve an 
approximate factor of an arbitrary degree m < n. The treatment is an algebraic 
one but the numbers can be easily substituted. As is easily seen the method 
can have several applications, The algorithm for such a method is given 
and its convergence is proved in this paper. 


$1. Consider the polynomial 


Pla) = agv^ Hartt -1 doaP 2+ 26. An 1-0 


of degree n with real coefficients do, %1, ... , @ and let an approximate factor 
of it be 
P(x) = ppal H pox-24- ... pm qt Pm 


where m < n. Then the problem is to improve the coefficients pj, ps. . -s Dm 
of P(x) so that the modified coefficients are better approximation to the factor 
of F(x) than the original ones. In other words, if the identity 

P(e) = (Mpeni Fo) ERA Ram Ri (1) 


holds, where 
Q(x) = borr mt byan MAL ... Hbn-m 


and the coefficients bo, bi, ... , bn-m, Ry, Ro, ... , Rm depend only upon 
fp Po, ---; Pm, then the problem is to find pj, Pa, ... , Pm, so that 


Rip pes... Pm) 50; t=1,2,...,m n m 


For this we use Newton's interative process. That is, we find the small 
corrections AP], Ape, ..., Apm, SO that 


Ri(pi+ AD Pot ADs,--- > PmF Am) = 0; GOA SUN nme (3) 
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Expanding eqns. (3) in Taylor series and supposing that the squares and 4 
products of Ap; are negligible as compared to Api, we truncate after the j 
first order terms and get a system of » linear simultaneous equations 


7 m 
| QR. : 
Ri(pys pa» +++» Pm)+ » Ap. =0; t=1,2,...,m ss (E) 
£ t=1 
in m unknowns Agy Apo «++, Apm. When this system has been solved, 


the procedure is repeated with the corrected values for pi, Pa, ---, pm. The 
problem, however, is to find the values of 9R:/0pi. For this, we proceed as 
follows. We use the identity (1) to obtain 


m 
by = 4 Y tan E z a (O) 
j=l 
where 
k=0, 1, 2; e, NM, AA San) 
and 


m 


Ri = Un-m+i— > Din -m4i-k 


/ is 


i-i 
= bs maid P DkÜn-m4A-k 50 auo 00 (6) 
hal 
T: where i = 1, 2,..., m and 0n m1; Un-m+2, - » -> bn are also defined by equ. (5) 
> for k = n—m--l,...,m. 
Differentiating eqn. (5) with respect to p; we obtain 
br _ S Pbr- 
27 == (i 2 8p, Dj 5d ad Da (7) 


It may be noted that a are constants and, therefore, their derivatives are 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


zero. Putting = = —Cg-1, it is obvious that 
i 
Obr = Obro E > Obi cm Se 
Opn) Op; °° Pm T 


for all &. Inserting cx in (7) we find 


Ck = bi > 303 Ls os s oo E) n 


ES 
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now differentiate eqn. (6) with respect to p: and write Ry for 0.R¿/0p; to obtain 


1-1 
aa ob 
Ru = n-m E 5 (v Ron URN Bes ee 


Opt =A 
m 
DkCn- m4i-t-ks ifti 
Š- k=t 
= Te 00 sc (9) 
= X DkOn-m4i-t-k (po = 1) if t>7 
k=0 


Thus the values of Ru are known and the system (4) can, therefore, be solved 
for Aji, Aps,.--,; Apm. The whole procedure is summarized in the following 
algorithm. 

Algorithn—Given an approximate factor 


P(x) = a+ pam... +m 
of the polynomial 
P(x) = aga®+aya"-14 ... Han 
determine for each y=0, 1, 2,... a sequence of coefficients 2m M= 12, 
., m) of P(x) as follows. Determine 
(i) the sequence (00?) from 


m 


b” = a,— » jun p” 


j=1 
where k = 0, 1, 2,..., n and 07 =o” =... =b” =0. 


(11) the sequence (otn from 


m 


Y (Y ORPA Y 
=L Bew ) 


where k = 1,2,...,2—m and em -—0 
(iii) the sequence (RP?) from 
0) p” (po? 
R; b, cup > Pr ona =mti-k 


kel 
where i = 1, 2,..., M. 


(iv) the sequence (RU?) from 


MU : > 
> Pr Ch —m+i-t-k if i< 
koi 


me (Y) : ` 
o Pr C m+i-t—k ifi>t 


k=0 


es 
Ri a 


where i, t = 1, 2,..., m. Then 
4 


TT pp T 
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+1) 


( (7) 
(v) the set (pf); k=1,2,...,m) is obtained by p? = pP Ap?) where 


(^2; D). ki 9... mp i is determined by 
m 
> RY ap” ec (Sil sue 
t=1 
$2. To prove the convergence, we give below the conditions sufficient to im- 
ply that there exists a neighbourhood around a real solution aot (Pr, Do +. 25m) 
of the system (2) such that for an initial choice p = (p, p Y > D.) in that 


neighbourhood, the above process converges to the solution. 
Let D(p) denote the determinant of the coefficient matrix of system (4), 


i.e. 
Ry Ris ose Rim 
Ror Rog geg Rom 
Dip) = 
Rmı Ring ooid Rmm | 


and let (81, S2,..., Sm) =s be a solution of the system (2). If the function 
Ri(p) have continuous first partial derivatives in the neighbourhood of the 
solution p — s, and if the determinant D(p) is non-vanishing in the neigh- 
bourhood of this solution, then Newton’s process of solving the system (2) 
converges (Householder 1953). Also if D(s) z 0, then since D(p) is continuous, 
it follows that D(p) 0 for all points p sufficiently close to s. Thus in order 
to show the uniqueness of the solution of the system (4) and the convergence of 
the above method we have to show that D(s) 4 0. To prove this for general 
case seems to be difficult but if the zeros Qj, Vo, ..., Um Of the factor 
omt siwi .. . Sm of F(x) are simple and distinct the cool can be proved as 
follows. tao F(x) is independent of D» f» ..., Pm, differentiating the 
identity (1) with respect to Pi, 92; : - - , Pm respectively, we obtain 


0 = am-10(5) + P(e) Feta Ry ana, y soo Bon 
1 


0 
0 = qm-20 (a) +P (0) seem LI» +2M2 Roo + Re Ryo 


0 = Q(u)+P(a) XO Fn Mos bo" Ran eo ERinm 
Setting p; = s; and «= a (6 1, 2,.:., m) we obtain m systems 
Gye Rate” Rp co ER = qr (CSU, 0:00 20) 
v Rott” Rat.. a ES = =x, Es (0 = ll. 
D . . . . (10) 


„m-l SUL . 
E Rim+t;' Boyd... dE = — Qi (1 = 1, 2,..., m) 
4n 
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each of m linear simultaneous equations in m unknowns with the non-vanishing 


determinant 

m 1 12 

E v, l 
m-1 m-2 

La v. b l 

AS 

m 1-9 

„m-l m-2 1 


= IH (uy), 1<<i<m 


i>j 
Æ 0 (because the zeros v, Ta, ... , Un 
of F(x) are distinct) 
of the coefficient matrix A. Here Qi = Q(a:) and the derivatives Ry are 
taken at p — s. Now we shall show that the solution vectors of the systems 
in system (10) are linearly independent. For this we have to show that 


(Ry, Ro... Rmi) +0 (Bio, Rev, + Rmo) F : - - Faml(Rim Rem ..., Rmm) =0 


(11) 
is true only for a, = 4$ = ... = Qm = 0. Let the identity (11) be true of 
any constants @, (la, ..., dm. Then 

MGRi+0Rj9+ Ws +amRim = 0 \ 
Roy 4-0 RH oos HamERem = 0 | (12) 


. ó o c || 
Q4 Rin 4 Go Bg d- ++. HamEmm = 0 | 
From (10) and (12) we obtain 
"TQ jay 4-27 Qe 4- so OS) (Uem ls aos a0) de (018). 
But Qi = Q(x) Æ 0, because the zeros vı, e; ... , Um of P(x) are simple zeros 
of F(x). 
Dividing (13) by Qı we obtain 
aaa a+ Et = 0 (US LA) .. (14) 
Since the determinant | 4 | of the coefficient matrix of the above system (14) is 
non-zero, this system bas only trivial solution, NUTS een sro UA 
Hence the solution vectors (Ri, oj .... Rmi) (Riz, Ro, ..., Rmg),..., 
(Rim, Rom, «+ + > Rmm) are linearly independent from which it follows that 
Ry Ry ... Em 
Ro, Ros ... Rom 


D(s) = 0 


Rmx Rng SiGe Rim 


Ds 
— BHARGAVA & SORAL: IMPROVING A FACTOR OF A POLYNOMIAL 
"s i 


Thus for all p sufficiently close to s the system (4) has a unique solution 
1, APz :--, APm and the point (pit+Ap1, pot Apos - «Dm Am) is 
more close to s than p. i 
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ON INITIAL DEVELOPMENT OF AXISYMMETRIC WAVES IN 
FLUIDS OF FINITE DEPTH 


by L. DEBNATH,” Department of Applied Mathematics and Theoretical 
Physics, University of Cambridge, Cambridge 


(Communicated by P. L. Bhatnagar, F.N.I.) 


(Received 20 April 1968) 


An initial value approach to the axisymmetric wave problem in a fluid of 
limited and very shallow depth is presented. An asymptotic solution of the 
problem has been worked out in a simple and rigorous way by the similar 
method developed by the author (Debnath 1969) in an earlier paper. Finally, 
the principal behaviour of the axisymmetric wave motions has been discussed 
in considerable detail and some physical arguments has been put forward 
in order to explain the strange character of the solution related to the con- 
centrated pressure distribution with the forcing frequency w. 


1. INTRODUCTION 


Debnath (1969) has suggested various reasons in favour of the initial 
value investigation into the general wave problems and explained the relation 
between the steady state and properly posed initial value problem in consider- 
able detail. Usually, it is assumed that the solution of the steady state prob- 
lem obtained by imposing a radiation condition at infinity is generally the 
limit as time T' — co of the solution derived by considering a well-posed initial 
value problem. Nevertheless, there appears to be no generally rigorous proof. 
that this is, in fact, the case and it is usually necessary to examine the limiting 
procedure in each individual configuration. In the same paper, the author 
has introduced a linearized theory of axisymmetric wave motions in a fluid of 
infinite depth due to local disturbances with the forcing frequency c acting 
on the undisturbed free surface of the fluid. The asymptotic solution of 
the initial value problem has been derived in a simple and rigorous manner. 
Finally, the principal features of the wave motions have been described in 
some detail. 

In this connection, reference may be made to the original works of Stoker 
(1957) who considered the unsteady two-dimensional surface wave problem 
due to a periodic concentrated pressure distribution. With due investigation, 
he suggested that steady state wave problem with the Sommerfeld radiation 
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condition at infinity should be derivable as the limit of the corresponding 
initial value problem. l 

The corresponding axisymmetric wave problem in a fluid of limited and 
very shallow depth is also interesting and important on its own merit, hence it 
deserves a special attention. So far, no such work has been reported to the 
current literature at least available to the author. 

The object of this discussion is to present a linearized theory of initial 
development of axisymmetric surface waves in an inviscid, incompressible and 
homogeneous fluid of limited and very shallow depth generated by an arbitrary 
periodic disturbance with the forcing frequency w acting on the undisturbed 
free surface of the fluid. An asymptotic solution of the problem has been 
derived for sufficiently large times by the similar method developed by the 
author (Debnath 1969) in his earlier paper. Finally, the limiting behaviour 
of the solution has been examined in some detail and the main features of 
the wave motions have been predicted. This is followed by some physical 
arguments in favour of the strange behaviour of the solution related to the 
oscillating point-pressure distribution with the forcing frequency w. 


2. MATHEMATICAL FORMULATION OF THE PROBLEM 


We consider the axisymmetric wave problem in an inviscid, incompressible 
and homogeneous fluid of finite depth due to a periodic surface pressure distri- 
bution in the form 

P(R, T) =P(R)“TH(T), R< A) 


2 
=0 RSA (Est 


acting on the free surface Y =0 of the fluid initially at rest, where P(R) 
is an arbitrary function of R, w is the fixed frequency and H(T) is the 
Heaviside unit function of time 7. 

We assume that X-Z plane be the undisturbed horizontal free surface 
and Y-axis be vertical positive upwards. We take cylindrical polar coordinates 
(R, O, Y) with cylindrical symmetry about the Y-axis such that R is equal to 
V X2--Z? from the Y-axis and the origin of coordinate is taken on the free 
surface. 

In view of the fact that the motion starts from rest, there exists a velocity 
potential G(/t, Y; T) which is governed by the Laplace equation 
0?b 1 0D ab 
OR?" ROR OY? — 


—h« Y«0,0« R« c 


Vv?o = 0 T 5 2. (2.2) 


where h denotes the depth of the fluid. 
The condition at the bottom surface is 


dy=0 on Y=—h .. ae ex TEE (273) 
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Neglecting surface tension, the linearized free surface boundary conditions 
are given by 


3 l 
D/+9B8 = — 5 PUT) y=0 e 006 
DEN j^ DX) MEM os 


where suffixes denote partial differentiation and # = E(R, T) represents the 
vertical surface displacement, p the density and g the acceleration due to 
gravity. 

The initial conditions are given by 

B=E=0 at T=0 
] dee 5c so (A 
[Or(K, 0; Dlr- = — PR 0)| (249) 
We further assume that the functions (R, Y; T) and E(R, T) possess the 
Hankel transform with respect to R in the generalized sense. 

Remarks: The problem formulated above can be investigated as a steady 
state just by omitting the initial conditions (2.6). For the steady state analysis, 
one has to impose Sommerfeld's radiation condition at infinity which is, in 
fact, essential in order to derive a solution of physical interest. 


3. SOLUTION OF THE PROBLEM 


It is convenient to introduce dimensionless variables 7, x, Y, 2, a, é, y and 
¢ defined by the following relations: 


w? . Pot 
r, @, Y, 2, @) = —(R, X, Y, Z, A), toT, 7=— HH 
(r, 2, Y, 2, a) zu ) ae 
and 
Pod 
= ogi 


and we introduce a non-dimensional parameter D defined by D = w*h/g. 
Substituting these relations into the basic eqns. (2.2)-(2.6), we obtain 


deed duy = 0 25 bo Sd .. (3.1) 
—Dzy:0,0«r« o 

$20 ony=—D .. ae ox .. (3.2) 

$e» = —p(r.t)) y=0 or Eo . (3.3) 

by — n | 1>0 TEL ee QM 

$2720  att-0 ) es 


d= —p(r,!) att=0.y=0) 


As in our earlier paper (Debnath and Rosenblat in press), we introduce 
the Laplace transform (Tranter 1966) $, 7 of the functions ¢, y respectively 
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with respect to tand then we introduce the Hankel transform $, 7 (Lighthill 
1958) of the functions 9, 7 respectively with respect to r. This enables us to 
simplify eqns. (3.1)-(3.5) and to write down the solutions for ¢ and 7 without 
any difficulty in the form 


E : (I, s) cosh k(y+D) "ma 
Bk, 3 9) = — hD (3.6) 
z 2(k) B(k, s à 
n(k, s) = —* E (3.7) 
where a? is given by the relation 
02 = a2(k) = k tanh kD. T: T .. (3.8) 


Making reference to the inversion formula for the Laplace transform 
combined with its convolution theorem and then to the inverse theorem of 
Hankel transformation, we readily obtain, 


es) 


sy; t) = i —a sin ad —tett) ESO (oy osh Ir D) y, ‘ 
$(r, y; t) = J. (ù cos at—a sin at—iett) k p(I)J (kr) COED (3-0) 
n(r, t) = f. (? sin a£--« cos at—aet) k (I) («2 — 1)-1 p(Ic)J'g(ker) dk. .. (3.10) 
0 


The above integrals (3.9) and (3.10) are the exact solution of the problem. 
But, in general, they cannot be worked out exactly. Hence asymptotic 
methods are utilized to evaluate them for a deeper understanding of the 
transient wave motions. 

In the case of infinite depth of fluids, i.e. when D — co, the integral for 


the velocity potential ¢(r, y; t) as well as the surface elevation n(r, t) has the 
following form: 


00 i, -- — 
$(r, y; t) — I, (i cos Y EL— ^/ k sin Y It — iett)eku lo d(E)Jo(kr)(E— 1)71 dee 


. (3.11) 
n(r, t) = I (VE cos / kt--i sin Y EL— Y Kett) 1312 51) To(Ier) (e —1)3 dk. 
. (3.12) 


These results are exactly identical with those already obtained by Debnath 
(1969). 


4. ASYMPTOTIO ANALYSIS OF THE SOLUTION 


As already mentioned that the function p(r) involved in the non- 
dimensional pressure distribution may be arbitrary, hence it would be 
sufficient for the determination of the principal behaviour of the wave 
motions to take the simple non-dimensional form as 

P(r, t) =eltp(r), r<a 
=0 n PERA 
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with the prescribed form of p(r) as 

(a) p(r) = o(r)/r 
where 8(r) is the Dirac function of distribution and 

(b) p(r) 2 1 

(c) p(r) = exp (—r?/r%) 

(d) p(r) = Jo(Ar) 

(e) p(r) = (a$—72), n>—l. 
We proceed to find out an asymptotic representation for the surface elevation 
n(r, t) in the general form and the above particular cases (a)-(e) can be used in 
the end. 

Let us consider the integral (3.10) for the surface elevation (7, t) which is 
given by l 
nl, t) = | (i sin ala cos ad —actt)(02— 1)71 k a(k) B(k)Jo(kr) dk. 
0 


It may be noticed that this integral has no singularities on (0, co), there- 
fore, its path of integration can be deformed into a path M (say) in the com- 
plex s = k+-ip plane as indicated in Fig. 1. The path M coincides with the 
real axis (0, oo) except that it is diverted round the zeros of the denominator 
of the integrand. We can then write down the integral as a sum of two in- 
tegrals which are now singular at the zeros of («2—1). 


^ 


sik, >, Y 


X Fra. 1. The complex s = k +îp. 
Hence 


nlr, t) = |... (ù sin at+a cos at—celt) s a(s)(x$— 1)71 p(s)Jo(sr) ds .. (4.1) 


where « = a(s) = Vs tanh sD and s = ko is the only real root in (0, co) of the 
equation «2(s) = 1. 
We thus obtain 
aro = AE T us 23 
where J; and J. are given by the integrals in the form 


Jh em | p(s) s a(s) (a2—1)71(0 sin «f--« cos at)Jo(sr) ds 
M 


Img) I, (a2—1) 


where —7 < arg s <7. 
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It may be noticed that the wave integrals I, and I, obtained above are 
very much similar in form obtained earlier (Debnath 1969), therefore, similar 
methods, as before, with a slight modification can be exploited to evaluate 
them without any difficulty. 

Making reference to the method used (Debnath 1969) for the evaluation 
of the steady state wave integral, we readily write down the result for J, in 


the form Y 
ps. x [ze e rm E 2 (L3) 
si rko W'(keo) 
where the function W (k) is given by 
MO esci A (4A) 

Remarks: It can be easily seen that in the case of infinite depth (i.e. 
D — œ), the pole s = kọ = 1 and, consequently, the value of Zə obtained 
above is in perfect agreement with the result of Debnath (1969). 

In order to evaluate the transient integral /1, we first replace the Bessel 
function Jy(sr) by its standard integral formula and next from a simple 
rearrangement, it easily follows that 


where the integrals £, Lo, L3 and L; are given by 


sa(s)pis) , 
L = A e- ttf (8) ds 
IN PECES ES 


fi, = Í Pie) sate etf -(5) ds 


u «(9)— 


- _ | g5se(s) AO Fh 
I= | Bee cos 


Ly = | Bo) sete) e-tf+(st ds 
y «(s 
with 
OE (ue To cos 0. 


It may be observed that these integrals are similar in nature to those already 
obtained earlier by Debnath (1968, 1969), Debnath and Rosenblat (in press). 
Therefore, they can be evaluated exactly in a similar way as before with a 
little modification. However, a few important points may be outlined here 
so that one can easily understand basic ideas behind the asymptotic evaluation 
of these integrals stated above for large values of t. Making reference to the 


works cited above, it is easy to show that the integrals Lg and L, are ofi): 


/ 
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And the stationary (or saddle) point and pole related to the integrals Z, and 
L are respectively at s= kı and s= kg. In fact, s= k, is the root of 
the equation > 

f-(s)=0. 
A similar calculation, as before, leads us to the asymptotic answer to the in- 
tegrals L, and Lə. Then substituting their values to the 0-integral (4.5) and 
making appeal to the method of stationary phase (Jeffreys and Jeffreys 1946) 


for large t in order to evaluate the 0-integral, it turns out that 


_ Blko)iko e (Eo) (eco) + 1) (2) E AS de < ko 
rko 


I,~ 21W*(ko) 
0 A 
$0). E gir —ta(14)) ecl (s 
| + halk) ay 
pp Para) 


(4.6) 
Finally, combining the integrals J, and I, together, we arrive at the asymp- 
totic representation for the surface elevation «(r, t) in the form 
“LND 97 c be RET, » s 
n(r, t) EA DEAN) aP (ko) (E) er roi) ! MES Al 
W'(kg) rko 1—Ha(ko) +1}, ky < ko 


hey (ley) (hy) Bes OE ofi) " 
far dee) e «REI yA guest 


Remarks: The quantity kı = ky renders the solution (4.7) invalid. So, a 
special method is needed to obtain a valid solution for y(», t) near the point 
kı = kg. Since our range of interest is large values of t, it seems to us that 
the derivation of the solution at the above point is not so important here. 


5. ASYMPTOTIC SOLUTION IN THE CASE OF UxLrurrED DEPTH 
In the case of infinite depth of fluid, that is when the parameter D —- co, 
the function «2(s) tends to the limit s. ses the pole s = ky = 1 and the 


saddle or stationary point s = kı = os - (Copson 1965) which are in perfect 
agreement with our previous analysis (Debnath 1969). Therefore, the above 
solution takes the following form: 


OS Nd (AE ANT 
aur) ena, > 
nr, 0) ~ 


r 


0 , t«&2r 


2 re L Ea 
5 COS (e: sin (ei) sou 


vend 
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This result is in perfect agreement with those obtained by Debnath (1969) in 
the particular cases (a)-(d) stated before. In order to compare this result 
(5.1) with those derived in the earlier paper, one has to substitute the value 
for the function p(k). 

We would like to mention here the Hankel transform (k) of the function 
p(r) related to cases (a)-(e). The functions p(k) associated with the cases 
(a)-(e) are respectively given by 

2 ja 
l, z Tit). 9m2e- Em [me = (2) | I, J (Ax) J o(kx) da, 
and 
2nT'(n-4- 1) (ak) - Vat) y (al), 
where Jn(x) are the Bessel function of the first kind of order n, and the finite 
integral is a standard one (Watson 1922; Sneddon 1951). 

Thus the solution of the problem is explicitly known in both finite and 
infinite depth cases. And for infinite depth case, we can refer to our earlier 
paper (Debnath 1969). In the case of limited depth, we intend to summarize 
our conclusion later on. 


6. AxisvMMETRIO Wave MOTIONS IN SHALLOW FLUIDS 
In the other extreme case, that is when D <1, we can replace o?(s) = 
stanh sD by s?D to the first approximation. Consequently, the integrals 
involved in the vertical surface elevation y(r, t) are simpler in nature. And 
a, procedure similar to that advanced in the papers (Debnath 1968, 1969) can 
be adopted to obtain the asymptotic solution for y(r, t) in the following form: 


wf 2m j ( lx dE) r 
Aa) MM 2 > tS —= 
mir, t) ~ (vz v») VD o GD 
0 5 b«& 5 


This solution for ņ(r, t) suggests that the waves are advancing with the phase 
velocity 4/gh which is independent of the wavelength and the amplitude 
of the waves decays like »-!?. Evidently, the waves are non-dispersive in 
character. 


7. DISCUSSION AND CONCLUSIONS 


From the salient and inherent features of the properly posed initial value 
formulation, it follows that the axisymmetric wave problem investigated above 
achieves a unique solution of physical interest. An asymptotic analysis of the 
problem enables us to derive explicitly both the steady state and transient 
wave solutions. Hence solution (4.7) as well as (5.1) for the surface ele- 
vation n(r, f) describes the behaviour of wave formation at all times instead 
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of only in the limit. And the solution for the two extreme limiting cases 
(i) h — co and (ii) & — 0 has also been discussed in some detail. 

The above analysis reveals that the transient component of the asymptotic 
solution associated with cases (b)-(e) does die out very rapidly in the limit 
t—> œ, for fixed values of r. Therefore, the ultimate steady state is attained 
in the limit. In other words, the steady state is set up in the limit provided the 
applied pressure with the forcing frequency w is distributed over a ‘finite’ 
circular region, however small, of the free surface of the fluid. On the other 
hand, classical steady state is not attained when the harmonic pressure distri- 
bution concentrated in a single point of the free surface. These conclusions 
are in excellent agreement with those obtained earlier in connection with the 
axisymmetric wave problem in a fluid of infinite depth and the two-dimensional 
wave problem (Debnath 1968). Although this disagreement is strange to some 
extent, but not really unexpected for mathematical and physical reasons whieh 
have been suggested earlier (1969) and would be stated below. 

In all the above cases (b)-(e), the classical steady state solution of 
special interest is obtained in the form 


qr, t) $ (4.67, 
Vr 


where A is some known constant depending on the actual functional form of 
p(k) involved in the cases under consideration. 

This solution represents the axisymmetric surface waves propagating in 
the medium with the amplitude decaying like r-!7. And in the limit D> co, 
that is h => oo, this solution confirms our previous solution obtained earlier 
(1969). 

The author (1969) has suggested some mathematical arguments to ex- 
plain the strange behaviour of the asymptotic solution (5.1) for the case of the 
pressure distribution concentrated in a single point with the forcing frequency 
w. From the physical point of view, the following argument may be put 
forward to resolve the strange situation. ‘The sudden harmonie point pressure 
with the high frequency content generates a continuous spectrum of circular 
waves of all possible wavelengths and the continuous spectrum is very rich in 
short wavelength content indeed. ‘Thus the sudden onset and the enrichment 
of short wavelength content are both needed if there is to be sufficient energy 
distributed among the slowly varying travelling wave trains so that more and 
more intense accumulations occur at fixed (large) distance r as the time ¿ 
approaches to infinity. That is why the transient component of the 
asymptotie solution (5.1) related to the delta function pressure increases 
like O(t2). So the classical steady state is not attained in the limit. 

On the contrary, the situation is well behaved when the pulsating pressure 
is distributed over a finite region of the free surface and, in particular, 
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when the Gaussian representation of the delta function in the form | 


8(r) = lim l =g Pla? 
a—0 UN T 

is used as a wave-generating mechanism and the limit ¿> co is taken first. 

Therefore, in any physically realistic situation where the surface pressure 
is distributed over any finite region of non-zero dimension, the ultimate steady 
state would be attained in the limit and the attainment of the steady state is | 
directly attributable to the double limit process 4 — 0; ¿> co in the proper | 
order as suggested earlier and a similar but appropriate energy argument as 
outlined above. 
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This paper deals with the morphology, origin and significance of the different 

varieties of tourmaline, occurring in the Lower Kaladgi quartzitic sandstones | 
outeropping around Saundatti and exhibiting authigenie growths. It is | 
concluded that the diagenetic changes resulting in the growth of minerals 
belonging to the complex isomorphous series of the tourmaline group are 
controlled by neither the chemical composition of the detrital nor the host 


INTRODUCTION 
The Lower Kaladgi (Pre-Cambrian) formations of Saundatti (lat. 15° 46”, 
long. 75° 7’), Mysore State, comprise of quartzitic sandstones, shales and lime- 
stones resting unconformably over the basement granites. The general 
classification and their stratigraphical position are as follows: 


tourmaline grain. 
| 


Basalts Deccan Traps 
I 

Sandstones A Badamis 
| Limestone | 
| Shales | 
| Pebbly quartzitic sandstone ; : , | 

i A M . Kaladgis (Pre-Cambrian) | 
Grey quartzitie sandstone ( x 


Pink quartzitie sandstone 

Basal conglomerate 

Granites, dolerites, biotite schist, hornblende- | | 
schist, chlorite schist, epidiorites, cale granu- } Archaeans 
lites, and banded ferruginous-quartzite 


The Kaladgi sediments of Saundatti, being the southernmost fringe of 
the main Kaladgi basin, have a general east-west trend with their dips ranging 
from 7° to 20° N. The quartzitic sandstones show current bedding, ripple 
marks, sun cracks and rain prints. They are fine grained to medium grained, 
being almost made up of quartz (90 to 95 per cent), and can suitably be termed 
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as orthoquartzites, except for the pebbly quartzitic layer which contains 
felspars making it subarkosic. Tourmaline, zircon, rutile, chlorite, sericite 
and hornblende are among the heavy accessories. 

The study of heavy mineral mounts and micro-sections of these quartz- 
ites reveals the presence of tourmalines with distinct authigenic growths. 


DESCRIPTION 


The detrital nature of most of the tourmaline grains is clear. Depending 
on the shape and degree of roundness, they can be termed as very well rounded 
and angular, according to Krynine and Tuttle. The detrital grains and the 
new growths are distinct from each other in dichroism exhibiting varying 
shades of yellow, blue, green, brown, etc. On the basis of dichroism and 
birefringence, the tourmalines under study have been roughly categorized 
into the respective members of the isomorphous series as follows: 


Heavy mounts 


Schorlite usually black, purplish brown, 
dark green greenish brown 
Elbaite pink, pale green, colourless 
greenish blue 
Dravite dark yellowish pale yellow 
brown 
Chromium dark green pale green, 
tourmaline greenish yellow 
Micro-sections (0:03 mm.) Birefringence 
Schorlite yellowish brown, ^ pale to dark 
black, dark yellow, 
green, blue of brownish, pale 
various shades violet, green, 
yellowish 0-024—0-035 
Elbaite colourless, pink, colourless 
pale blue 0:015-0-023 
Dravite pale yellow colourless 0:022-0:025 
Chromium green to bluish pale yellow 
tourmaline green 0:035-0:046 


Alty (1933) and others have described the authigenic portion to be colour- 
less and the detrital grain as strongly dichroic. But the overgrowths on the 
tourmalines described here show distinct dichroism, and in some may even be 
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considered to be strong. The average birefringence of the host is 0:022 while 
that of the overgrowth is 0-013 (as determined by Berek’s compensator). 
They are optically negative. 

Among the tourmalines, apart from the authigenic development at the 
termination of the ‘c’ axis (Aliy 1933—Figs. 1, 3 and 7), growths are also 
observed on the prismatic face of the host (Fig. 6), and in some, the over- 
growths completely envelop the host grains (Figs. 2, 4, 5 and 8). 

The overgrowths greatly vary in size. A few of them have attained the 
size of the host while there are some whose development is so much stinted 
as to appear as mere rims around the detrital grain or as small protuberances. 
Some others do show an exaggerated growth which may be twice or even 
thrice the size of the nucleus (Figs. 4, 5). The largest overgrowth noticed so 
far is that of elbaite measuring 0-300 0-017 mm which has developed on a 
detrital grain of schorlite measuring 0-1500-117 mm. The sizes (in mm) 
of the tourmaline grains and the overgrowths are tabulated below. 


Maximum Minimum Average 
Host Js 0:8x 0:3 0:045 x 0:03 0-153 x 0-096 
Overgrowth .. 0:3 x 0:067 0:024 x 0:003 0-088 x 0-027 


Among the host grains schorlite is the most abundant contributing about 
66-2 per cent of the total, followed by dravite which is about 14-3 per cent. 
Elbaite and chromium tourmaline are the least, accounting for the rest with 
9-6 per cent each. 

On the other hand, among the overgrowths, elbaite accounts for 61:9 
per cent, followed by chromium tourmaline 23-9 per cent, schorlite 9-6 per 
cent and lastly dravite which is about 4:9 per cent. 


Nuclei Overgrowths 
percentage percentage 
Schorlite a 66:2 9:6 
Elbaite a 9:6 61-9 
Dravite N 14:3 4-9 
Chromium 9:6 23-9 x 


tourmaline 


A ——M—M——————————— 
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12 
9:05mm 0:05mm 


FiGs.1,39and 7. Authigenic development at the termination of ‘ce’ axis. Fra. 0. Authigenic 
growth on tho prismatic face of the host. Fias. 2, 4,5, 8. Authigenie growths completely envelop- 
ing the host grains. The zone of the roots or reorganization clearly observed in Figs. 5 and 8. 
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The host grains usually contain inclusions whereas the overgrowths are 
generally devoid of them and invariably clear. It is not possible to identify 
the inclusions due to their submicroscopic nature. The overgrowths have a 
relatively lower relief than the host and are normally found to be angular, 
with striations and slight fractures. The zone of roots or reorganization is 
clearly displayed by some of the grains (Figs. 5, 8). In most of the grains 
perfect optical continuity is observed between the overgrowth and the host, 
but in a few, slight deviation in optical continuity may be observed as in- 
dicated by slight difference in the extinction position of the two. This may 
be attributed to slight variations in their optical orientation. 


DISCUSSION 


Alty (1933) reported authigenic tourmalines from the Lower Devonian 
sediments of Canada and Krynine (1946) from various geological units. In 
India they have been reported from the Sathyavedu stage (Gundu Rao 1952), 
the Gulcheru stage of the Cuddapahs (Viswanathiah ef al. 1963), the Lower 
Kaladgis of Badami area (Viswanathiah and Govinda Rao Sindhia 1966), from 
the Lower Kaladgis of Jamkhandi (Govinda Rajulu and Nagaraja 1968), and 
from the Vindhyan formations of Sone Valley, Mirzapur district (Awasthi 1961). 

Various authors in the last few years have drawn attention to some of 
the factors connected with the authigenic tourmalines. Special mention 
should be made of Krynine's (1946) excellent work on tourmalines which is 
highly informative and has proved to be extremely useful in the present study. 
Based on the colour and morphology of the grains he has recognized five 
types of major provenances. ¡Making use of his concept, the source for the 
tourmalines under report can be traced to granitic and pegmatitic provenances. 
Thus the small idiomorphic tourmalines of dark brown, green or pink colours, 
with cavities and bubbles are assumed to have been derived from a granitie 
terrain. Similarly the bigger bluish to greenish blue grains, with fewer in- 
clusions, are suggestive of a pegmatitic parentage. 

There is à marked fluctuation in the degree of roundness of the detrital 
tourmaline grains. Various shapes ranging from angular to well rounded can 
be observed. Apart from the grains with authigenic growths, there are also 
prismatie, acicular and euhedral grains which appear to be considerably fresh. 
This contrast coupled with the different degree of roundness of the grains are 
the supporting evidences in favour of different provenances. 

The authigenic tourmalines speculated to be rare and belonging to the 
anadiagenetic phase (Fairbridge 1967) have developed new growths by selec- 
tive absorption of certain chemical constituents from the percolating solutions. 
They are the probable products of a reaction between the ions of a boron 
and alumina-rich interstitial water and the primary (detrital) particles, which 
have acted as nuclei during lithification and diagenesis of the parent rock in 
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the depositional basin. Each of the three diagenetic phases (syndiagenesis, 
anadiagenesis and epidiagenesis) are characterized by certain suite of authi- 
genic minerals. - 

Apart from tourmalines, authigenic growths are also observed in quartz, 
zircon and sericite. The development of these new growths as a result of 
diagenesis breeds an interesting study and the process itself is known to be 
influenced by the nature of the overlying basin water and their fluctuations of 
oxidation-reduction potential (Eh), hydrogen ion potential (pH) and concen- 
tration of the concerned ions. The interstitial waters acting as agents of 
these ions, on saturation, might have led to the precipitation of diagenetic 
minerals. Itis more likely that the minerals and ionic solutions which existed 
in the fresh Kaladgi sediments were in a state of chemical inequilibrium. To 
fabricate a state of equilibrium, a chain of events could have been set in motion, 
whose rates and directions were controlled by the environment. During this 
resuscitation of equilibrium, selective absorption by ageing gels have led to new 
minerals and can, therefore, be considered as a vital process in authigenesis. 


[^] 
o 


Aufhigeric 
Tourmaline Percent 
D 
o 


o 


A B c D 
— — Classified units in ascending order —»=— 


Fre. 9. Graphical representation of the average percentage frequency of authigenie tour- 
maline for the arenaceous horizons (A, basal conglomerate; B, grey quartzite; C, pink quartzite; 
and D, pobbly quartzite). 


Furthermore, evidences like current bedding, ripple marks, sun cracks, 
rain prints, ete., found in the parent rock indicate that the authigenic growths 
have developed under shallow water conditions. 

In order to study the horizontal and vertical variations of tourmalines 
in all the stratigraphic units, 56 representative rock specimens were selected 
from 14 traverses at regular intervals. Thus each traverse had four speci- 
mens representing all the four horizons. Heavy minerals were separated from 
all these specimens and the tourmaline grains were counted. The average 
number of tourmaline grains for each traverse works out to be 566 grains, 
out of which 46 grains show authigenic development and the rest are non- 
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authigenic. Lateral variations of tourmaline grains are completely neglected 
as they hardly show any noticeable fluctuations, and each individual unit 
maintains a uniformity from one end to the other. On the other hand, the 
vertical variations are well defined when all the four units are taken into 
consideration, as represented in the graph (Figs. 9 and 10). The average 
percentage frequency of 46 authigenic and 520 non-authigenic tourmaline 
grains for each arenaceous horizon are calculated and plotted on the graph 
accordingly. 

In Fig. 9, the representative percentages of authigenic tourmalines in 
each horizon are calculated for the average 46 grains and represented on the 
graph. The least number of authigenic growths is in the basal conglomerates 
amounting to about 4-5 per cent and gradually increases to the maximum in 
the topmost pebbly quartzitic sandstone beds to about 45-4 per cent. This 
may also suggest a gradual increase of the suitable chemical ingredients in 
the percolating waters during deposition, cementation and diagenesis of the 
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Fic. 10. Graphical representation of the average percentage frequency of non-authigenio 
tourmaline grains for the arenaceous horizons (A, basal conglomerate; B, grey quartzite; C, 
pink quartzite; and D, pobbly quartzite). 


sediments. This evidence may help to unravel some of the aspects of the 
palaeogeographical conditions. 

In Fig. 10, the percentage distribution of the average 520 grains per 
traverse of non-authigenic tourmalines is calculated and plotted in a similar 
manner. The basal conglomeratic bed contains 26:9 per cent of the total 
number of grains. Here, most of the tourmalines are prismatic, euhedral, 
pale green or pale blue in colour and sometimes colourless. There is à sudden 
drop to 5:8 per cent in the immediately overlying pink quartzitie sandstone. 
From here onwards there is a gradual inerease in the successive beds which 
reaches the maximum of 38:6 per cent in the topmost horizon of pebbly 
quartzitic sandstone. The tourmalines from the pink quartzite and upwards 
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are mostly well rounded, dark brown, green or pink in colour. But in the 
last pebbly horizon there is a mixture of well rounded, angular and also small 
prismatic euhedral grains. 

From these evidences it is surmised that the tourmalines of the basal 
conglomerate might have been derived from the underlying basement granites 
and pegmatites. The supply of tourmalines from this source seems to have 
stopped during the formation of pink quartzites. The nature of the tour- 
maline grains present here suggests a different percentage, while those that 
occur in the uppermost horizon of pebbly quartzite appear to have been derived 
from more than two terrains. 


CONCLUSION 

The clear-cut differences in the varieties of tourmalines between the host 
and the overgrowth, like elbaite on schorlite, dravite on schorlite and the 
converse, etc., may probably indicate that the detrital grain of tourmaline is 
not dependent on the identical chemical composition of the complex iso- 
morphous series of the group, to develop an authigenie growth. 

The investigation has also proved the utility of both authigenic and non- 
authigenic tourmalines in certain sedimentological studies like tracing the 
relationship between the different horizons of a basin, viz. the Kaladgis, and 
thus may be applied to similar formations elsewhere. Further the concepts 
discussed above may also prove to be helpful in solving other problems like 
classification and correlation of sediments, and reconstruction of palaeo- 
geographic conditions. Similarity of depositional environment can also be 
traced as between the four stratigraphic horizons of the Lower Kaladgis 
exposed around Saundatti and Jamkhandi respectively (Govinda Rajulu and 
Nagaraja 1968). It is significant to point out that among the authigenic 
tourmalines observed in the Kaladgis of Jamkhandi, schorlite and elbaite 
constitute 60 and. 50 per cent respectively among the host grains and the over- 
growths. Almost similar data, viz. schorlite (66-2 per cent among the hosts) 
and elbaite (61-9 per cent among the overgrowths), are obtained in the present 
investigation where the Lower Kaladgi sandstones of Saundatti are involved. 

The development of authigenic minerals in general throws light on the 
physico-chemical environment which existed in the sediments upon deposition 
and also helps in defining the physico-chemical boundary limits of diagenesis. 
We are sure that similar data will be of great use when the tourmalines of the 
entire Kaladgi basin and also its equivalent formations are studied. 
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In this paper our aim is to prove the analogues of the results of Cassels (1952a, 
1952b) in the field X£(t) of formal power series in an indeterminate ¢ over the 
finite field K and equipped with the usual non-Archimedean valuation. Denote 
by P, the set of n-tuples of polynomials in 4 over K. Then we prove the 
following theorem. 


Theorem 1.—Let K be a finite field in which 224+m= 0 is not soluble for 
some integer m and for which the characteristic is not equal to 2. Let Ly, 
Ly, Lg, L, be four linear forms in variables x, Y, z, &.with coefficients from 
K{t} and determinant A ~0. Then we can find elements Tos Yor Zos So in 
K{t} such that 


[A^] 


| Z9 E 


AX) +mLi(X) | > 


for all X = (v, y, 2, 8) = (vo, Yo, 29, 89) (mod P, 1). 

Following similar methods we can also prove the following theorems. 
Menor 2.—Let K be a finite field. Let Li, La be two linear forms in the 
variables x, y with coefficients in K(/) and determinant A 240. Then Lo € 
K(t) and yy e K{t} can be found such that 


| (X) bs (X)] > 


for all X = (v, y) = (a, Yo) (mod oe 

Theorem 3.—Let K be a finite field in which a?-|-m is irreducible for some in- 
teger m and for which the characteristic is not equal to 2. Let Li, Lo, Lg be 
three linear forms in variables v, y, z with coefficients from K{t} and deter- 
minant A 250. Then there exist to, yo, 2; in K(t) such that 


a 


| Za(X) |. | 22(2x)-+mz2(x) zu 


for all X = (v, y, 2) = (2, Yo, 20) (mod Jy 


f * This is a part of tho author's dissertation approved for the Ph.D. TR of Ohio State 
University in December 1965 written under the guidance of Professor R. P. Bambah. 
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$1. Let K be any field and ¿ be an indeterminate. Let K[t] be the 
ring of polynomials over K, K(t) the field of rational functions and K(t) the 
field of power series in £^! over K. ‘The elements of K{t} are of the form 


U = Amt” ou MT... (up to — co). 
In the field K{t}, we define the usual valuation || as follows: 
(i) [0]=0. 
(ii) If x = apio 1171... tm FO then |v] 2 en. 
We also define 
[lel] = læ- Hæ- ... |. 
Write Py, = (K = (X = (ay, ..., Gn); a € K[t]) and 
Rn = (K{})" = (X = (z1, ..., €n); ve K(t)). For 
X = (%,..., Un) € Rn, we define 
|X| = max (12,1, ..., [n]). 


It is well known that the fields K(t), K(f) and the integral domain K[f] 
are similar in their relative structure to the real field, the rational field and the 
ring of rational integers respectively. ||x|| corresponds to the fractional 
part of the real number x, P; corresponds to n-dimensional fundamental 
lattice and Ry corresponds to n-dimensional Euclidean space. 

$2. Let Ly, Lo,..., Lryos be linear forms in r+2s = n variables such 
that Lı,..., Lr have real coefficients and Ly4;,..., Lr; have complex 
coefficients and Lyis41,..., Lryos are complex conjugates to Lr, ..., Lrss 
respectively. Let A 250 be the determinant of these forms. In the cases 
r= 2, s= 0; r= l, s = l andr=0, s = 2, Davenport (1950a, 19505, 1951) 
has proved the existence of positive constants C, s and reals £j ..., £4 
such that 


|DiLo. . . Ln] > 5 


for all (%1,..., $5) = (én ..., £n) (mod 1). He proved the theorems with 
Co, 9 = 128, Cy, = 8X 1018 and Co, o = 10132, 

Cassels (1952a, 1952b) gave an alternative proof with Cy, y = 45-2, €, , = 
420 and Cy,» = 5300. The best values of these C's, however, are not known. 

The above-mentioned results have their obvious analogues in connection 
with the field of formal power series. 

Following the techniques of Davenport, Armitage (1957a@) proved the 
analogues for finite K in the case of quadratic and cubic forms, but limiting 
himself only to the cases, where the product Lı Lg... Ln does not represent 
zero non-trivially. 

Following the methods of Cassels (1952a, 19525) in the real case, the author 
proved the following theorems in his dissertation. 
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Theorem 1—Let K be a finite field. Let L,, La be two linear forms in 
variables v, y with coefficients in K(f) and determinant A 7&0. Then there 


exists (Zo, Yo) € Re such that 


IZ;(3. 3] LS) 


for all X = (v, y) = (vo, Yo) (mod Pe). 
Theorem 2.—Let K be a finite field in which a?-4-m is irreducible over K 


for some integer m and characteristic of K is not equal to 2. Let Li, Lo, Lg be 
three linear forms in variables v, y, z with coefficients from K(t) and deter- 
minant A < 0. Then there exists (vo, Yo, Zo) € Rz such that 


[LO (2300 4-29) | > LA 
for all X = (v, y, 2) = (Xo, yo, 20) (mod Pg)*. 

Theorem 3.—Let K be a finite field in which a?-|-m = 0 is not soluble for 
some integer m and characteristic of K is not equal to 2. Let L4, Lo, Lg, L4 be 
four linear forms in variables v, y, z, s with coefficients from K{t} and deter: 
minant A 750. Then there exists (a, Yo, 29, So) € Ry, such that 


Lal 


| EXC mE (X)|.| DX) +mL(X) MS 


for all X = (a, y, 2, s) = (Lo, Yo; Zo» So) (mod P,)*. 

REMARK: The results of all these theorems are best possible. 

The three theorems mentioned above may be proved by similar tech- 
niques. Since Armitage (19574) proved parts of Theorems 1 and 2, we shall 
give the details of Theorem 3 only. However, whether we can choose 2, Yo, 
29, So in K(t), when the product has coefficients in K(t), remains open. It 
may be remarked that Armitage solved the problem for quadratic and cubic 
cases, when the product does not represent zero non-trivially. 

$3. Before proving the actual theorem, we need a few preliminaries. 

Since we are given 22-+-m = 0 is not soluble in K, it is also not soluble in 
K(tj. Adjoin «, the root of «?+-m = 0 (in the algebraic closure of K{t}) to 
the field K{t} and we get the field K{t} (x) (to be denoted by .K(" {t} hereafter). 
The elements of Ktm {t} are of the form z = a--ay where v, y e K(t). We 
define a function 

f: Kim (t > the set of non-negative reals by 


$) =J(+ay) = max (a, jy!) 
It can be checked that fis a valuation in K™ (f. Since the restriction of 
f to the field K(t) is the original valuation, without ambiguity we shall write 


Fe) = i21. 


* In my dissertation tho theorems are stated for m — 1. I am grateful to Professor Bate- 
man for pointing out that tho same proof works for general m. 
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Definition.—Let z = x--oy e K™ {t}. Then we define 
|| || = max (| 211, lg il- 

Definition. —We shall call the elements 2, = v;--oyi and 22 = v,—oyi 
‘conjugates’ of each other and we shall write 2, = £s or zs = 21 

REMARK: If 2, = Zo, then |21] = |22|. 
Suppose the given linear forms are 

LX) = ajz+bjyytcz+djs (j =1, 2, 3, 4). 
Write 
L = L¡+aLa, L, = Lb, 


L, = L3+aL,, L, = L4—a«LD,, (02 —m). 


Let A’ be the determinant of Ly, Lo, Ly, Lj. It can be checked very 
easily that A’ = —4mA. Thus A' e K(t). 

Since 224+m = 0 is not soluble in the field K and the characteristic of 
the field is not 2, 

[A =|—4mA|=|Al. 

Let M, Ma, Mg, M, be forms in variables u, v, w, r such that the coeff- 
cients of M; are cofactors of the corresponding coefficients of L; (j= 1, 2, 3, 4) 
respectively. Then 

L,M,+L,My+L¿My+L,M,= A'(xu+yo+2w+sr). .. (3.1) 

It can be checked that Mı = Ms and My = M, and det (My, Mo, Ma, 
M,) = A”. 

Now we prove a few lemmas. 


Lemma 1.—Let Ao, M, ..., An, ... be a finite or infinite sequence from the 
field KK {t} such that 
[Ansa | = elAn |; `£ 0. 


Let po ui ... be arbitrary given elements of K(). 
Then there exists é e KW (f) such that 


| | An£--ÀnE in | | Ego (m = 0, 1, DE 
Proor: First we note the following. Suppose 
Anf = o-bar. 
Then È 
Ane = 0—0T 
and zA 
Ané+Ané = 20 € K(t). 
So that for A&£-l-3a£, we have to fix our attention on ‘o° of Ané only. 
Without loss of generality, we can assume that the degree of un is at most 
—1. Suppose : 
An = A ena M ASS sh 
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, , (n) (n) +. 
where a’s and b’s belong to K and at least one of aj, and bj, is not zero. 


Let M 
- )),-2 
pn m rper Res. 


Take l = er. M 
E= aq al Ipaa, al o p. By at o +B_4,-2¢ HM 
where «’s end f’s belong to K. We shall prove that with suitable choice of 
ws and B's, we have the required result. We choose these by induction. 
Coefficient of 71 in £4-Xo£--- o = 2 (ajo. y. 3 — mb) Bog, 1) 10. 

Since at least one of e b is not zero, we can choose & 5, 1, B_j-1€ K 
such that coefficient of 2 in Ag£-FAs£--uo is not zero. Suppose we have 
selected already 

o o y Po o Pa 
Since je >jx-11, choose o-(j, 42). -+o -ip B-(jp-142) ris B- jj, arbitrarily 
and 0-(4+1) B-(j,+1) such that t7! has a non-zero coefficient in 
ME PARE x 
k k 19 k) MU 
A A m8- eene - HB tia) er 


This is possible, since (a0, UP) 4 (0, 0). 

Thus the induction is complete and lemma is proved. 
In the subsequent discussion X stands for a finite field. 

Lemma 2.—Let L,(U),..., En(U) be n linear forms over K{t} in n vari- 
ables 4;,..., Un with determinant A 750. Then for every constant C> 0, 
there exists only a finite number of elements U of P; satisfying 

max (|L,(U)|,.--,|Zn(U)|) < €. 
Proor: We define a new function 
F(U) = max (|Z,(U)],--.,|LZn(V)|) 
and by a result of Mahler (1941), we get 
F(U) z v|U| 

for a constant y > 0 and for all U e Ry and our lemma follows immediately 
from this. 
Since K is a finite field, K{¢} is locally compact and we have 

Lemma 3.—Every infinite bounded sequence in  (/) (in the sense of valua- 
tion) has a convergent subsequence. 

Lemma 4.—Let p, q be integers such that 

enta > e73| A |8. 
Then there exists a non-zero vector U = (u, v, w, r) e P, such that 
[M,(0)|=|M,(0)| < e», | M((U)| =| M(U)] < et (3.2) 


FORMAL POWER SERIES 689 


Further the inequalities in (3.2) have only a finite number of solutions from 
P,. 
Proor: Since M, = jf; and Mg = jf, let 
M, = Np +N, Ma = N,—«Ns 


Ms = N3+oN4, Ms = Ng—oN4 


where Vy, Ns, Ns, N; are linear forms with coefficients from K{t}. Then we 
can see 


0% [det (M, Me, Ms, My)|=|det (Ny, Ne, Ng, N4)| | AP. 


Since e?7t®¢ > e73| A|3 by a theorem of Mahler (see Armitage 1957b), there 
exists a non-zero U = (u, v, w, r) e P, satisfying 


[W1(0)]| <e?, |Na(U)| xe», | N«(U)| <e2, |.N4(U)|<e2 e (3.3) 
But, by definition, for any U e P,, 


| M,(U)] = | Mf(T)| = max (| N:(0)), [N2(0)) 
and 

| Mg(U)| = | M4(U)| = max (| N(U)], | N«(U)]- 
So non-zero U e P,, which satisfies (3.3), must also satisfy (3.2) and this proves 
the first part. 
To prove that (3.2) has only a finite number of solutions U € P,, we first 
observe that any solution of (3.2) is also a solution of (3.3). But, by Lemma 
2, (3.3) has only a finite number of solutions U = (u, v, w, r)eP,. This 
proves the lemma. 

Remark: Let K be a field (not necessarily finite), in which a?--m = 0 

is not soluble. Let Li, Lə, Ly, L4 be four linear forms in variables (xv, y, 2, 8) 
= X with determinant A 240. Then Lemma 4 incidentally proves that 
there exists a non-zero vector X e Py, satisfying 


| L2(X)-+mL5(X)| |L X) 4m LX) | <e] A |. 


Now we distinguish the following three cases: 
(i) M, 40% Mo, My 407% M, for any non-zero Ue Py. 
(ii) One and only one of the pairs (M, My), (Ws, Wa), say (Ms, Ma) 
assumes the value (0, 0) for some non-zero vector U e P,. 
(iii) M, = 0 = Ma for some non-zero Ue P, and Mg = 0 = M, for 
another non-zero vector U e P,. 


$4. Proof of Theorem 3 in Case I 


First we shall prove a few more lemmas. 
Lemma 5.—Let M > 1 be any given integer. Then there exists a finite 
sequence of values on, Bn = an, Yn, 9n = ya(0 € n < N) of My, Ma, Ms, M, 
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respectively corresponding to a sequence Qj; of elements (Un, Un, Wn, Tn) € P, 
with ged (Un, Un; Wn, rs) = 1 and satisfying the following properties: 

(i) |; |< e7?1 Aj. 

(ii) [an| < e Jan-11, (1 > 1). 

(iii) b? Y lA (n > 1). 

(iv) [Yo] € e%, Jay] > e". 

(v) lay] < eV. 

(vi) [Yn] > |Yn-1]; (v > 1). 

Proor: Let | A|=e%. Consider the solutions of 


3d 
|M (u, v, w, r)| = | Wolu, v, w, r)| < max (ari, a 5 ]-2) 2 (41) 
|Mg(u, v, w, r) | =| Milu, v, w, r)| < e- M 

By Lemma 4, the above inequalities have a non-zero solution U = (wu, v, w, 

r) €.P4 and the number of such solutions is finite. Consider 
min | JM s(u, v, w, 7)| 
as (u, v, w, r) runs over non-zero solutions from P, of inequalities (4.1). Since 
the number of these solutions is finite, this minimum is attained and by hy- 
pothesis, it is not zero. Let this minimum be e. Then T > M—1. 
Consider the solutions of 
| 3 a 

(Milu, v, w, 1)]=|Mo(u, v, w, r)| ze +2 i n HA 

|Mslu, v, w, 7) | = | Ma(u, v, w, 1)]  e- 771 
By Lemma 4, the inequalities (4.2) have at least one non-zero solution U = 
(u, v, w, r)eP,. Further the number of such solutions is finite. Hence min 
| Mg(u, v, w, r)| taken over non-zero vectors in P, satisfying (4.2) is attained, 
say at (Uo, Vo, Wo, ro). Obviously, the ged (uo, vo, Wo, rg) = 1. 
Let % = Mi(Uo, vo, Wo, Yo), Bo = Malto, vo, Wo, 1o), 
Yo = Malo, Vo, Wo, 70) and 8g = M (tto, Vo, Wo, ro). Then 
1%] =|8| x e-7-1 x e (since T > M—1). 

Also for every non-zero U = (u, v, w, r) e P, satisfying (4.2), we must 

have |M,(U)|>e%71. Otherwise, we have | M,(U)| <e*-1 and then this 
U satisfies (4.1), which contradicts the minimal nature of T. Hence, in 
particular | «¿| > e”. 
Since (uo, Vo, Wo, ro) satisfies (4.2), (us, v, Wo, ro) also satisfies condition (i) of 
Lemma 5. Suppose there exists non-zero U = (u, v, w, r) e P, satisfying 
|41,(0)|<e ||. Then we claim |M3(U)| >|Y0]. For suppose we have 
|-Ms(U)| «|vo|. Then this non-zero U also satisfies (4.2), and this con- 
tradicts the minimal nature of |», |. 
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Now we shall construct by induction a sequence Jar satisfying the condi- 
tions of Lemma 5. (wo, vs, Wo, rg) satisfies the conditions of the lemma. 
Suppose (to, Vo, Wo, To), . ++ , (Un=1) Unt, Way, *».1) have been constructed. 
Let 

a = M (ui, 04, wi, 14), Bi = Malut, 4, wt, 14) 


yi = Mg (us, vi, we, ri) and 8; = M (Ui, vi, wy, ri) for all 4. 


Consider the solutions of 


| M(u, v, w, r)] =|Mo(u, v, w, 7)| < elena], 
| (4.3) 


3d 
|Ma(u, v, w, r)| = |.Ma(u, v, w, r)| < le al! 


By Lemma 4, there exists at least one non-zero solution U = (u, v, w, r) 
€ P, satisfying (4.3) and there are only a finite number of such solutions. As 
before, min | My(U) |, where minimum is being taken over all non-zero U eP, 
satisfying (4.3) is attained, say at (wn, Un, Wn, ra). Conditions (i), (ii), (iii) of 
the lemma are obviously satisfied. It must satisfy also condition (vi), for 
otherwise, it would have been selected at the previous stage. Since |an|< 
l| an1] for all n, we get integer N’ such that 


| &n | c eM 


for all n > N’. We stop at such an AN". 


REMARK: 1. Clearly Sy; is not unique. 
9. If M, > Me, every Oy, is an Jy also. 


Lemma 6.—There exists 


(a) (an) an) 
(% > Y > % >» % Ro 


such that 
[sunt Un twn sin | = eat 


for all (a, Y, 2, 8) = Gs Js Vsus am shu) (mod Pj) 


and for all n(0 < n < N). 
Proor: Consider 


Was By =0 (n=0,1,..., N). 


Owing to condition (ii) of Lemma 5, 
|An| = elaa-il (GS Wy choy AN 
We can, therefore, apply Lemma 1 and get ¿ e KW (2), such that 


anët ČnÈ 


meth | ca Iss e 
A' A' 
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Let ¿=x9+0L; n, Le K{t}. Define 
eS, yee el SD) € R, by the equations 


1 MI g gar 
Le, y yen, 200, UD) — s, L(a, ye EP ) (An — 


L ales AM) Us 205 a) Ow (a M) y Ed es Ge) 


Then Gey e A S) also satisfies 


(M M 200, (M) (3) yen, AM) AM E 
L(x Jl jn son) 6 Lia, , ay o m JE 


o 


D 


r^ p. (ML M) M) s00 va (M) yn JM) — Qu) 
L(x , uh > 2 )=0= L, o s Yo» % > 90 ) 
Substituting x = n y= yen z= aan s= Dd 5 


U = Un, V = Vn, W = Wn and r = rp in (3.1), we get 


uy Fyon 2p, Ot es ROMS 
Hence 
E atu nV ontzi wn ts pn LA |l 
Santon =el0<n<N). 
A 
Further if (2, y, z, s) = EE, yp, AM, ap (mod P), 
we have 


(ttg - Jn tHzwn+srn)— — (af un yo, +20 +5; n) 
= (e—a ent =p ent (= 25) yn (5759 n 
and the right hand side is an element of K[t]. Consequently 
|] en +40, +20p + 875 || = Il aU Du, y vna wte rn | =elO<n<N). 


Hence 
| $n -d-y?n-d- zwn -8ra | > || vn+yen+2wn+s7n || = e410 < n <N). 
This proves the lemma. 
Lemma 7.—There exists a doubly infinite sequence of polynomials (Un, 
Un; Wn, rn) (00 < m < co) with ged (Un, Vn, Wn, r4) = 1, such that conditions 
(1), (ii), (iii) and (vi) of Lemma 5 hold and 


lim | Mi(us, Up, Wns ;)|- 0 = lim |Ms(un, Un, Wn, tn) |. 
n—> o n=- —00 
Also 
lim |14g(Un, Vn, Wn, Tn) |= œ = lim |M (un, Vn, Wn, Tn) |- 
no n- -o 


Also there exists (%o, Yo, Zo, So) €.R4 such that for all (x, y, 2, s) = (Yo, Yo» 
Zo 89) (mod P4), the condition 


s [tun +Y0n 42007 +srn | > e (— 00 <n < oo) 
is satisfied, 
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Proor: Let L be any given integer. First we assert that there are 
only a finite number of pairs of non-zero elements (u, v, w, r), (w, 9, w, r’) € P, 
satisfying 

(a) | Mw, v, w, 7) P| Mg(u”, v, w^, i) E <| Aj, 

(b) [Mi(u, v, w, 7r) 2] Mg(u, v, w, 1) 2. < e? | A B. 

(c) Mi, v, w, r')| <el < et] My(u, v, w, 1) | 
From (a) and (c) we get 

[Ma(e”, v’, w', 1") 2 < en?2L2] A [. 
Therefore 

| Malu’, v, w, r')| < e-L-1+134/2, 
But for (w,v,w,r)eP,, 

| My(u', v, w, ')| = | Me (a’, v, w, 7°) | 
and 

| Dfs(u’, v, w, 2’) | =| Uw, o, w, )]. 
So we have 

| M (u^, v, w, 2’) | =|Mo(w, v, w, 2’) | < el [from (c)] 
and 

(Ulu, v', w', 7’) | = [.Ma(u', v', w', 1) | < e- F-14821, 
By Lemma 4, there exist only a finite number of (u', v', w', r') e P, satisfying 


the above. 
Further given non-zero (w', v', w’, 1) e P,, we assert there are only a 


finite number of non-zero (u, v, w, r) eP,. Since by assumption, M (w’, v', 
w', 1') 40, from (a) we get 
(Miu, v, w, 1) 2 | AP. |M, o, w, 2”) 152, 
From (b) we get, 
| Mu, v, w, 1) £ < e? AP] Mi(u, v, w, r) 2 
< e-A-?L| ^ B. 
Now by using Lemma 2, we can see that there are only a finite number of 


(u, v, w, r) € P, satisfying these. 
Now suppose M > 1 is a given integer. Construct a sequence Jar satis- 


fying the conditions of Lemma 5. By conditions (iv), (v) and (ii) of Lemma 5, 
there exists a unique integer N, such that N z N, > 1 and 
jan] < 1 & e! | axa]: 


D D 


; : 
Let 9;, be the sequence of (tin v. t^, 7,) defined by 


PR a X 
(e Vas Ons n) = (tana Unt CntNy Pan) for 


—Ni sus Ns = N—N. 
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Let æ, B Y, ò, be the values assumed by Mi, My, My, M, respectively at 


(uj, vi, wh. 7,). Then conditions (i), (ii), (iii) and (vi) of Lemma 6 hold for 
(e De Ww. r) and 
ES | ES etja] , les. | SO EA ege e 


We now construct Y by a diagonal process. By the observation already 
made, there are only a finite number of possibilities for 


D / , 


(do pr) oo woa) 
and so at least one pair of these must occur infinitely often, as M runs 
through the set of natural numbers. We fix one such possible pair, say 
(do, Do, Do, Fo), (G1, 9.1, -r 7.1). We now consider only those Sip ’s which 
contain these. 

For these M, conditions (4.4) imply that for all M large enough Ne = 1, 
for otherwise My (to, o, Do, Fo) = 0, which contradicts the hypothesis. 

Denote by än, By, Yn, 5n the values assumed by My, Mo, Mg, M, respectively 
at (Un, Un, Wn, fn). Now 

[M (ui vi o, 7| < e ag] < et | M (tto, Do, Do, Fo)| 
for each Syy. 

Again, there are only a finite number of possible choices for (uj, N Ww, i) 
eP, We pick out (ŭi, 01, 4, 7,) e P4 which occurs infinitely often. Let 9, 
be the infinite subsequence of 9, consisting of those 5,,’s to which (i4, 2, à, 
7) belong. We now consider only these 9,,s. For these 9;,'s, if M is large 
enough, N, > 2. 

Consider 


D 


]Mi(&.1, 9.1, 01, 1) € |& a] E e^! | Mw. s, v» W_os 222) 

for all those 9,,'s. As before, there are only a finite number of choices for 
(Qe 9 o) w 2 4 2) € P,. We pick out (à. », 0_2, t9 9, F_2), which occurs infinitely 
often. Let 95 be the infinite subsequence of 9, composed of those 9;,s to 
which (. », 0_2, 0.2, ?.«) belongs and we consider only these 9,,'s. For these 
Iis, No > 2 if M is large enough. 

Continuing this process, we get a doubly infinite sequence (ün, Un, n, Fn) 
(— c0 « m « oo), which satisfies conditions (i), (ii), (iii) and (vi) of Lemma 5. 

Obviously, by construction 


lim |g,|— 0, lim fän] = oo. 
n> 0 n-» —00 


Due to condition (i) of this lemma, we must have 
lim bz3 = n 


no -o 


6n 
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Also lim |¥n|= co, for otherwise, there exists an integer k such that 
n— co 


[7n] < ek 
for all n. Then 
: LM (u, v, w, 1)|=|Mo(u, v, w, r)| 1 


[Ma(u, v, w, r)| =|M,(u, v, w, r)| € er 


has infinitely many solutions (u, v, w, r) e P4, which by using Lemma 2, we 
can see, is not true. "Thus this sequence has all the required properties. 
Since there is no danger of confusion, we now suppress the bar on (Un, Un, 
Wn, Tn) and Gn, En, Yn, Òn. 
During our construction òf double sequence, we have used the subse- 
quences 9,, Dj, ... of the sequences 9;,'s which have the properties 
(i) Dip is a subsequence of 9; 
(ii) Given any finite set of (uy, vj, wj, r;)'s, this is a subset of ©; for j 
large. 
Let 9; be the first member of 9;. Then the sequences 9; has the 
following properties: 
(1) 8; = Tm, for some mj; m; — oo as j > co. 
(2) Each finite set of (uy, vj, wy, 14) is a subset of 9; for all j large enough. 
By Lemma 6, there exists a sequence 


(009, yn, alm), 0 j= EID aR 


and we may assume that 
(my) mj) (my) (17) | 
lci Je «n fe [ <1, [e] <1 
for all j. 
By using Lemma 3, we can see that 


(Gem s yd, apu gut). j= = l, ma 


D 


has a convergent subsequence tending to a limit point (xo, Yo, Žo: So) (say). 

Without loss of generality, assume this itself is convergent. This choice of 

(Lo, Yo» Zo» So) does what is required. For take (x, y, 2, 8) = (tos Yo, Zo» So) 

(mod P.). Let (Un, Un, Wn, Tu) be any member of this double sequence. For 
, * > 

all j large enough, (tn, Un, Wn, n) € 9; = Jy (say), 


and 
[zun t 0n 20 ESPA | 
= | («wo +20) un + (y= Yot Js Non (z— —20 +2 De nt(s—sots rat 


(«o eunt (uoy wn + (292 o (s SoS) Pral. 
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k 
(eat, y— Vo, aah) sse) 
(ES, yP, z MR ET mod P), 


| Gr rote Jus 4-(y— BN ose cei | > e! for all n. 
If j— œ, k=> oo; therefore, if j is large enough 


| Min +(Yo— yP jont (20-2 wat (80-8) rn | <el 
and hence 
(eunt yvn+zwn+ srn] > 671 for all n. 
This completes the proof of Lemma 7. 

Proof of Theorem 3 in Case T.—Let (vo, Yo, 2o, So) and the doubly infinite 
sequence (Un, Uy, Wn, Tn) (0 <n < co) be given by Lemma 7. Take any 
(v, y, z, 8) = (Lo, Yo, Zo So) (mod P,). Suppose a, b = à, c, d = č are the 
values assumed by L,, L,, L;, L, respectively at (x, y, z, s). 

First we assert neither of a, b, c, d is zero. For suppose a = 0 = b. 
Then from (3.1) we get for any n, 


Qa d Up d ZR H- ST y = ente 
So [vunt yt 3-20, +s | = dl S E 
A A 
(Since | cy, | =|d8»| and | A^| =| A].) Making n= —oo, we get 


[ta EY 0 +20) +87 n |>0 
and this contradicts the fact that |xun+yv.+202+sr»] >e! for all m. 
‘Hence a £0X0. Similarly c Æ 0 Æ d. 
Since |Y»| > co as n — œ and |Yn|—>0 as n—>— co and |Yn311 > | Yu; 


we can choose an integer n 
such that 


estar HER PR PO) 
From condition (iii) of Lemma 7, we get 


les LA Plat? s [1a p r M .. (4.6) 


p=lab||¢nBn|=|aP lan), 9 =|cd][Yn8,]=1c]2]7. P 


E? = |abed|| A PB = [ac] A P. 
d (4.6), we get 


k, 
k. 
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From (3.1), we get 

| A! (25 t Yvon Zn - s'n) | = | dan 2-08 -- c - d8;, | 

< max (| ax], |cYn]) 
i.e. ] £x (rts + Yn zw 87s) |? < max (| den 12, | cya B) < k. 
But from Lemma 7, we get 
[tts Yn +20n +8rm)| > em} 
| A (tts + yvn+zwn+srn) |? > e7?| A Fg. 


Hence 

k>e2|aP 
or |abed|*| A |] > e-2| af 
or |abed| >| A ]/et. 


This proves the theorem in this case. 


$5. Proof of Theorem 3 in Case II 


Here we prove the following lemmas. 

Lemma 5'.—Let M > 0 be a given integer. Then there exists a finite 
sequence of values &n, Bn = Gn, Yn, 9n = Ya (0 < n < N) of Mı, Ms, Ms, M, 
respectively corresponding to a sequence Sy, of polynomials (wn, Un, Wn, Tn) 
(0<n<WN) with ged (us, Un, Wn, Tn) = 1 satisfying conditions (i), (ii), (iii), 
(v) and (vi) of Lemma 5 and 
(iv^) : Yo =0,% 250 for O0<n<N. 

Proor: By hypothesis there exists («, B, Y, 8) € P4, (a, B, Y, 9) #0 and 
with ged (a, B, y, 8) = 1, such that 

Mya, B, Y, 8) =0= AM (a, B, Y, ò). 
Consider. the solutions of 


[M (u, v, w, r) | =|Mo(u, v, w, r| <| Wy(a, B. Y, è) | | 


5.1 
| Mg(u, v, w, 1) | = | Mau, v, w, r)|=0< t1. (5.1) 


By hypothesis, (5.1) has at least one non-zero solution («, B, Y. 8) € Py. 
Using Lemma 2, we can see that (5.1) has only a finite number of solutions 
U = (u, v, w, r) e Ps. Consider 

min | M4 (u, v, w, r)i 


the minimum being taken over all non-zero (u, v; w, r) e P4 with ged (r, v, w, r) 
— 1 satisfying (5.1). As there are only a finite number of solutions, mini- 
mum is attained, say at (uo, Vo, Wo, ro). Then JM (uto, vos tw, ro) Æ 0- 
Now suppose there exists non-zero (w, v, w, 7) € P, with ged (w, v, w, r) 
= 1, such that 
|M,(u, v, w, r)| < el [Mi(Ug. Vo: Wo, ro). 
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Then Ma(u, v, w, r) 0, for otherwise, this (u, v, w, r) contradicts the choice of 
(Uo, Vo, Wo, ro). In particular, if there exists (u1, Uy, 1, ri) satisfying condi- 
tion (ii) of Lemma 5, then y, 40. 

With this choice of (ug, vo, Wo» rp), we proceed as in Lemma 5 and de- 
termine the sequence Jy satisfying conditions (i), (ii), (iii), (v) and (vi) of 
Lemma 5. 

Since | Y» | > |Yn_1| and y, 40, we get condition (iv’) of Lemma 5. 

Lemma 6'—For the sequence obtained in Lemma 5’, there exists 


(M) (BB) (af) 
Gra > % > 9) E Ra 


such that 
| ttg +Y07 T 2205 Sn | EXON 
for all 


M M M M 
(a, y, z 8) = ($^, Ys 2 8) (mod Ps) 


and for all n (0 < n < N). 
The proof of Lemma 6' is the same as that of Lemma 6. 
Lemma 7’.—There exists an infinite sequence Y of 
(Uns 9n, Wns rn) € P4 (0 < n < co) with ged (Un, Un, Wn, 4) = 1 
satisfying conditions (i), (ii), (iii) and (vi) of Lemma 5, and 
lim |Y4|2 œ, lim |&,| 2 0 
n— 0 n -> 00 
and Yo = 0, Yn 25 0 for 0 — n< oo. 
Also there exists (vo, Yo, Zo, So) € Ra such that for all (v, y, z, s) = (y, 
Yo Zo So) (mod P,) the condition 
[aun tyvntzwn tsin] >e! (0<n< 0) 
is satisfied. 
Proor: Let (u, v, w, r) € P, be a fixed element of P,. Let L be any 
integer such that 
el < 1| M(u, v, w, r)]. 
Then we claim there are only a finite number of non-zero (w’, v', w', r') e P, 
with ged (wu, v', w', r’) = 1 satisfying 
(a) |M (u, v, w, 1) 12] Ma(w', v', w, r') E | A B, 
(b) [M (w^, v, w, r')| « el. 


The proof of this is exactly the same as in Lemma 7. 

Now choose (uo, vo, Wo, ry) as in Lemma 5'. We note here that (to, Vo, 
Wo, To) as chosen in Lemma 5' is independent of M and we use same (us, Vo, 
Wo, To) for all Jy. 

For each Jy, the corresponding (uw, vi, wy, rı) satisfies 


[Mg(uz, V1, 201, 71) |?] M 1 (wo, vo, £o; ro) | <| A IE 
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and 
LM (us, 01, Wy, 71) | < e! [ag]. 


By the above observation, there are only a finite number of (Uy, v, Wy, 
11) € P, satisfying above. So we pick out (a, Ön 01, 71) e P, which occurs 
infinitely often. Let O, be the infinite subsequence of those S,,'s to which 
(ty, 91, D,, F,) belong. Now we consider only these 89's. 
Then if M is large and Sy, € Dı, then Jy contains more than three terms. 
Now for each Jy € 9, 
| M (Uo, v2, Wa, r2) | e 1| M (à, 9, $1, 71) | 
and 
| Mg(uo, Va, Wa, ro) [| M, (ŭn 91, 01, F1) P. S LA P. 
As before, there are only a finite number of choices for (wz, Va, we, ə). 
So pick (ds, 9», Də, a) € Py which occurs in infinitely many Jy. Let Dz be the 
infinite subsequence of 9, of those Sar's to which (às, ?», 03, T) belong. 
Continuing this process, we have an infinite sequence (ün, Un, Un, Tr) (0 < m 
< co) which satisfies conditions (i), (ii), (iii), (vi) of Lemma 5 and further 
Melo, To, Wo, Fo) = 0, Malin, On, Wn, Fn) FO for n 4 0. 
Also by construction 
lim |M,(Un, Ön, Wn, Fn) | = 0. 
n— 0 


As in Lemma 7, we can see that 


lim | Mglün, Ön, Wn: Fn) | = ©. 
n-—» 00 
[n the construction of (ün, Un, On, Fn)'s, we have used subsequences 9, 
Da, ... which have the same properties as in Lemma 7. 
Following the same method as in Lemma 7, we can get (to, Yo» Zo; So) 
with the required properties. In the rest of the proof, we suppress the bar. 
Proof of Theorem 3 in Case IT.—Let (Un, Un, Wn, rn) (OS n < co) and (2%, Yo, 
Zo, $9) be given by Lemma 7’. Suppose a, b = d, c, d= are the values 
i mU au D D A 
assumed by Li, Lo, Ds, Ly respectively for some (x, y. 2, S) = (Yo, Yo: Zo» So) 
(mod P,). 
First we assert neither of a, b, c, d is zero. For suppose a = 0 =b- 
Then from (3.1), we get for every n(0 <n < 00), 
Cynt dôn 
Qa EY Un tz Wnt stn = RS 
So 
> : DOCS D 
aug -]-ytg-I-2--$r = 0 (since Yo = 0 = 89), and this contradicts Lemma 7'. 
Now suppose c = 0 — d. Then from (3.1) we get that for every n, 


ap +b 
ig tytn hewn ES Un = a Bn 
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aan +bÊn | GO, | 
or [tta +Y0n 20 d- Sra | = AGNES < lA] 


Making n — co, we get 
lim [m tytn 20-875 | =0. 


n> o 


This contradicts Lemma 7’ again. 

The proof of Theorem 3 in this case is exactly the same as in earlier case. 

$6. Proof of Theorem 3 in Case III.— M= 0 = Mo for some non-zero 
U e P, and Mg = 0 = M, for some other non-zero U e P,. 
Then we must have one of the following: 
(a) For each M > 1 there exists an Oy with the properties of Lemma 5’; or 
(b) there exists a finite sequence of values on, Bn = Gn, Yn, 8n = Fn of Mj, Mo, 
Ms, M, respectively corresponding to (Un, va, Wn, ra) € P4 with ged (un, vn, 
Wn, Tn) = 1 satisfying conditions (i), (ii), (iii) and (vi) of Lemma 5 and also the 
following: 

liv) Yo = 0 = ðo Yn 250 for0<n< N. 

(v) an =0= By, on £0 for0<n<N. 

The proof of this obviously follows from that of Lemma 5’. If we have 
(a), the rest of the proof is the same as in case II. If we have (b), we prove 
the following Jemma. 

Lemma 6".—For the finite sequence obtained above in (0), there exists 
(To, Yo. Zo» So) € R such that 

[en yn + 2n - $7] >g! v oe .. (6.1) 

for all (x, y, 2, 8) = (Lo, Yo, Zo» $9) (mod P,) and all n(0 < n < N). 

Pnoor: Since ay = 0, yy 40. So we can choose 
0 = a-4-«7; o, T € K{t} such that 

210470 || = e71. 
Owing to condition (ii) of Lemma 5, we can apply Lemma 1 with 
An = 4y_n, Pn = Yx-n0--Óy 4 0(1 sans N) and we get 


£ = 3 cale Kt); y, Le Kft}, 
such that E o 


llant+3n ¿+ 7004810 ll —e(0zn»-N). 
Since «y — 0, due to (6.1) we have 


ll aye +pye+7y0-+-5y8 li = li YNO--YNÀ ll =el 
and 


I ané +H Bnë Hyn + Gnd il SelO<ns N). a Be (6:2) 
Let to, Yo, Zo, So be the solution of 
Liv, y, 2, s) =A’, Lola, y, z, s= tA 
Lg(v, y. 2,8) S o A', Lv, y, z, 8) 2 TA. 


E. 


FORMAL POWER SERIES 701 

This (vo, Yo, Zo» So) also satisfies 

L(t, y, 2,8) =A Liu y, 2, 8) EA. 

L(x, y, 2, 8) = 0A', L(x, Y, 2, 8) = 0a’. 
From (3.1) we get 

Coun +H YoUn+zZoWnat Sorn = éant ZBa4- OYn + bòn. 
Then from (6.2) we get 
|| toun- Yon- ZoWwnF sorn MEN éant ZB 4- 8Y 53-08, |] =e! 
(0 x m x N). 


Take (x, y, 2, 8) = (29, Yo, zo, So) (mod P,). 
Then 
| Tun +Y0n +20 + Sr) = | | Um + YU +20 ES? | 
= || 292 +Y0Un +20 Sora || =elO<n<WN). 

Proof of Theorem 3.—Let (to, Yo, zo, So) € R, be given by Lemma 6”. 

Take 
(v, y, 2, 8) = (Yo, Yor Zo, 8) (mod P,). 
Let a, b, = à, c, d = be the values taken by 
Lo Lo, La, L, 

respectively at (v, y, 2, s). As in case II, we can suppose a 034b. Now 
suppose c = 0 = d. From (3.1) we get for every n 


aan + bf, 
G+ Yn ZU H- ST = end : 
So 

Ux d yUN +20y+sry = 0. 
This contradicts Lemma 6”. 


Now as in (4.5), if we can choose an integer n(1 <n < N) such that 


[Ym-11? <|A]! 


«| | a 
t s [Yn 12, 
the proof is similar as in case I. If not, we have 
a 
Ist <lar|S|. 


But from (3.1) and using «y = 0 = By, we get 


i Y 22 + __leyn+ddy| GA 
| wun tyon +zwn tsr] = TA 7] < Tal: 
So 
Jc2}jyw[? _ Jac] 
eta Ego beeen sty < TAE AE [ap 


ds 
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But by Lemma 6^, we get 
[vuy +yoy Tg Ter | > em}. 


Hence 
ac| A LA] 
i > e-2 or | a2c?| E Lal or |abcd | > Cm s 


This proves the theorem in this case. 
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RADIAL OSCILLATIONS OF COMPOSITE POLYTROPES— 
PART II 


by MANMOHAN SINGH, Department of Mathematics, University of 
Roorkee 


(Communicated by C. S. Ghosh, E.N.I.) 


(Received 22 January 1968) 


Small radial oscillations of four composite polytropic models with interfaces 
éi = 0:2094, 0:3915, 0:6160 and 0:8250, having indices 1-5 in the core and 3 
in the envelope, have been considered. ‘The anharmonic pulsations of these 
models have also been considered. 


l. INTRODUCTION 


We consider in this paper the small adiabatie radial oscillations of the 
composite polytropes having indices 1-5 in the core and 3 in the envelope. 
Four such models with interfaces at é: = 0-2094, 0:3915, 0:6160 and 0:8250 
have been obtained in the Part I [Manmohan Singh (1969). The radial 
oscillations of these models have been worked out for the first four modes 
taking « — 0-6, 0-5 and 0-4. We have also considered their anharmonic 
pulsations for the case « = 0-6. 


9. SMALL RADIAL OSCILLATIONS 


Using Emden variables £ and $, the equations of pulsation for the core 
(n = 1:5) may be written as 


dy l(, 5£d$ (5 5g) 

ases BE det Nd T3 8 dé Q) 
2 Rop, : 

where T = P oc Ee xe a 8) 


and other symbols have the same meanings as given by Prasad and Gum 
(1961). l i 

For the envelope (n = 3) the equations of pulsation in terms of variables 
é and 0 are 


dn 4 AF [5 ee) E X83) 
de tE (ue at 6 + £6 de)” -> (3) 
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Here 
Ro?p; bi 
yPe $i 
w 0; " 
=e eo a se m ss (6 
where 6, and 44 are the values of 0 and ¢ at the interface. 

Eqn. (1) has singularities at the centre (£ — 0) and also at the surface. 
But as the solution of this equation does not extend to the surface, we only 
require a solution finite at £ — 0. So we develop a solution in series near 
the centre which is non-singular at the centre £ — 0. 

With the help of well-known expansion [B.A.M. Tables (1932)] for 4 we 
have 


Qi 


A?£2 A4£4 ASES - 
eas 28 ae 80 do] n .. (5) 


The solution in series for the centre near £ — 0 is assumed to be 
n = Co--Co£?-F o4£5-F eo£0 E ... ac m .. (6) 

On substitution in eqn. (1) we can find the coefficients cs, cy, cs. ... by equat- 
ing the coefficients of various powers of £. 

The eqn. (3) for the envelope has regular singularities at £ — 0 and 1. 
A solution in series started from one extremity usually becomes divergent at 
the other. For the envelope we require a solution finite at ¿= 1 and not 
extending to the centre. To obtain it we put ¿= 1—41 in eqn. (3). We 
have 


h(i! E E de 1 dé a E 
de I= TIENDEN oo ce A 
To envelope a series solution we have taken 0 as follows: 
A i Luo 2 sis wi (8) 
where 
dé 

ey eta € 
wn (i5), E ea E jo De (9) 

and the series assumed is 
y = latta? -ast3 +... S. ss .. (10) 


On substituting (8) and (10) in (7) we can find the coefficients 44, do, (3, . . . by 
equating to zero the coefficients of various power of t. 

The solutions of eqns. (1) and (7), satisfying boundary conditions at the 
two extremities as well as the interface, have been found out in the same 
way as by Prasad (1953). 

The values of «/44 for the first four modes for three values of a (= 0:6, 
0-6 and 0-4) are given for the four models in Table T. 
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TABLE T 


Values of w| At 
a AAA 


Model with 


e [va zero-th mode Ist mode 2nd mode 3rd mode 
£ = 02094 0:6 6-023504 11:226409 19-105066 29:279197 
0:5 5:306420 10:560777 18:494241 28:689014 
0:4 4-481545 9-916510 17:899210 28-108813 
é = 0:3915 0:6 11:919778 21-318248 oe :279639 
0:5 11:402158 20-829913 32-793357 
0:4 10-903206 20:348435 32-309936 
¿é = 0:6160 0-6 14-232968 39-415017 
0-5 13:807355 38:919717 
0:4 13:382857 38:428582 
£ = 0:8250 0:6 3-498776 15-719071 29-998478 46-498867 
0:5 2-931580 15-239002 29:561448 46:048814 
0-4 2-358007 14-759121 29-124757 45:599013 


3. ANHARMONIC PULSATIONS 


For small radial oscillations the square of the amplitude is neglected in 
the equation of motion. This gives a symmetrical sine curve for the radial 
velocity of pulsating stars. But observation shows that the ‘Cepheid Vari- 
ables’ have velocity curves, which are far from being represented by a simple 
sine curve. They usually show a steep rise to a maximum and then a slower 
decline to minimum. This necessitates the consideration of the second order 
terms. 

The equations of anharmonic pulsations (Prasad 19495) for the com- 
posite polytrope may be written as 


2, Ad 
end = f D Pima a > >, Parmar a .. (1) 
wi Im) ji 
where 
r= RE Se v e A M w ad a, (12) 
T0 A 5n En S. m ss os .. (13) 
On 
fioc Y NS x vs oe rs S. .. (14) 
ei 
& 3/2 al " 
In= | doti, d+ j Bit dé .. ds S .. (15) 
0 a 
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and 
es 5[* 
Digg, = 25 ¿(9- *) (3y+1) H $3I29/ ESN many, d£ 
"4g 
go? RT 
T4—— | 630' E87 Mk hm d£ | 
Pi de 
3y—1 si T Dio TORA 
ar 9 PE (nj, +7, 1871, 12h) dé 
Ti 0 


5/2 1 
Duo | eie (ng, 7, JS PER In) dé | 


1 i 


5/2) 1 
es pre ty ded Í wen ac .. (16) 


1 & 


where j, k, m are the positive integral values starting with unity, q's are un- 
known functions of time, y's are the relative amplitudes normalized to unity 
at the surface of the star and 27/o, are the periods for small oscillations, 
£1, 0; and ¢; are the values of é, 0 and ¢ at the interface. 

Ps and D's are evaluated by numerical integration. This gives the 
following differential equations for the first three modes: 


(i) E, = 0:2094 
us E 


= 0-952050. q7-L-0-593014 q,9,+0-069846 9,93 

E EN :863767 q = 4-082256 q? +7:873256 q142+2-463878 qid, SC) 
g TEH 171753 q3 = 5:144444 q?--21-600000 q192+19-511111 qiqs 

(i) é = 0:3915 


dq, 


pa FI = 0874390 gi 0-370095 g1g2+0:005027 944 


d? 
pete 145947 gz = 6-514288 q;+9:596708 q192+2:338361 ggg ) .. (18) 


d2 
A+ 3-837979 ds = 0:809235 q2-.-21-908740. q192-+21:827020 KA 
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-1 


(iii) €; =0:6160 


d?q 
dz di = 0-906197 q;-+0-264662 q,99—0:029630 gags 


12, 
"us +3-237463 q2 = 11073776 q; +17:342202 9192+3:174770 q,93 (19) 


dq a Re ; 
TE 5846325 qda = —11:054628 q, +28-307066 1192+33-9788648 9193 


(iv) & = 0-8250 


12q 2 

ata = 0-700618 9?+0-319456 q19—0-014567 qq 

1242 

a +4 492734 do = 7379650 q?--18-095526 q192-+0:180129 q193 (20) 
12, D DJ 

a + 8573002 ds = —4:044655 q?-++21-650820 q199-+38-715086 qid; 


The method of the solution of these equations is that given by Prasad 
(1949a). The solutions of eqns. (17), (18), (19) and (20) have been started by 
taking a, = 0:06. The radial velocity curves for the four cases are shown in 
Fig. 1. 


ath lgi = 0.2094 
== -~ Eis 0,3915 
m&r 0.6160 
IV fix 0.8250 


Trc. 1. The radial velocity curves for the composite polytropes. 


4. CONCLUSIONS 


From the results of $ 2 we conclude that the composite model with inter- 
face ¢; = 0-2094 is quite near to the standard model in the pulsational be- 
haviour, the composite model with interface £; = 0:3915 has drifted away 
from the standard model and the composite model with interface $; = 0-6160 
is very near to the complete polytrope » — 2. The last composite model with 
interface £i = 0:8250 has shifted closer to the complete polytrope m = l-5. 


= 
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As for the anharmonic pulsations, we see from Fig. 1 that radial velocity 
curves of model with interface at é: = 0-2094 and 0:3915 lie very near to 
each other and they have skewness near that of the standard model. In the 
last two models with interface at é: == 0-6160 and 0:8250 the radial velocity 
curves show humps and so we regard these two models unsuitable as far as 
their anharmonie pulsations are concerned. 

The comparison of the pulsation characteristics of these models with 
some other models has been given in another paper by the author (1968). 
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ON MECHANICAL RESPONSE IN A PIEZOELECTRIC ANNULAR 
DISK TRANSDUCER OWING TO A STEP INPUT WITH 
A PRESCRIBED TEMPERATURE FIELD 


by B. CHAUDHURI, Department of Mathematics, Jadavpur University, 
Calcutta 


(Communicated by B. Sen, F.N.I.) 


(Received 18 April 1968) 


The methods of transform calculus are used as tools to calculate the mechanical 
response in a piezoelectric annular disk owing to a step voltage with a pre- 
scribed temperature field. 


INTRODUCTION 


It is sufficiently well known that the piezoelectric transducers play an 
important role in the generation and detection of ultrasonic waves, vide Mason 
(1950), Cady (1946). The problems of producing mechanical response by a 
voltage step in a piezoelectric transducer, in the form of a plate and bar, have 
been investigated by Redwood (1961a, b), Sinha (1962) and Giri (1966) and 
their discussions are confined to the plate and the bar transducers. Sinha 
(1965) in his recent paper has studied the mechanical response of a piezoelectric 
annular disk transducer owing to an impulsive voltage input. ‘The present 
paper is an attempt to study the mechanical response in a piezoelectric trans- 
ducer in the form of an annular disk owing to a step voltage input with a pre- 
scribed temperature field. The solution has been facilitated by the use of 
Laplace transforms and under suitable assumptions similar to those of Red- 
wood (1961a). 


PROBLEM, FUNDAMENTAL EQUATIONS AND BOUNDARY CONDITIONS 


For the present case we need an annular disk transducer of piezoelectric 
material. Let 7, and rg be the inner and outer radii of the boundary of the 
disk. Let V be the electrical voltage in the form of a step function. Let the 
temperature at the inner boundary r— 7, be T and the outer boundary 
+ = rg be at zero temperature. Now our problem is to investigate the nature 
of the mechanieal response in the above type of transducer. 
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To serve the purpose of our problem we first need a relation connecting 
the variables representing mechanical, electrical and thermal fields. We 
choose the mechanical (radial) displacement u as the mechanical variable, the 
electrical voltage V as the electrical variable and the temperature T as the 
thermal variable. To simplify the algebra, we assume the disk to be open- 
circuited. Referred to the usual coordinates, the equations of piezoelectricity 
as suggested by Mason (1950) have been modified for the present problem and 
they are given by 


or = C DEC = hD; -AT i E 
, du y u P 

og = Cre p Un z DAT .. .. m (2) 
A: du u\ Dr : 
E = O or O0 Go oe (3) 


where (er, c9) are the components of stress, D, is the radial component of the 


electric displacement D, E, the radial component of the electric intensity Z, 
C31, C ¡2 the elastic compliances, h the piezoelectric constant, e the permittivity 
and T the temperature. 

The equation of motion in the radial direction is given by 


ND 
or ^r 00 | 


Assuming the disturbances to propagate in the radial direction, we have 
from the last equation, 


The equation of steady heat flow to be satisfied by the temperature T is 


vm =0. OPE a 2 (0) 
7B 
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The solution of eqn. (6), subject to the boundary conditions 


T=T, on r= P. 


T=0 onrzers 


T Yo 

T = — log 2 

15 Ya r 
log = 
imi 


Considering eqns. (1), (2), (4), (5) and (7), we have 


du lOu u M Ti 1 (—AX) T Oria p 3u 


Lect EE : : : ue E (8 
AA r Cn yah © Cur log E r Oy, 0? ) 
n ri 
Taking Laplace transform defined by 
00 
ü(r, p) =| u(r, tje-?!di, p>0 
0 
we have 
ai 1 dt 1 ee , 2) dt 
HEB (Lene Ja- (zen; log Gera. °° .. (9) 
where 
= (é) 
Cu 
Ko =i a O sine 
11 log = 
1 
gù TJ 
pa log 2 
ri 
Solving eqn. (9), we get 
KC (X2-+K; log 19)C31 
d = Al,(apr) + BK1(apr)+ aa log r— a (10) 


where 7, and K, are Bessel functions of imaginary argument and of order 
one, A and B are constants to be determined from boundary conditions. 


From eqn. (3) 
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Therefore the potential V across the disk is given by 


vT. 
V= H Edr = —h[(&)r = r2— (@r 2 01]— a (rg—71) log r2 
p log — 
ri 
+ x (ra log rg—73—73 lograra] zdr an ax (UT) 
p log = "i 


where we have considered the fact that the disk is open-circuited (so that Q, 
the surface charge is zero and hence D, related to Q is zero). Now eqns. (10) 
and (11) constitute the fundamental equations of the problem. 

With the value of 4 from eqns. (10) and (11) ean also be put in the follow- 
ing form: 


i KO (K54- K; log r9)C11 
e lu qj A AL Av. ae 
V t| AZ, (pro) + BK; (apre) + FEED log rz prap? 
K Cn (K+K, log r9)C 11 
—Alj(upr)—BK,(«pr1)— prip? : log? ed 


TTA(rg—7) log re 

— n E, c 5 (r2 log rg—rg—7 log 44-71) 
p log 2 p E ze 
"n n 


e la Ber, pte «| 
T r 
Ti 


igen (os r log egre) q (K ,—Ko—K, log rs)hCsi (1 (: à 
pp? Nom T2 pp? ; 
In accordance with our assumption, we have 
y = VH(t) e oe (13) 


where H(t) is Heaviside unit function equal to unity when ¢ > 0 and equal 
to zero whent <0. We, therefore, have 


(14) 


In addition to the condition (13), we assume the boundary of the disk 
to be rigidly fixed, so that 


(Qy=g, 20; - COS RENT 
Because of eqns. (14) and (15), it follows from eqns. (10) and (12) that 
L 
Al, (apre)-+BK,(apre) =- p 00 00 56 (16) 
_ M (Vo4-N) 
4a314-Ba1, = — HU emu d S sa (17) 
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where 


n tero» Eu E a .. (18) 


49 = i [rae tera »i .. (19) 
Ti 


x K,O logre (K+K, log re) Cy 
pre pT2 


Pess 2AK C31 log ry = h(Ke+K, log ro)C —AK;0 log re 
pry pry y pre 


4 Es Ra log r)C31 ( 1 1) 


p Ty Te 


TTA(rg—71) log r T. 
N= EU = D (ra log 72 —72—r1 log r,-+11). 
log T log — 
1 


Solving eqns. (16) and (17) we get 


l 1 
[MK (apra) — Lar] — — (Vo +N)Ki(apro) 
e pec 
[aiI (apra) —a31 K3(«pro)] 


1 1 
[MT (apra) —La11] pup (Vo+N)I, (pra) 


B= : 
[aK (pro) —a1911(0pr2)] 


Substituting these values of A and B in eqn. (10) and then taking the 
inverse transform of the resulting expression, we get the mechanical displace- 
ment u corresponding to V given by condition (13). From this, the displace- 
ment of the inner boundary 7 = 7; can also be calculated. This displacement, 
from eqn. (10), is given by 


E W 
(Ur=r, = AL (pri) + BE (apr) +75 


where 


iS K,O logri (Ke+K, log raJCn 
% pry p" 
By eqns. (20) and (21) 
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Its evaluation for large values of time will be possible if we expand it 
in powers of p [see Carslaw and Jaeger (1959)]. p is small for large values of 
time. Therefore, in the expansion of the relevant quantities in powers of p, 
we neglect squares and higher powers of p. 


The expansions for I,(«pr) and K,(«pr) are given by Carslaw and Jaeger 
(1959). 


T. = (kapr)2m+1 
Heg) = Dy m! I'(m4-2) 


¿mi (m1)! (m+1)! 


I EE 
« 


jam mel 
K;(opr) = (log (4apr) + YH, (apr) -1> LEMA [2 a+ > | 


Restricting to the kind of approximation referred to above, we can write 
Iy(apr) ~ 4 (apr) 


Ky (apr) ~ a; (apr Mun 


c 

where a, and b, are numerical quantities. 
Thus 

r a 

Ti y 
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? Ky (apr) it bj(1 1 

[EP ar cto 2 nta ms 

So 
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Ty (apro) | sane) dr ~ jroby (- = 1) 
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Ya ry 
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Ty 
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Hence from eqn. (25) it follows that 
Gea 
M Tag ri 


ri ~ hp3 . (a E ra) yn -( I alk 
fe r T9 pute 
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n.i 

(Vo+tN) (2 n 

a (ee ree] 
is 2j" go mE MET 


or 


> Mn Vota ( a) L(Bryra—=19 15) W 
Oran gis (o) 2hp CH Hg 2ro(11—12)p3 ty 


Taking the inverse transform, the mechanical displacement of the inner 
boundary of the disk r = rı, for large values of time, is given by 


Mte h n L(2ryrg—1;—72) We 
(rer, ~ 4h Ta 4ra(",—r») 2 


H [ates (142) Ja 
me 


2h 


Thus to conclude we say that the mechanical response owing to a step 
voltage input under thermal conditions is partly step in character, for large 
values of time. The similar result can also be had from (24). In the absence 
of temperature the result agrees with some known results. 
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MAGNETO-THERMO-ELASTIC INTERACTIONS IN AN 
INFINITE SOLID DUE TO TRANSIENT HEAT 
SOURCES 


by SACHINDRA Kumar Baxsut, Department of Mathematics, 
Ramakrishna Mission Vidyamandira, Belur Math, 
Howrah 


(Communicated by B. Sen, F.N.I.) 
(Received 29 February 1968) 


The object of the present paper is to determine the distribution of stress, 
strain, temperature, electric and magnetic fields in an infinite solid if it is 
subjected to a transient heat source. The solutions are first obtained in 
terms of Fourier and Laplace transforms, The inverse transforms are studied 
in particular cases, 


INTRODUCTION 

The studies in magneto-elasticity which has largely been possible due to 
the merger of the theories of electro-magnetism and elasticity are of recent 
origin. Apart from these, attempts are also being made in recent times to 
accommodate the thermal field in the interplay of electro-magnetic and 
mechanical fields. These studies are considered to be important primarily 
because of their diverse applications in various branches of geophysics, optics 
and in the devices such as aero-magnetic flutter and also in plasmatrons, 
plasmaguides, ete. Though there exists in abundance a literature of magneto- 
elasticity, vide Kaliski (1960a, b and 1961), Sinha (1965, 1966, 1967), Paria 
(1961), Giri (1966), Das (1967), the subject of magneto-thermo-elasticity is 
still in its formative period. The initial attempts in this direction have been 
constituted by Nowacki and Kaliski (1963), Nowacki (1962, 1963), Paria 
(1962, 1964), Wilson (1966). The papers of Sinha (1964) and Murthy (1966) 
should also be mentioned in this connection. ‘The present paper is a similar 
attempt aiming at the determination of the effects of thermal and electro- 
magnetic conditions on the distribution of temperature, stresses as well as the 
amount of the induced magnetic field. The solution is achieved by making 
use of the equations of Maxwell, equation of elasticity and equation of heat 
conduction. It is found that the use of Fourier-Laplace transforms facilitates 
the solution of the problem. à 


FoRMULATION OF THE PROBLEM, FUNDAMENTAL EQUATIONS 
Let us consider an infinite elastic solid body in which an initial magnetic 
field is impressed and having a uniform temperature. There is no initial 
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stress or strain. The body is subjected to a time-decaying heat source placed 
on & plane extending in all directions. The plane is taken as y-z plane, z- 
axis being along the original magnetic field Hz. The z-axis is perpendicular 


to this plane. The displacement u has the components (t, 0, 0), w is a function 
of x and t. 

The problem being one of magneto-thermo-elasticity, we have to solve 
Maxwells electromagnetic equations, Fourier equation of heat conduction 
together with the equation of elasticity given by Hooke-Duhamel. 

Maxwell's equations reduce to 


curl H = div B=0 


in absence of displacement current. 
These equations should be supplemented by Ohm's law given by 


a nar "T s vo AC) 
The Hooke-Duhamel law can be written as 
Tjj = 21eg4-- (Ae — BT)81 m m 5n .. (3) 
where 
eg = Hu Hu, 1), e— div u | Es .. (4) 
8,—1 fori—jand8g—0 forixj 


The Fourier law of heat conduction is 
oT Q > 
k V2T4Q = p +ToB 5 Emo div j RES) 


The equation of motion of the solid is 


CHIA Ori 


Pop = ae + (x5). aye Aes... (9 


The symbols used here have their usual meanings as in Paria (1962). 
@ represents the intensity of the heat source, To the initial uniform tem- 
perature of the solid, and T the perturbed temperature over and above To. 
From eqns. (3), (4) and (6) we get 


ou > E > > 
Pas = HV ?uRA+u) grad e—B grad TH (GxB) .. .. (7) 


SOLUTION OF THE PROBLEM 


We can assume for this problem from symmetry that all the physical 
quantities must be functions of x and £. 
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Let us take Q = qo8(z)e- (2 > 0) where $(x) is the Dirac-delta function 


and go the strength of the source. We get the following equations as given by 
Paria (1964): 


ppm Re H.20, 4220 2. a” G 
Od us — uu Hs dece Hy=0. CO MEE 
H,-0; 3:20 A 
LN A Cos qu 

Q4 30) oe — = — polly = po "EDGE c (9j 
Eg Ho = pe, T mp P TN C. 8 


js EE s : c 
where vy = zg" the magnetic viscosity and Hz = H3-+hz. (hz is so small 
» 


that its square and higher powers and its product with u can be neglected). 

Equations (11), (12) and (13) will determine hz and u, T and hence from 
eqn. (8) we get Ez and from third equation of (9), jy can be determined. 

We now proceed to obtain the solution of the problem by means of 
transforms. 

Let us define the Laplace and Fourier transforms of a function f(a, t) by 
f and f', respectively where 


He f(a, t) exp (— pt)dt, F VS zl. f. exp (¡¿x)dx. (Re p > 0) 


From equations (11), (12) and (13) we get 


(p--vu£)&, = Hs ip £u p: E .. (14) 
well sith, = (pp? -- (A--240)£2)u' — Bi£g" e +. (15) 
= {pcyp(p+2)+ (p+ 2)k£2yp' — ToBi£p(p-- Qu .. (16) 


VAS 


Solving these equations we get 


n do "n 
OS Mei i. 


o Bi£(p-Eva£?) 
1 — 2s. Lo 


MIÉ "Y ma 


722 SACHINDRA KUMAR BAKSHI 


where 


Dy = vg(ME-2pY£*--p(vupp +A +24 pe H3)£? - pp? 
Lo = {pcyp(p+2) --(p--Q)k£2)L7, +T Epp 4-QCp--vu£?) 


Putting 
R poH Dom TB 
OS e vi 

and a ^s ree MEE A) 
ky = A 


The above relations can be put in the form 


= do My 
E + A E 
peu Y 2m M, UE 
| pa- Br 
| ü AAA (19) 
(o. 9» — HspBE? 20 
k, EXE SMS : (20) 
t where 
» My = (vut (2 4 ) | 
(e f= M(vup--c;4- VN) 
3 vuo és = M(vup4-c) —/ Ñ) 
E N = vyp? —2vgp(20)—6,) -6; e »» (21) 
t 
iG A+2, 2 
== d -—d Rn) 
and 


Mo = (p+ ké?) M, +cipé2(p+vné2)er 


Using the second and third equations of (9) and the first equation of (8) and 
taking transforms, we get 


=__do Hippi 
5 POw/9g pMa 
ORE RA A 09 


mec Ko do. i£M, 
Z o "PC / gr pn 


ji taking the inverse transforms in eqns. uos a and (22) we can 


7, U, hz, Ex, Ey, jy. The other components of H, E and j j have already 
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been proved to be zero. The only non-vanishing stress component is rzz 


which can be derived from eqn. (3) in the form 
Qu 
Tzz = (A+2p) En —pT 
Taking the transforms we obtain 


e. EUER My {é (vupti Ra) +2?) 


For a perfect conductor we have o — co so that vy = 0. Therefore, 
M, = p(p+2) [ces £t +p (o 4-15 p) +p] 
M; = p(p4-Q)(cHs£ + £)(£2--£)) 


where 
= = 0% ¿hos = = o(L+R, +e, ) 
E= ip(lp+o+ VR) 
E = 5 (kipti -V R) 


2kyc, 
R = ken? —2kyp(2c2—c2) +6, 
Therefore, it follows from eqn. (18) that 


= do A B 
T = —=>"——__ |= 3th 2 
pey 2 .(p4-Q) (a EE, 
where 
o kyp—6 
A= x ]1———-— 
cand 


nva) 


Taking the inverse Fourier transform. of T', we geb 


a do 
IET 2pcypl(p--Q2) vm ez a Ale Z5) Bé exp (— -&e| 


Similarly from eqn. (20) we get 


;' qoHsB 1 [: e 5 | 


b= — a/a pora) VE ERE RÀ 


The inverse Fourier transform then gives 
qoHsB 


n a ala? ev 


(30) 


(31) 


(33) 


; -& exp (= &e)] (33) 
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From eqn. (23) we get 7,. Putting vg = 0 and then taking the inverse 
Fourier transform we have 


E doB ZA mi E cs € 
ra [zv esi) TE ow te] (30 
where : 
c 5 
A1 = oy Bevis + Rd] .. oe .. (35) 
1 ies a 


TEMPERATURE FIELD FOR SMALL THERMO-ELASTIC COUPLING 


er is called the thermo-elastic coupling factor which is very small for 
common metals. We shall expand various expressions containing er retaining 
only its first power. 
From eqn. (28) we have 
hyp-rc; 


VE= OR hres cier 2s ss so (GN) 


Hence from eqns. (26) and (27) we get 


2 
I C ET 
é =pvVk [Hse] 9 3s (SS 
3 1 kp) (38) 
i= Jz 1 eer 
ds xp T e .. (89) 


Also from eqns. (30) and (31) we obtain 


2.2 
BEC VE kıp 


re (40) 
(kxp—c;)* 
2 
1 c kiper 
B=>—|1--— 
la | ue s 
With these approximated values, we have from eqn. (32) 
= ? 2 = Pe "x 
2pcyT! A 1 oe Ap, GT | ege 6 " a —u E 
— = — a. ao + — A 
d — kp.(p4-Q) PHR | (kyp—cs)? " 2 (ipd?) 
| 1 1 ] -aE 
Sninga D E lx == n e e (42) 
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Taking the inverse Laplace transform, vide Erdelyi (1954), and introducing the 
ze; 


dimensionless quantity 7 = 7 
ei 


, We obtain 
y — 0,95 
2k Va vt mty? 
2pCyk M c= 
2poolt y -A "e| A i 
doi NONE 


where 


dr Fer (9, UN Ry) .. (43) 
0 


Ru =1+4+Ry T T oo .. (44) 
and 


cay Qu, n 
e- 9 Qh. Viv! ( e A in) 
F(n, t, Ry) = =>. 712.0 t dr 
Vs e 0 


vh E í E 
p dnd | E de 1 (n—724/ Rae Bana Ru 


X orfc (2 -1V Eu) +(n-+72\/ Ru). QR ou Ryu 


E 


MILI gr 
X erfe (Ltv 2x)? dr fort< m o 


E = 


—= Vt = AAA Am 
ky g [VA E e i EVE UVR orfo (2 zum Rs) 


72 == ACC y nki * 
ade tv Ra Ry Ry . erfc (+v 2n) } | for t > 8C Le .. (46) 


From eqn. (43), it follows that the first term on the right-hand side which 
decays with time exponentially represents the solution of the classical heat 
conduction equation while the function F(n, t, Ray) is the perturbation due to 
the thermo-elastic coupling factor er. The initial magnetic field Ry occurs 
only in the perturbed function. : 
V ROSE I vives the values of F(», t, Rar) for different values of y, Ry when 
Es o 


Qi = 0-72, the conductor being aluminium. 
8 
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TABLE I 


Values of F (7, t, Rm) in units of 10-1 


SE 0-0 0-6 1:2 1:8 2-4 3-0 3-6 
_ Rm EE AE E — 4. EO NEN . 
1-0 —3:487  —1043  —0:179 0-016 0-015 0-004 0-003 
2-0 —8:278  —2:526  —0:648  —0:089 — — 0-002 0-012 0-015 


The continuous curve as in Fig. 1 represents the perturbation when 
the magnetic field is absent while the dotted curve the perturbation function 
when the magnetic field parameter Ry, = 2. It is observed from the graph 
that the magnetic field reduces the temperature near the end (7 = 0) of the 
solid while, in the rest of the region, temperature is not much affected by the 
presence of the magnetic field. 


7 -y0 -40 -90 -ço -T9 - go 


F Opt Rus > 


Fra. 1. Effect of magnetic field parameter on perturbation 
due to thermo-elastie coupling factor. 


DISTRIBUTION OF STRESS AND INDUCED Maqnetio FIELD FOR 
SMALL Tore 


The stress and induced magnetic field are determined considering small 
values of time. In this case p is very large. Therefore, we can take 


2 
= = C- 
VR = kyp, f= VR.» (14; 2) 

4l p 
COS Ds 
a Cg 2k 


in place of relations (26), (27) and (28) respectively. 
Relations (35) and (36) reduce to 


5n 
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” Equation (33) reduces to 


2 : 
jy a M NE: ye Shot) 
2p?cyk PU PF). c dkg 


= E Ai) p ; 
EE sles) , nue T 5g T .. (47) 


Taking the inverse Laplace transform, we get 


— _ Hk 
i he jc 2p2CyleyCg $in, t, Rar, er) 
where 
203 
Tm A 
and 
gln, t, Ry, €r) 
2k Tr m ; 3 
ET. al V Enter) et us Is: 72. exp (o z(t i l 
o T e T? "TA! 
Qk 
E 1 €T g^ 4 7" 1 er lee fortc Ki (48) 
2t pi owe N A T es 
and 
n 9 
lf o 2529 - >= (Ruter) œ ERUNT. 
=9)1 v9 S l o VRm tapp . e AERE (Ruter) 
Q 
nk, 1l - Qu A P a) 
x C1Cg +9 t Q 


T 
pA afi 2 Qk 2 
eds 2%. 2. r2exp ) SL p E ll 
e V Ryter'* | I C "i c A ES j 


E 1 


Ee 
Los e % EE jE Jar tor > S EO 


To study the effect of the magnetic field parameter on the induced field 
we can put ey = 0 since er occurs with Ry which is very large in comparison 
with ET: 

Table II gives the values of ¿(7) when Qt = 0:72, er = 0, Ry = 1-9. 
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TABLE II 


Values of (9) in units of 10-12 


5 :6 


1 
= 


ln) 1-043 0-79 0:904 0-911 0-943 0-944. 1-060 


From the graph (Fig. 2) it is easily seen that the effect of Ry, at a certain 
instant of time gradually increases along the length of the solid beginning 
from a particular point of its length. 


109 
1:04 


0:99 


2 
só 
= 


$ (xo > 


—b 
0:5 o4 0:5 OG 
Mi — 
Fia. 2. Effect of magnetic field parameter on induced magnetic fiold 


Taking into 


consideration the approximations for small-time 
(34) reduces to 


equation 


2 
Cc 
taa 


dob. 015 LR GE Ar 4hieg 
TAC) | Vep \" ceu 


1 c, —2R ci) “Bro [Da 
TIRE Ires CO pend CACH H' e c3 2h 3h ^ (50) 
V Tp 4kyezp 


Wh a 1 
The inverse Laplace transform then gives 


Sb: doB 


PAUL t, Ry, er) 


where 
Qk, 
b(n, t, Ray, er) = m 3 e FE 
A nn (2—R xf JE Ly y 
Vi 2Ry ae 
0 
A wj 
8 T? MAE Te 
xe = T Da 
OR A oto ( Fucker 
and 
E 
EE. Ge MEC E UT 
p EEO) °° as 2 


Jiy—l Or (Ry er)! ; Qkin 
C2 2Ry' k? E T 


1 pan a _ Ruter 
+e ES Ry — gà ON aJRy 
Qk Ryte 2 
afi -9t amu nd 28 
.dr 
ders apa 


C1c3 


~2t + m - 
A EST noe. jas i Q9 evi. G- 


O Vh Ru 


Eu+ep Y 


O SR A NER: 
xe | EN aio m . erfe Roe 


to find the effect of magnetic field parameter we pu 
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STABILITY OF CERTAIN TYPES OF FORCE-FREE 
PERIODIC FIELDS 


by M. R. RAGHAVACHAR, Department of Applied M athematics, Indian 
Institute of Science, Bangalore 12 


(Communicated by P. L. Bhatnagar, F.N.I.) 


(Received 21 June 1968; after revision 3 September 1968) 


The stability of a certain type of periodic force-freo fields having neutral 
lines, investigated by Schatzman on tho basis of energy principle, is discussed. 
The normal mode technique is employed restricting to neighbourhood of 
neutral lines. ‘The neutral lines contribute to instability in agreement with 
Schatzman. Repetition of stability analysis employing energy principle 
revealed that every cell appears to be stablo contrary to Schatzman's con- 
clusion that instabilities can be present under certain conditions. Moreover, 
the present study provides the growth rate of the instabilities. 


1. INTRODUCTION 


It was pointed out by List and Schlüter (1954) that cosmic magnetic 
fields occurring in regions of low density should be force-free. As the ionized 
matter is highly conducting large currents can flow giving rise to large Lorentz 
force predominating over other forces such as arising from pressure gradient 
and gravitation or inertia. Consequently the Lorentz force should vanish to 
maintain steady states over considerably long time scales. Such types of fields 
have been characterized by Chandrasekhar (1956), Bhatnagar (1957) and Menzel 
(1964). In particular, Menzel has given a class of periodie force-free fields 
which could support quiescent prominences. Stability of the force-free fields 
has been investigated by Trehan (1957), Woltjer (1958) and Chakraborty 
and Bhatnagar (1960). Schatzman (1960) has discussed the stability of the 
periodic force-free fields of the type 


H = n|- : cos «x sin ly, = sin Ke Cos ly, COS KX COS w| .. (11) 


Menzel, where «2-L/? = o3, « being the constant of proportion- 

e magnetic field. We note in particular that these fields have 

neutral lines. He has employed the energy principle given by Bernstein ef al. 
(1958) and obtained cerbain necessary conditions for instability. 

In the present study we discuss the stability of an ideally conducting 

homogeneous medium in presence of fields of the type (1.1) by normal E 

technique, restricting ourselves to neighbourhood of neutral lines, as the 
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general discussion is mathematically very complicated. For the above fields 
the neutral lines are given by 


z = (2m+1), y=5 "+1; m;n =0, +1, £2,... .. (1.2) 


We find that neutral lines contribute to instability which agrees with the 
result obtained by Schatzman. But the conclusion drawn by him that there 
«2 B? sin? 2w 
32m  YPo 
be correct as can be verified by repeating the Schatzman stability analysis on 
the basis of energy principle. Contrary to Schatzman's conclusion we find 
that in this case there is always stability. 


does not seem to 


can be instabilities present along z-axis if 7? < 


2. LINEARIZED HYDROMAGNETIC EQUATIONS 


We shall work through the physical quantities rendered dimensionless 
according to the following scheme: 


(2, y”, 2) —a(z, y, 2), t = («V4y 


se y FF UT i 
Sp’ , Uu , 
la Po YPo a Va? a Ho 


B? 
where V, = = am is the Alfven speed, po, pọ being the pressure and 


density of the medium in equilibrium. 
The linearized hydromagnetic equations in dimensionless form are: 


0 Wut 

ay e) +div v = 0 ob T zi s 55 (alb) 
PE oo og SUA uus is 4) 
2v > > > > 

a^ —8? v (8p) + Ho x -]- (curl A) x Hy ‘ye oo (83) 
ah s S 

m curl (vx Ho) Ji 2s ne 2 299 (2/4) 


2 
where 82 = (e/V? 2» H, is dimensionless steady state magnetic field, and 
Cs = (YPo/po)* is the sound velocity in the undisturbed medium. 


3. DISPERSION RELATION 
Taking perturbed quantities proportional to 
exp (twt-+-ikz) 
w and E being dimensionless frequency and wave number respectively eqns, 
(2.1)-(2.4) reduce to after eliminating h, 9p and 8p: 
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(S24 H+ Hus H? hue? (He A Hy od Jue [on esi ee 
9H, aH 
EE e ty ay + [at —Hy +H, PT 9Hy Lat rr. ex dug Eie 
oy a Ox Ox 


02H OH, 07H. 
Hy 0x aa tA (s dy ae) le- HoH oss + (824 Hi) U2y — HH yr "yy 
9H, 0H, 
— tt, H,— Hy +H, eat- HyH,—ikH,H;--H; E us, 
OHy| .,, 0H. 02H 0H, H 
H e qn (i km axi Ha E (i 2H: ER) 
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E e 2 9H; 
—HzH guzz— HyH ote + ikS2 — H,Hy—Hy => yum KH, —2H; 2; 
0H, 0H; OH,  . yy OHy 
—ikH¿HByuy+ [za —Hy (1:2. an m) oy +ikHy ES 
—ikH; (ira. A —H, a u— HyH ¿0yy — HzH vey + |a- iH, 
=) ee Va + | 1824 HH, —9Hy He Ha 2H opel z | ikHz:— pus ea) 
oy oy Ox “Oy 
E Hz . : OH. oH 0277, 2 
TH, = —H, ~~ ay —ikHy (un. ^" iH Tr Hs oH 
4 2H gH ytozy + Hwy |u, ES 24H, p | agi |n, A epee 22 vm 
RH HAE (8.8) 
where Hz, Hy, H; are given by 
Hz = —b cos ax sin by, Hy =a sin ax cos by, Hz = cos ax cos by (3.4) 


with a = xfa, b = lja and a2--02 = 1. 
As we shall limit ourselves to the neighbourhood of neutral lines, we put 


m 
curo y= tu a 2 .. (8.5) 
where é and 7 are small. We then have 
HM Sip hg. 
Hy = abn +S Ve Bo i. hd .. (3.6) 
a3b 


Hz = ab£n— eL eno eee 


Substituting for Hz, Hy and H, in (3.1), (3.2), (3.3) and retaining up to 
quadratic terms in £ and y we get 


(054m?) tig Py + tty + (2 — ibn ut Emus, 
+O Vey + Ena +2nve tiknu- (iy — 2Ey-+ikO?)w ; 
eaten] = 0 IQ V. 7o 
Enget O ue Hénin H 2u — ike? 
térve t (CPE?) nn + Reve + (W242ikén)v 
+lilén—E2)we + (10 ¿+ 26 4-i£)u,, = s ca G9 
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ik(O? +12) ug+ihenuy +Hiktu+ikenoe 
+ (œ+ £2)yikvs +ileno+ ¿Eu —2WENWe, 


py Hwet ent (W2—k O — kE —kn)w=0 (3.9) 
where : 


Consistent with our approximation we assume 
u = a9-- 3€ -4-099--a3£?-E a 489 Han? . . . | 
v = bo b,£-E 5n +38? +5489 + bsp? .. - 
w = €o4-C1€ 4-02 4- 03£? E e4£m 052? - . . . 
Substituting the above expressions for u, v and w in eqns. (3.7)-(3.9) and 
equating to zero the coefficients of érg5, 0 < r-+s < 2, we obtain the following 
linear homogeneous equations determining a*s, b's and c's: 


.. (8.10) 


2C%a3+ Wa +O Dd HC — 0 .. .. s. (QUI) 
Wart ikog m0  .. A IA Sa) 
(W2--2)a 4-204 +ikC%, = 0 ee p. 3 .. (3.13) 
Waz — 0 A se css (8g 
(W?--2)a4 —2ikag 4-605 4-25; —2c; +ikcz = 0 e .. (3.15) 
2ag--(W2-E6)ag--9b,--ikbg--ike, +e — 0 .. .. «. (3.16) 
Cia, +20 bs + W?bo--ikO)es — 0 — .. NE Xe. .. (3.17) 
2a 4- (W2-4-2)0, 4-ikC704 = 0 E E: a .. (3.18) 
Webg4-2ikO]eg =0 .. .. ets AO) 
2a, —ikag4- (W2+6)b3+2b¿=0,+ikc> = 0 ac 20) 
2g + Bag (W2-+2)b4-+ 2ikby Hike, 4-26 = 0 SAN 
W2b, — 0 lo n 
ikO2ay 4-ikC?bs - (W2— k267)o, = 0 EE S e EY 
SikC as -ikag--ikO?b,-E(W?—X202--Y)ey —2 0 .. — .. (8.94) 
ikO a, +2ikO bs +ikbo+(W2—k°0?+1)co=0 — .. — .. (8.25) 
¿hay +Hikby+(W2—10744)og=Hk09=0 0. .. ~- (3.26) 
2ikas--24kb,-F(W*—2209)6, — 0... Ea) 
ika,-Rikbs--(Wi2— 420?--4)eg— ky =0 — .. — .. (3.28) 


Equations (3.14) and (3.22) give ag = 0 and bs — 0. We are left with 16 
unknowns and 16 equations. These form three mutually exclusive sets, say, 
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A, Band C. Theset A consisting of the three equations (3.13), (3.18), (3.27) in 
the unknowns a5, bı and c4; the set B consisting of the five equations (3.12), 
(3.19), (3.23), (3.26), (3.28) in the unknowns dy, 05, Co, c; and cg, and the set C 
consisting the eight equations (3.11), (3.15), (3.16), (3.17), (3.20), (3.21), (3.24), 
(3.25) in the unknowns ao, «4, ds, Do, Da, by, c; and cs. 

Each of the three sets leads to dispersion relation when we apply the 
condition for the non-trivial solution of the equations. 

The set A gives the dispersion relation 

W2 = k?0?—4 
or in physical parameters 
A E (8.20) 


T2 — 5— gaps 


This relation represents the sonic modes modified by the presence of inhomo- 
geneous magnetic field. ‘These modes are stable or unstable according as 


le o 
S a Cs 


The set B breaks up into the following two relations: 
(i) W-k0?44=0 as E. .. (3.30) 
(i) W*-2(9—420) W?--R401 —5 0 .. — .. — .. (3.81) 


The first relation is the same as obtained in set A, while the second gives on 
solving 


/ 9212 272 
w? = (ie! — a Aje” mE POHAR .. o (8,82) 


72 
x22 Y 
A y B G_o . 
yz 2» the two values of o? are real and negativo, giving rise 
c 
s 


Now, if £2 < 


2 
: : A. ele V. 

to growing and decaying modes while if 52 > 2 —#, the values of o? are 
nd C» 


s 


complex conjugate, again leading to instability. 
Finally, the set C gives the dispersion relation which breaks up into the 
following two relations: 


(WADWA (1L— 1205) W?-- 1207 ]- W2W241) 4400 (1243) =0 (3:33) 


We can easily show that each of these factors, cubic on JP2 possesses at least 
one negative root < — 3, thus giving growing and decaying modes. 

From the above discussion, we conclude that there is always instability 
in the neighbourhood of neutral lines. This result is in agreement with the 
conclusion drawn by Schatzman that neutral lines contribute to instability, 
However, by employing normal mode technique, we are able to obtain magni- 
tudes of growth rates of instability present around neutral lines. 
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4. Srupy or STABILITY THROUGH ENERGY METHOD 


In this section we again investigate the stability of the system (1.1) on 
the basis of energy principle given by Bernstein ef al. (1958) and also em- 
ployed by Schatzman. 

'The energy principle states that if 


ow = —4 | E. Fé) ar EXER C ay 


where E denotes an arbitrary displacement from the equilibrium, 
HE) = yp grad (aive Jn (a à) x Q— B, X curl Q, 


> > \ 
Q = curl (Ex Bo) 
and the integration is performed over the smallest cell bounded by the neutral 
lines, we have stability or instability according as ôW > or < 0. 
Following Schatzman, we assume 


Se 2 X (Gmn, bmn, mn) exp ince nly + rz) nS .. (4.2) 


m n 


The integral (4.1) reduces to 


E S » » [5s : T (nn) enn : F(Enn)| a 0 (4.3) 
mo n 


where K denotes a positive constant. 


If R = SK Bo the contribution to (W) from magnetic part is 


] (=> c > > =. > 
2UR* o li curl R—curl curl Raz. le curl Z*— curl curl r2 .. (4.4) 
Taking Rmn = (Amn, Ban, Cmn) exp i(m«a--nly-+-rz) e (4.5) 


the m, nth component of (4.4) reduces to, after dropping subscripts from 
A, B, C, A*, B* and C*, 
mnkl( A.b* + A* B) --nlr(BC*-4-B*C)--mkr(AC*-4- A*Q) 
— (n2 4-72) 4LA* — (mr? 4-72) B.B* — (m? --9202)0 0 
-Eai[r4B* — A*B)--m(BC*— B*O)g-nI(CA*— C*4)] ..— (4.6)] 
We ean express the components of Rin terms of the components of E to 


obtain A, B, C as follows: 


Ho 
Amn = a (bm-1, n-1-F Um 3, n+1F0m+1, n-1-F Uma, n+1) 


Ko, ; 
XE (Cm—1, n-17FCm-1, n1— m+, n-1—Cm34, «| e (4.7) 


T This expression differs from that of Schatzman. Perhaps some error has crept in his 
work due to some slip. 
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Ho 


l 
Brn = — E (Cm—1, n-1—Cm-1, 1+1Cm41, 1-17 m1, n+1) 


4 Lo 
+(4m-1, n-1- 4m-1, nti FH Um41, n-1 + Um+1, z (4.8) 
Ho 
Onn = Hm [&(5 1, n-1+4%m-1, n4-1— 6m41, n-1— 6m41, n41) 
jl Q 
by a, n-1—bm-1, n+1+9m+1, n-1—0m+1, n41)] (4.9) 


Substituting for A, B,...in terms of a, b, ... and collecting the terms with 
index mn, we finally obtain the following expression for (BW): 
E 6; 2. en) lal? 
—=3 (8) mn = Sa SoA ee aU nx?) |a |? 
A 


KH, 


2 2/2 
+ 5 (r+ mr?) [6]? + lemen SE fet 
2 2 
+ [rxa-pnlb]245| roma] 24 E |ro-4+n0b 2 .. (4.10) 


where Cs, V are the sound speed and Alfven speed respectively. As the 
expression (4.10) shows (8W) is positive definite and, therefore, the system as a 
whole is stable. This is in contradiction to the result obtained by Schatzman 
where he has shown that instabilities are present in z-direction if 


«2B sin? 2w 


¡PEA 
32m ypo ' 


w —ian-l-, 
K 


5. CONOLUSION 


By the energy principle which in a way provides a gross treatment, every 
cell appears to be stable but the contribution to èW by the neighbourhood of a 
neutral line is negative as proved by Schatzman or by the present investigation 
in $3. This explains the motivation of the present investigation. Moreover, 
the present study provides the growth rate of the instabilities. 
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